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SECTION  I 


INTRODUCTION 

Widening  flight  envelopes  coupled  with  increasingly  stringent  opera- 
tional demands  on  military  aircraft  have  forced  the  use  of  more  complex 
stability  and  control  augmentation  systems.  It  is  recognized  (Refs.  1 
and  2)  that  this  trend  has  created  problems  in  interpreting  the  applicable 
military  specification  for  the  flying  qualities  of  piloted  aircraft,  MIL- 
F-8785B(ASG)  (Ref.  3).  For  example.  Section  3-2. 2.1,  which  relates  to  the 
short-period  response,  makes  no  provision  for  the  dynamics  of  the  aircraft 
plus  flight  control  system  (FCS)  if  these  are  different  in  dynamic  form 
from  those  for  an  unaugmented  aircraft.  The  specification  assumes,  in 
essence,  that  the  longitudinal  dynamic  response  will  have  only  two  modes 
(phugoid  and  short  period)  in  the  low  and  mid  frequency  regions,  when  in 
practice  there  may  be  additional  significant  modes  in  these  frequency 
regions  due  to  equalization  of  the  FCS  feedbacks . Because  of  these  addi- 
tional modes,  an  augmented  aircraft  may  appear  quite  different  to  the  pilot 
than  an  unaugmented  aircraft  which  appears  to  have  the  same  short-period 
dynamics  (Ref.  U),  That  is,  it  may  develop  that  the  short-period  mode 
characteristics  meet  the  specification  and  yet  the  aircraft  actually 
responds  as  if  it  had  short-period  characteristics  which  do  not  meet  the 
specification.  This  opens  the  possibility  that  it  will  receive  low  ratings 
from  evaluation  pilots. 

How  should  this  situation  be  treated?  Are  new  requirements  needed 
which  account  directly  fer  the  additional  dynamics,  or  is  it  possible  to 
develop  lower-order  mathematical  models  of  the  higher-order  system  that 
can  be  reliably  compared  against  the  MIL-F-8785  requirements  as  presently 
stated?  Moreover,  the  criteria  of  MIL-F-8785  presume  an  analog  augmenta- 
tion system.  There  is  now  an  additional  concern  for  response  characteris- 
tics which  are  unique  to  use  of  digital  augmentation  systems.  It  is  the 
purpose  of  this  study  to  identify  those  characteristics  which  are  unique 
to  the  digitally  controlled  system,  and  to  review  the  quantitative  tools 


available  which  permit  an  assessment  of  t tie  el't'eets  of  these  unique  eh arse  - 
terlst  los  relative  to  the  Ml  h-F-t'/'b'  requirements.  The  larger  Issue  of 
what  additional  flying  qualities  criteria  are  required  Is  outside  the  scope 
of  this  research. 

The  results  of  this  study  Indicate  that,  two  import. ant  characterist  ics 
arc  Introduced  by  digital  control  laws.  They  arc: 

• The  effective  delay  Introduced  hy  the  A/D  and  D/A 
process  and  the  delay  introduced  by  the  digital 
algorlttuns  and  computational  frame  time  of  the 
computer  cycle. 

• The  "control  roughness"  or  inter-sample  ripple 
Int  roduced  when  t.he  digital  computer  Is  coupled 
to  t.he  control  actuators  using  data  holds. 

The  first  eharaeterlst !c  is  of  great  concern  to  the  flying  qualities  com- 
munity since  even  relatively  small  delays  are  potentially  important  in  cer- 
tain closed-loop  piloting  tasks  involving  motion  cues  and  p l lot -induced 
oscillation  possibilities.  Measures  of  the  effective  time  delay  Introduced 
by  the  a/D,  computation.  D/A  conversion  process  can  be  computed  using  a 
vartety  of  analytical  techniques  that  apply  at  the  sampling  Instants.  For 
example,  w’-,  w-,  or  ..-domain  analyses  and  discrete  frequency  response  tech- 
niques can  he  applied  effectively  to  obtain  both  a quantitative  and  qualita- 
tive evaluation  of  the  delay  Introduced. 

The  second  characteristic  Is  of  concern  since  aircraft  response 
result  ing  flaw  control  roughness  Is  In  effect  an  .additional  disturbance 
source.  Moreover,  t.he  step-like  control  deflect. tons  caused  by  the  use  of 
.•.ero-order  data  holds  c.an  produce  excessive  actuator  wear  and  fatigue 
damage,  reducing  service  1 l fe  and/or  reliability  of  the  aircraft,  system. 

Tills  second  characteristic  Is  more  difficult,  to  assess  since  It.  involves 
the  response  of  t.he  cont  inuous  system  during  the  inter-sample  interval. 
Therefore,  analysts  tools  such  as  the  "sampled  spectrum"  (sampled  frequency 
response)  arc  of  little  value.  It.  Is  this  observation  which  prompted  the 
development. , during  this  study , of  a basic  analysis  tool  for  computing  the 
spectral  content  of  the  continuous  responses  of  a discretely  controlled 
system.  We  are  tempted  to  describe  this  method  of  analysis  as  "frequency 


L 


response"  evaluation  since  it  applies  equally  for  continuously  controlled 
systems  as  well  as  multi-rate  discretely  controlled  systems.  However, 
this  must  be  further  qualified.  This  is  necessary  because  of  possible 
confusion  with  the  concept  of  the  "sampled  spectrum."  In  computing  the 
sampled  spectrum  one  finds  the  lowest  frequency  sine  wave  that  fits  the 
sampled  response  at  the  sampling  instants.  The  sampled  spectrum  is  com- 
monly, but  incorrectly,  understood  to  be  the  frequency  response  of  a dis- 
cretely controlled  system.  For  this  reason  it  is  preferable  to  describe 
the  new  method  more  precisely  as  the  "continuous  frequency  response  of  a 
discretely  excited  system."  Once  a method  is  in  hand  for  computing  the 
spectral  content  of  the  continuous  response,  then  one  has  the  means  for 
assessing  both  of  the  digital  characteristics  discussed  above.  The  fre- 
quency response  magnitude  data  can  be  used  to  quantify  control  roughness, 
while  the  phase  data  can  be  used  to  quantify  the  effective  time  delay. 

There  is  another  main  objective  of  this  study  effort  — that  of  giving 
the  practicing  engineer  working  familiarity  with  three  analytical  tools  we 
consider  to  be  well  suited  for  the  analysis  of  digitally  controlled  systems 

1)  Analysis  (and  synthesis)  in  the  w'-domain. 

2)  Multi-rate  transform  domain  approaches. 

3)  The  continuous  frequency  response  of  a discretely 
excited  system. 

The  w’-domain  is  related  to  the  well-known  w-domain  by  a scalar  trans- 
fbrmation  and  to  the  z-domain  by  a bilinear  algebraic  transformation.  We 
believe  that  those  engineers  skilled  in  frequency  domain  design  procedures 
will  immediately  feel  at  home  with  analysis  in  the  w'-domain,  since  all 
analog  control  system  design  technology  transfers  completely  for  digital 
control  system  design. 

The  second  item  deals  with  a multi-rate  transform  domain  approach  that 
has  been  developed  into  an  effective  tool  for  analyzing  the  transient  inter 
sample  response  of  discretely  excited  systems.  That  is,  it  yields  recur- 
sion equations  describing  the  inter-sample  performance  to  any  degree  of 
fineness  desired  without  increasing  computer  storage  requirements.  A basic 
understanding  of  this  approach  leads  to  the  development  of  the  third  tool 


5 


which  gives  the  continuous  spectrum  of  a discretely  controlled  system. 

It  is  this  third  tool  which  we  perceive  as  being  of  the  greatest  value  to 
those  engineers  working  at  the  flying  qualities,  digital  flight  control 
system  interface. 

In  Section  II,  we  will  first  review  the  fundamentals  of  sampled  data 
control  theory  to  the  extent  necessary  to  "refamiliarize"  the  engineer 
with  its  terminology  and  background  mathematics.  No  proofs  will  be  given; 
we  will  simply  state  results  and  give  illustrative  examples  to  demonstrate 
viewpoints  which  will  be  needed  at  a later  point  in  order  to  develop  an 
understanding  of  the  three  analytical  tools  mentioned  above.  This  section 
will  also  serve  to  introduce  a multi-rate  terminology  and  give  the  reader 
n chance  to  gain  familiarity  with  the  notation  (always  a crucial  factor  in 
how  easily  a report  re  is ) . 

Section  III  is  devoted  to  modest  extensions  of  classical  analysis  and 
synthesis  techniques  for  digitally  controlled  systems.  Emphasis  is  upon 
analyses  conducted  in  the  w' -domain.  The  w’ -domain  offers  the  advantage 
that  non-minimum  phase  effects  of  the  sampling  and  data-hold  operations 
and  of  sampling  rate  can  be  directly  accounted  for  without  approximation 
while  using  conventional  frequency  domain  design  tools  such  as  root  locus 
and  Bode  plots.  These  conventional  frequency  domain  design  tools  can  be 
used  to  considerably  greater  advantage  in  the  w'-domain  than  in  the  w-  and 
z-domain,  because  several  more  powerful  analogies  between  the  '-domain  and 
tire  w’-domain  exist.  These  analogies  are,  in  a sense,  the  key  to  exploit- 
ing the  w’-domain  for  design  purposes,  making  direct  design  in  the  w'-domain 
more  attractive  than  either  optimal  procedures  or  the  numerous  approximate 
methods.  First,  the  basic  properties  which  make  the  w'-domain  preferable 
to  the  z-doma.in  or  w-domain  are  reviewed.  Following  this,  illustrative 
examples  are  used  to  highlight  the  analogies  between  a and  w'.  Finally, 
we  demonstrate  the  clear  manner  in  which  tile  noiuninimum-phnse  effects  of 
sampling,  data-hold  operations,  and  sample  rate  are  evident  in  the  w'-domain 
analyses . 

Almost  all  existing  or  contemplated  digital  flight  control  systems  are 
multi-rate  in  nature.  The  predominant  reason  for  tills  is  the  necessity  to 
work  within  the  basic  restrictions  imposed  by  computer  word  length,  memory 


capacity,  frame  time,  etc.  Thus,  it  is  typical,  for  example,  to  update 
inner  loops  at  a faster  rate  than  the  "slower"  outer  loops. 

A basic  transform  domain  approach,  the  "T/N  method,"  which  eliminates 
a variety  of  dimensional  and  indexing  problems  associated  with  state  tran- 
sition (Ref.  5)  and  switch  decomposition  (e.g.,  see  Refs.  6-8)  methods,  is 
developed  in  Section  IV.  The  approach  is  very  efficient  for  computing  the 
response  of  discretely  controlled  continuous  systems  at  both  sampling 
instants  and  at  equally  spaced  times  during  the  inter-sample  interval.  The 
methods  also  requires  a minimum  of  information  to  define  the  system  for  the 
computer  program.* 

The  T/N  method  discussed  in  Section  IV  is  used  as  a departure  point  in 
Section  V to  extend  frequency  response  concepts  to  include  the  continuous 
frequency  response  of  a discretely  controlled  system.  The  discrete  fre- 
quency response  concept  (e.g..  Ref.  7)  has  hot  been  particularly  productive 
for  the  analysis  of  discretely  controlled  systems,  since  it  is  limited  to 
determining  the  amplitude  and  phase  of  the  single  sinusoid  that  fits  the 
output  samples  of  a single-rate  system  at  the  sampling  instants.  Develop- 
ment proceeds  by  first  removing  this  restriction  for  open-loop  systems  and 
then  extending  the  results  to  single-rate  closed-loop  systems.  Finally,  the 
solution  for  the  multi-rate  closed-loop  case  is  developed  using  the  results 
of  Section  IV. 

Although  each  section  contains  a representative  number  of  illustrative 
examples,  an  application  pertinent  to  the  flying  qualities/flight  control 
system  interface  is  not  taken  up  until  Section  VI. 

Section  VII  contains  the  conclusions  and  recommendations.  The  appendices 
contain  support  material,  much  of  which  will  be  familiar  to  practicing  engi- 
neers. 


•But  the  computer  program  is  somewhat  restricted  as  to  the  form  of  the 
problem  that  can  be  accommodated. 
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SECTION  II 


i 

i 


T 


I 


MATHEMATICAL  PRELIMINARIES 

A.  PROBLEM  SCOPS 

Systems  governed  by  linear  differential  and  difference  equations  with 
constant  coefficients  are  reviewed.  This  is  a framework  which  can  ade- 
quately deal  with  aircraft  nonlinear  equations  linearised  about  an  opera- 
ting point.  In  the  treatment  of  discretely  controlled  systems,  "impulsive" 
sampling  operations  are  assumed.  All  developments  will  use  vector  nota- 
tion, although  scalar  examples  will  be  used  frequently  for  illustrative  pur- 
poses. In  the  interests  of  brevity  we  will,  wherever  possible,  pursue  the 
direct  and  remarkably  effective  methods  of  Aseltine  (Ref.  9)  in  bringing  the 
basic  issues  to  the  fore, 

B.  • -DOMAIN  PRELIMINARIES 

The  partial  fraction  expansion  table-lookup  approach  for  transferring 
information  from  one  domain  to  another  will  be  used  frequently  in  subse- 
quent mathematical  developments.  Although  this  is  a familiar  topic,  our 
viewpoint  and  approach  is  sufficiently  different  to  warrant  a careful 
review. 

A main  objective  of  this  section  is  to  produce  a table  of  abbreviated 
transforms  that  will  permit  one  to  transfer  between  the  t,  3,  z,  and  w’ 
domains  in  a relatively  easy  and  straightforward  manner. 

We  start  by  developing  a transform  table  (from  the  time  domain  to  the 
s-domain)  using  the  Laplace  transform  of  an  exponential  time  function  and 
several  properties  of  the  Laplace  transform.  We  will  then  mimic  this  pro- 
cedure in  order  to  augment  the  table  with  z-trans form  and  time-advanced 
z- transform  entries. 


P 


* "‘V 


To  begin,  compute  rC[e~Ht]: 


£ [e""1 ) 


dt 


e-(s+n)t 

» n 


e-(s+a)t 
— (s  + a) 


to 


O 


1 

s + a 


It  is  convenient  to  use  the  notation 


-at 

e 


I 

s + a 


0) 


to  indicate  that  e corresponds  to  1 / ( s + a)  and  vice  versa.  This  trans- 
form pair,  coupled  with  two  other  properties  of  Laplace  transforms,  provides 
the  basic  information  needed  to  develop  a table  of  transform  pairs.  The 
first  (trivial)  property  is: 


f,(t)  + kf2(t)  F-| (s ) + kF?(s) 


(2) 


The  second  property  of  interest  is: 

tf(t)  <=*  - ~ F(s)  (*) 

Note  further  that  "a"  in  Eq.  1 is  a parameter  which  is  not  restricted  to 
have  real  values.  For  example,  let 

a -*■  a + jb  (U) 


in  Eq.  1 and  obtain 


1 

(s  + a)  + jb 


(?) 


p-(a+jb)t  = e“at,fcos  bt  - j sin  bt] 


Rationalize  the  right-hand  side  of  Eq.  5 to  obtain: 


—at  , . . —at  . , , 

e cos  bt  — je  sin  bt 


,s  + a)  - ,1b 


(s  + a)2  + \>c 


Using  the  property  of  Eq.  2 now  gives  two  additional  transform  pairs  (one 
for  the  damped  cosine  wave,  the  other  a damped  sine  wave): 


-at  , . 
e cos  bt 


-at  • 

e sin  bt 


(s  + a)2  + b2 


(s  + a)2  + b2 


This  approach  for  developing  the  table  leads  to  a frequency  domain  descrip- 
tion of  a damped  sine  or  cosine  wave  which  describes  the  poles  in  terms  of 
the  wave  damped  natural  frequency  and  an  exponential  damping  envelope  rather 
than  damping  ratio  and  undamped  natural  frequency.  That  is,  the  descrip- 
tion is  in  terms  of  [(s  + a)2  + b2]  rather  than  [s2  + Ft^s  + o>q].  For  the 
purposes  of  this  report,  the  form  given  in  Eqs.  7 and  8 will  prove  more 
useful  than  the  alternative  description. 

We  may  set  a = 0 in  Eqs.  7 and  8 to  obtain  the  transform  of  sine  and 
cosine  waves. 


cos  bt  <=-> 


s2  + b2 


sin  bt  <=>  — r- 

s'  + b' 


To  continue,  setting  a = 0 in  Eq.  1 produces  the  transform  pair  for  the 
unit  step  function. 


u(t)  <=i>  — 
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Applying  the  property  given  in  Eq.  3 to  Eq.  11  yields  the  transform  pair 
for  a rump  Amotion: 


and  obtain  the  transform  pair 


Repeat  the  application,  this  time  to  Eq.  1 
for  t--': 


In  a like  manner,  the  application  of  this  property  to  Eq.  1 gives 


By  proceeding  in  this  fashion  we  have  now  built  up  Table  1 . This  abbre- 
viated table  is  sufficient  for  most  of  our  needs  since  the  occurrence  of 
multiple  poles  in  aircraft,  stability  and  control  work  is  somewhat,  rare. 
Later  we  will  use  the  same  approach  to  develop  a table  of  z-trans forms . 
This  will  result  in  a third  column  for  Table  1. 


Consider  next  an  efficient  method  for  expanding  a given  rational  poly- 
nomial function  F(s)  in  terms  of  the  entries  of  Table  1.  Since  the  entries 
of  Table  1 constitute  the  set  of  elementary  time  functions  of  interest.,  it 
makes  good  sense  to  expand  F(s)  in  terms  of  these  time  functions.  To  ill  us 
t.rat.e,  given  the  proper  rational  F(s), 


then  it.  is  convenient  to  expand  the  right-hand  side  in  terms  of  partial 
fractions,  viz.. 


* 


1*'S*  » S 4 


(s)(s  f l)(s-  4 !.s  t I'M 


P(«>  (n,) 

* « + i (s  4 ) ■ * (■■  y. 


s + ts  > r*  - (s  + p)*-  + (^)i;  (17 

Otu'o  the  partial  traction  coefficients  A,  D , C.  and  D are  known,  the  time 
function  can  be  written  immediately  as: 

f(t)  Au( t ) + Be~  4 e-'^  ( C cos  H 4 p sin  M. ) (18' 

There  is  a variety  of  ways  to  find  A,  B,  0,  and  P.  For  example,  multiply 
each  side  of  tv}.  In  by  s and  lot  s - 0 to  find  A 1.  Next,  multiply  eacli 
saie  by  (s  * )\  and  let  s -1  to  find  P I.  Next,  multiply  each  side  of 
t lie  equation  by  s and  let  s «.  to  find  out  that  A + B 4 C = % and  there- 
fore 0 1 • Finally,  pick  any  convenient  value  of  s and  evaluate  Eq.  1o 

for  the  remaining  unknown,  P (D  l also).  This  approach  is  neat  but  not 
very  general.  We  need  a method  better  suited  for  machine  computation.  Such 
a method  is  described  next. 

Let  F(  s i be  a proper  rational  function  of  the  form 


with  a simple  pole,  (s  t a).  Obviously, 


4 a)N(s) 
D(s) 


However,  if  (s  + a)  is  not  explicitly  factored  out  of  the  denominator,  an 
indeterminate  form  will  result,  when  s is  sot  equal  to  -a  (0/0).  To  resolve 


this  difficulty,  without  explicitly  dividing  out  (s  + a),  use  L'Hopital's 
rule: 


A _ (s  + a )N ' + N1 

= — d 7F5 1 


where  D - dD/ds.  Equation  Z1  treats  all  simple  real  poles  and,  since  only 
polynomials  are  involved,  is  ideally  suited  to  machine  computation  using 
"nesting  procedures."  Moreover,  the  same  polynomials  are  used  for  all  par- 
tial fraction  coefficient  evaluations  involving  simple  roots. 


Again,  an  indeterminate  form  results  if  [(s  + a)2  + b2]  is  not  explicitly 
divided  out.  L'Hopital’s  rule  permits  one  to  avoid  the  division  step: 


m 


. Lll  WHi 


- "WPII  w 


This  simplifies  to 


B — jA 


SN 


s = -a+jb 


(26) 


and  clearly 


B = 2 Re 


N 


s = -a+jb 


(27) 


A = — 2 Im 


N 


s = -a+jb 


(28) 


Other  expressions  can  be  obtained  for  multiple  roots.  For  example,  if: 


/ v w 

>(s)  = — 


(s  + a)^  (s  + a) 


(29) 


One  can  show  (see  Appendix  D): 


A - ^ 
A " D"  i 


B = 


s = -a 


2N’ 

D" 


AD' 

W 


(?0) 


s = -a 


Some  illustrative  examples  demonstrate  the  simplicity  of  the  process. 
First,  let: 


F(s)  = 


3s  + 2s  + 
s(s  + 1 


^sr’  + 2s  + 1 
s'  + 2s2  + s 


(3D 


Therefore , 


3 s + 2s  + 1 
s + 2s^  + s 


(s  + 1 (s  + 1 ) 


_C_ 

S 


(?2) 


Since  D’  - 3s"1  + Us  + 1 and  D"  = 6s  + U,  D'"  = 6. 


+ 2s  + 1 


s = -1 


2(6s  + 2) 
6s  + 4 


(-2)6 

.1  (3)(6s  +4)  s =_1 


As  a second  example. 


F(s)  = 


s'1  — 5 

s3  + 3s2  + 7s  + 5 


(s  + 1)(s2  + 2s+5)  (s+1)[(s+1)2  + ('2)2] 


A + B(2)  + C(s  + 1 ) 
s + 1 (s  + I)2  + (2)2 


Since  D’  = 3s2  + 6s  + 7, 


A - i 


s2  — 5 


= -1  3s2  + 6s  + 7 s = _i  3-6+7 


2 Re  7^ — - — -Tv-*- = 2 

(3s  + 6)s  + 7 s=_1+j2 


B = -2  Tm 


'3s  + 6)s  + 7 s = - ! + j 2 


Those  readers  with  pocket  calculators  having  rectangular  to  polar  conver- 
sion features  will  have  little  problem  in  verifying  the  coefficient  values 
for  C and  B. 


Thus  fur  we  have  dealt  only  with  "proper"  rational  functions  wherein 
the  order  of  the  numerator  polynomials  is  at  least  one  order  less  than  the 
denominator  polynomials.  In  the  event  F(s)  is  not  proper,  it  can  be  con- 
verted (using  synthetic  division)  into  two  parts  — the  first  is  a poly- 
nomial in  s,  the  second  part  a proper  rational  function.  For  example. 


F(s) 


s'1  + + bs‘  + 5>s  + !( 


S'  + 


+ 


p 

+ s + 


s -t  2 

(s  + l)-'  4 (l)^ 


(4l) 


There  fore , 

f(t)  - b"(t)  + b’(t)  + &(t)  + e_L[sin  t < cos  tj 


(4f) 


We  assume  the  reader  is  familiar  with  Laplace  transform  of  impulse 
to  the  extent  that  1 is  recognized  as  an  impulse  in  time,  &(t);  s 
nized  as  a "doublet,"  B’(t);  and  sy  is  recognized  as  a "triplet," 


functions 
is  recog- 
5" ( t) . 


Before  proceeding  to  the  development  of  a table  of  z-transiorms  (in  the 
next  section),  we  pause  to  review  the  concept  of  the  system  transfer  func- 
tion. Given 


f(t)  <=>  F(s) 


(4?) 


we  presume  familiarity  with  the  transform  pair  (or  refer  to  Appendix  A) 

<=>  sF(s)  - f(0)  (44) 


which,  through  repeated  application,  enables  us  to  transform  the  differen 
tial  system 


x + 2x  + 2x  = f(t) 


(4f>) 


in  time  to  an  algebraic  system  in  s; 


(4b) 


1 


This  equation  highlights  three  basic  parts  characteristic  of  all  linear 
systems: 


A:  The  part  related  to  the  system 

B.  The  part  related  to  the  response 

C The  part  related  to  the  input  and 

initial  conditions 


Dividing  through  by  Part  A gives 


X(s)  = --KS) + IiL.+.  2Mo)  , _ x(0) 

+ 2s  + 2 s^  + 2s  + 2 s?  + 2s  + 2 


X(s)  = W^sjpfs)  + W2(s)x(0)  + W5(s)i(0)  ( 

Wt,  VS,  and  W3  are  the  "transfer  functions"  (or  system  functions)  which 
relate  the  transforms  of  the  output  components  to  the  transform  of  the 
input  components.  For  example,  let  the  initial  conditions  x(0)  and  x(0) 
be  -ero.  Then  W^(s)  can  be  written  as; 


W, 


(*9) 


It  is  in  this  sense  that  transfer  functions  characterize  the  system  inde- 
pendently of  the  input.  The  Laplace  transform  of  the  system  time  response, 
on  the  other  hand,  can  always  be  expressed  as  a sum  of  products  of  system 
transfer  functions  with  the  Laplace  transforms  of  the  various  corresponding 
input  time  functions  and  initial  conditions. 

C.  z -DOMAIN  PRELIMINARIES 


An  abbreviated  table  of  transforms 
starting  with  an  exponential  transform 
taken  with  z-transform  pairs.  Let  z = 


was  built  up  in  the  previous  section, 
pair.  Here  a similar  approach  is 
e for  the  purpose  of  constructing 


In 


this  table.  Then  the  following  equation  may  be  regarded  as  the  definition 
of  the  z-transform: 


fT(t)  = f(0)B+(t)  + f(T)6+(t-T)  + f(2T)8+(t-2T)  + •••  (50) 

The  Laplace  transform  for  delayed  time  functions  (refer  to  Appendix  A)  can 
be  used  to  transform  Eq.  50  into  the  frequency  domain. 

F^s)  = F*(s)  = £ fCnTje-8"1  = £ f(nT)z“n  = F(z)  (51) 

n=0  n=0 


sT 

where  z = e . 

Some  comments  on  notation  are  in  order.  We  use  f^(t)  <£=>  F^s)  to  denote 
that  the  time  function  f(t)  is  sampled  at  1/T  samples  per  second.  F^(s) 
denotes  the  Laplace  transform  of  the  sampled  continuous  time  function,  given 
that  the  continuous  time  function  has  the  transform  F(s).  While  the  tradi- 
tional notation  for  F^(s),  which  is  F*(s),  is  perhaps  better  established, 
we  have  elected  to  use  the  superscript  notation  because  it  permits  an  expli- 
cit statement  of  the  sampling  rates  involved  in  multi-rate  systems  (the  main 
concern  of  this  report).  This  point  is  discussed  further  in  Appendix  B. 

Let  us  now  proceed  to  develop  an  abbreviated  table  of  z-transforms.  Let 
f(t)  = e-*^ , so  that  Eq.  51  becomes 

fT(t)  <->  F(z)  = f)  e— anT  z-n  (52) 

n=0 


1 - (e~aT/z ) 


<=> 


z - e 


-aT 


(53) 


We  now  Know  that  this  continuous  time  function  sampled  at  1 /T  samples/oecond 
gives  the  transform  pair 


(e^)T 


z — e 


■^aT 


(*>M 


Just  as  in  the  case  of  the  Laplace  transform  we  can  let  "a"  take  on  conveni- 
ent complex  values.  For  example,  let  a a + jb. 


[e 


-(a+jb)t|r 


_ e-nT"e-|bT 


(55) 


or 


(e”flt  cos  bt.  — je”*^  sin  bt)^  <$=$ 


— a7 


z - e (cos  bT  - j sin  bT) 


(*•«’) 


Rationalization  of  the  right-hand  side  of  Eq.  5o  results  in: 


cos  bt  - Jc-“‘  sin  bt)T  ■?(■«  -jTf  c°s  b\l  r feT*  ”, 

(z  - e cos  bT)"  + (e  1 sin  bT)c 


-aT 


Therefore 


/ —at  V1.\T  , „ z(z  — e cos  bT) 

(e  cos  bt)  <=>  -* -nfr • 

z — 2e  cos  bTz  + e 


(57) 


and 


(e~Q,T  sin  bT)T  <=£ 


— aT 

ze  sin  bT 


2 _ ■— aT  , m 

z — 2e  cos  bTz  + e 


We  now  have  the  transform  for  the  sampled  damped  sine  and  cosine  pair.  Con- 
tinuing, let  a = 0 in  Eqs.  57  and  58  to  obtain  the  transform  for  sampled 
sine  wave  and  cosine  waves. 
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m z(z  — cos  bT) 

( COS  bt  ) 

z‘  — P cos  bTz  ■*  I 


M) 


(sin  bt)1  <:-r>  — K 


7.  sin  bT 


(60) 


7.c  — 2 cos  bTz  + 1 


Further,  let  a 0 in  Eq.  51*  to  obtain  the  z- transform  for  a sarapled-unit 
step. 


/e-^^o)  MOlT  <=*  T^T  {61) 

To  complete  the  table,  first  find  the  transform  oi  tt(t): 


[tf(t)f  = [uTf(nT ) 1 <=>  £ (nT)f(nT)z 

nO 


— n 


— zT  £ f (nT ) [-(n)z(  n_l)] 


n 0 


We  recognize  the  term  in  brackets  as  the  derivative  of  z'n  with  respect 
to  z.  Therefore: 


[tf(t)|T  -aT  aj  £ f("T)z  n 

n 0 


(60) 


Equation  6.°  is  completely  analogous  to  the  s-domain  result  (Eq.  31.  With 
it  we  can  complete  our  abbreviated  table  of  z-trans forms , since 


. -zT 


-aT 

z — e 


zTe 


-aT 


(63) 


to 


Setting  a = 0,  the  z-transform  of  a ramp  is  obtained: 


Tz 


(64) 


We  may  now  complete  Table  2,  showing  the  relationship  between  time  and 
the  s-  and  z-domains.  Strictly  speaking  there  should  be  a "f(nT)"  column 
as  well,  but  t.he  manner  in  which  the  table  is  presented  is  sufficient.  It 
forces  us  to  remember  that  f(nT)  can  be  thought  of  as  the  continuous  f(t) 
sampled  at  intervals  of  T seconds  starting  at  t = 0. 

Given  an  F(z),  one  can  expand  in  partial  fractions  in  exactly  the  same 
manner  as  was  done  for  F(s).  We  will  defer  demonstrating  the  process  until 
the  additional  considerations  introduced  by  data  holds  have  been  reviewed. 

D.  FUISE  TRANSFER  FUNCTIONS  AND  DATA  HOLDS 

The  correspondence  between  a given  F(s)  and  the  associated  F(z)  from 
the  transform  table  is  depicted  in  Fig.  1 from  the  transfer  function  point 
of  view.  From  Fig.  1 and  the  results  of  Appendix  B: 


C(z) 


CT(s ) 


z e 


sT 


GT(s) RT(s J 


z e 


sT 


G(z)  R(z) 


m 


Figure  I . Correspondence  Between  Transfer  Function  and 
Pulse  Transfer  Function 
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AN  ABBREVIATED  TABLE,  F(z)  ADDED 


F(Z)  = 2[^(t)] 


f(t) 

F(a)  £[f(t)] 

u(t) 

s 

t 4, 

s* 

f 

a* 

,^'t 

1 

s * a 

l 

-Ht 

te 

1 

(s  + a)" 

sin  bt 

b 

s*  -f  b* 

That  is,  given  that  C.(s)  is  isolated  by  samplers  (and  it  is  in  Table  : 
then  G(z ) can  be  written  directly.  However,  if  G(s)  is  the  product  of 
several  of  the  tabulated  functions  of  s,  then  a partial  fraction  expansion 

ot  G(s)  must  be  made  before  G(z)  can  be  found.  For  example,  if  G(s)  -1  s;\ 
then 


G(s) 


( 66 ) 


However,  if  0(s)  l/[s?(s+a)],  then 


G(a) 


_A_  ( C 

s2  s s + a 


(67) 


or 


G(  2) 


8Tz  Bz  Cz 

“ ■ • + 


(z  - O' 


(68) 


Moreover,  the  situation  depicted  in  Fig.  l does  not  represent  the  form  of 
the  problem  which  is  usually  of  interest  since  no  data  hold  is  included. 

The  form  of  interest  is  shown  in  Fig.  2, 

In  iig.  2,  M(s)  represents  the  transfer  function  of  a data  hold.  Using 
the  results  of  Appendix  B,  we  can  write 


C 


T 


(gm)t  rt 


(69) 


Figure  Typical  Data  Hold  Configuration 


The  data  hold  is  a physical  device  that  takes  the  samples  of  the  signal, 

rn 

R , and  constructs  a continuous  signal  which  in  turn  forces  the  continuous 
system  represented  by  G(s).  It  is  convenient  to  think  of  the  data  hold  as 
a coupler  between  the  digital  computer  and  the  control  actuator  input.  This 
coupler  is  usually  inherent  in  the  digital-to-analog  (D/A)  converter  device. 

We  will  assume  that  the  reader  is  familiar  with  the  characteristics  of 
various  data  holds.  (Their  properties  can  be  reviewed  in  Refs.  7 and  8.) 

The  presence  of  the  data  hold  in  a single-rate  system  is  only  a mildly 
complicating  factor.  This  is  especially  so  if  we  are  using  a zero-order 
hold  which  has  the  transfer  function: 


M 


1 — e 


— sT 


(TO) 


For  example,  let  G(s)  = a/(s  +a),  M = ( 1 — e~sT)/s,  then 


(GM)J 


— sT 


0 -e-sT)a]T 


s( s + a) 


(TO 


But,  1 - e s (z  - 1 )/z  and  is  unaffected  by  the  sampling  operator.  There- 
fore, 


(GM)1 


- 1 


z — 


1 — e 


K-H’-W 


-aT 


U - e-*T) 


(7?) 


As  we  shall  see,  this  simple  computation  becomes  more  complicated  when 
the  data  hold  and  input  sampler  are  working  at  a different  sampling  period 
than  is  the  output  sampler.  This  situation  will  be  encountered  later  in 
the  analyses  of  multi-rate  sampled  systems.  For  single-rate  sampled  systems 


t 


it  cun  be  appreciated  that  the  partial  fraction  technique  is  quite  ade- 
quate for  computing  overall  pulse  transfer  functions  (from  the  output  of 


a sampler  to  the  output  of  the  previous  one).  This  result  pertains  whether 
or  not  there  is  a data  hold  involved,  even  though  the  data  hold  alters  the 
form  of  the  answer  dramatically. 

Transfer  functions  for  representative  data  holds  are  given  in  Table  % 
The  data  holds  of  greatest  interest  are  the  zero-order  hold  and  the  slewer 


data  hold.  Their  characteristics  will  be  discussed  at  appropriate  points 
in  the  text. 


TAME  *.  SOME  REPRESENTATIVE  DATA  HOLDS 


TRANS PER  FUNCTION 


E.  AN  INTRODUCTION  TO  MULTI-RATE  SAMPLING 


Sufficient  aspects  of  single-rate  sampled  systems  have  been  reviewed 
to  this  point  that  wo  may  now  consider  analysis  of  a restricted  class  of 
multi-rate  sampled  systems.  A basic  multi-rate  sampled  system  element  of 
interest  is  shown  in  Fig.  A. 


Figure  3.  Slow  Input /Fast  Output  Sampling 


It  is  required,  in  Fig.  3,  that  N be  an  integer  so  that  the  output 
sample  rate  is  N times  the  input  sample  rate.  Write  the  output  equation 
for  tiie  continuous  variable  C as 


C GRT 

so  that 

cT'»  roRTiT/I) 

To  further  reduce  Eq.  ?4,  we  merely  notice  that  Fig.  4 is  equivalent 
in  an  input-output  sense  to  Fig.  '*>.  The  T input  sampling  operation  in 
Fig.  ’ is  replaced  by  two  sampling  operations  which  are  equivalent  to  the 


(73) 

(74) 


Figure  4.  Equivalent  Model  for  Fig.  3 

T /N  ft 

original  one.  The  time  functions  [rm  and  RT  are  obviously  equivalent, 
since  the  intermediate  samples  of  the  new,  but  fictitious,  input  [rT  ] 1 N are 
zero.  This  "trick"  does  show  that  we  may  write 


which  is  the  correct  result  (see  Appendix  B).  The  rule  is,  for  sampling 
operations  such  as  those  shown  in  Fig.  7’,  that: 

"The  outer  sampling  operator  [♦]T  N 'operates  through'  the 
inner  ones  if  the  ratios  of  the  inner  sampling  periods  to 
the  outer  s;intpling  period,  T/N,  are  integers." 

Consider  an  example  application.  Lot 


r(  t ) u(  t. ) 

Therefore 


Referring  to  Table  2,  write 


sT  'N 


( 1 


(77) 


(78) 


Note  the  important  points.  Since  z has  been  defined  as  o * esT  N,  the 
pulse  transfer  function  corresponding  to  G(s)  1 /( s + a)  must  represent 
the  newly  defined  sample  period,  T 'N.  Moreover,  the  s -transform  for 
rT  = u^(t)  must  be  adjusted  to  reflect  the  resampling  with  period  T;N. 

This  adjustment  is  also  required  by  the  fact  that  z is  defined  as  esT  N 
as  a result  of  the  resampling."  The  z-transform  for  [ rT J r n r"1  is 
obtained  in  two  steps: 

• Obtain  the  z-transform  for  RT  in  the  usual  manner 
with  z = es^. 

• Replace  z ^ es^  everywhere  it  occurs  in  K^(z)  bv 
sN  - (esT  N)N> 


’In  this  regard,  some  readers  might  find  it  helpful  to  read  the  review- 
material  of  Appendices  A and  B. 


Suppose  the  problem  is  to  find  a continuous  time  response  which  when 
sampled  results  in  (A  A1  for  the  previous  examples.  The  use  of  partial 
fractions  now  requires  an  expansion  of  the  right-hand  side  of  Eq.  78  into 
N + 1 terms,  since  the  denominator  now  has  N+  1 roots.  This  is  not  a 
pleasant  prospect,  especially  if  N is  large,  say  on  the  order  of  TOO  or 
more.  When  N is  large,  partial  fraction  expansion  is  not  the  best  way  to 
proceed.  A recursion  equation  approach  is  preferable.  Rewrite  Eq.  78  in 
terms  of  the  input  RT: 


CT/N 


z = esT/N 


(79) 


Divide  Eq.  79  by  z and  write  the  corresponding  recursion  equation. 


where 


= e 


iT/N 


"n-1 


+ R. 


n,N 


rt 

Kn  ,N 


0 if  Frac  (n/N)  / 0 

R(nT/N)  if  Frac  (n/N)  = 0 


(80) 


(81) 


The  notation  R‘,  N is  interpreted  as  follows.  Form  the  ratio  of  the 
index  of  the  recursion  equation,  n,  to  Nj  and  take  the  fractional  part  of 
the  number  that  results.  If  this  is  zero,  then  nT/N  is  a point  in  time 
for  which  an  input  sample  is  taken.  However,  if  the  fractional  part  of 
this  ratio  is  non-zero,  we  are  at  an  inter-sample  time  point  for  the  input. 
Consequently,  no  input  satrple  is  taken,  giving  a zero  input  value  for  this 
sample  point  in  time. 

— aT 

To  illustrate,  let  T = 1,  e = 0.5  (a  = 0.693 2),  and  compute  the 
transient  response  to  a unit  step  input  using 


/ 0 Frac  (n/N)  / 0 

Rn,N  ~ j (83 

' 1 Frac  (n/N)  = 0 

The  transient  response  is  shown  in  Fig.  5.  Notice  that  the  absence 
of  a data-hold  circuit  forces  the  continuous  output  to  be  comprised  of  the 
summation  of  a train  of  weighted  impulse  responses  of  the  first-order  sys- 
tem l/(s  4 a). 

The  following  is  accomplished  by  implementation  in  terms  of  this  recur, 
sion  equation: 

• The  inter-sample  response  can  be  computed  with 
any  desired  degree  of  time  resolution  desired  by 
increasing  N. 

• Computer  storage  requirements  are  unchanged  as  N 
is  increased.  The  order  of  the  recursion  equation 
(i.e.,  the  number  of  states)  is  not  affected  by  N. 


P.  SLOW  INPOT /PAST  OUTPUT  SAMPLING 
WITH  DATA  HOLD 

The  inclusion  of  n data  hold  in  the  slow  input/fast  output  sampling 
problem  introduces  an  additional  consideration.  Adding  a data  hold  to  the 
block  diagram  of  Fig.  3 gives  Fig.  6.  The  output  equation  is 

CT/N  = (GM)T/NRT  (84) 

where  M represents  the  transfer  function  of  the  data  hold.  It  is  important 
to  notice  that  the  data  hold  is  configured  for  a T second  sampling  period, 
whereas  the  sampling  operation  on  GM  is  for  a T/N  second  sampling  period. 

An  example  will  demonstrate  the  nature  of  the  problem. 


Figure  6.  Slow  Input /Fast  Output  Sampling 
with  Data  Hold 


Let 


— sT 


M 


, G 


sT/N 


s + a 


(85) 


Therefore , 


(GM) 


T/N 


- e-T)a 


;( s + a V 


T/N 


(86) 


F 


s 

t 


Since  z = esT/N,  it  follows  that 

e-sT  = e-(sT/N)N  = z-N 


and 


-sT  zN  - 1 


1 — e 

Moreover, 

a 

s(s  + a) 

so  that  the  transform  table  gives 


1 1 


s s + a 


n _ 1 

\ S s + a/ 


2(1  - e-^T/"' 


Therefore, 


CT/N  = 


Z“1  z-e^T/N  = 7z'-,)(z-e-aT/N 


zN  - 1 zh-e^T/N 


z(l  - ) 

zN  (z  _ l)(z  - e"^T/N) 


or 


.t/n 


. - e-01/"  t 

, , -1  A — (N-1 )\ 

8 - e^T/»  ' 

1 + z + . . . + z J 

(87) 


(88) 


(89) 


(90) 


(91) 


Comparing  Eq.  91  with  Eq.  79,  it  is  seen  that  the  inclusion  of  the 
zero-order  hold  has  introduced  an  additional  multiplicative  function: 

(,  - .-/»)  . I’-f17”)  [,  + Z-1 )j 

z”(  Z — 1 ) z 


i i 


1 


Methods  for  treating  this  additional  complication  when  N is  a given  large 
number  will  be  discussed  in  Section  IV. 

0.  FAST  INPUT /SLOW  OUTPUT  SAMPLING 

The  next  fundamental  multi-rate  configuration  of  interest  is  the  fast 
lnput/slow  output  sampling  configuration  of  Fig.  7.  This  presents  a more 
involved  situation  from  a computational  viewpoint  than  does  the  slow  input/ 
fast  output  case  studied  previously. 


Figure  7.  Fust  Input/Slow  Output  Sampling 


As  always. 


so  that 


C 


gmrt,/n 

(90) 

[gmrt/n  ]t 

(9?) 

In  En.  9?,  the  T operator  can  no  longer  "operate  through"  since  the 
ratio  of  the  inner  sampling  period  to  the  outer  is  not  an  integer.  However, 
it  is  legitimate  to  replace  Fig.  7 with  the  equivalent  block  diagram  shown 
in  Fig.  8. 


Figure  8.  Equivalent  System 


i 


It  is  readily  apparent  that  the  output  is  not  affected  by  intro- 
duction of  the  "phantom"  T/N  sampler,  since  the  T output  sampler  rejects 
all  the  unwanted  samples  of  C^/^.  One  can  write  the  following  equation 
directly,  from  Fig.  8, 


CT  = [GMRt/n]T  = [(GM)t/nRt/n]T 


(9M 


On  the  surface,  this  exercise  does  not  appear  to  have  been  too  fruitful, 
since  the  T operator  still  cannot  operate  through.  Nevertheless,  Eq.  9t 
offers  a significant  computational  simplification.  This  is  so  because  the 
following  routine  procedure  can  now  be  followed: 

• Compute  (GM)T/N  and  rT/N  z- transforms  in  the  normal 
manner,  using  z ^ es^/N  . 

• Expand  the  product  [(GM)T/N rT/N ] in  partial  fractions 
and  use  the  transform  table  to  find  the  continuous 
time  function  which,  when  sampled  with  period  T/N, 
results  in  [(GM)^/N RvN  | # 

• Find  the  z-transform  of  this  sampled  continuous 
generating  function,  for  a sampling  period,  T,  and 
z * esT  to  obtain: 

[(GM)T/NrT/N]T  = [gmrt/njT 

This  process  is  in  the  spirit  of  the  convolution  approach  described  in 
Appendix  B.  The  main  difference  is  that  we  avoid  the  use  of  a complex 
inversion  integral  by  substituting  the  procedure  of  going  from  the  z-domain 
for  a T/N  sampling  period  (i.e.,  with  z * esT/N)  to  the  time  domain  and  then 
back  to  the  z-domain  for  a T sampling  period  (i.e.,  with  z * esT). 

Notice  that  in  the  previous  section  for  the  slow  input /fast  output 
problem,  it  was  possible  to  separate  out  the  transformed  input  from  the 
pulse  transfer  function  for  the  system.  Having  done  this,  it  is  then  pos- 
sible to  write  a recursion  equation  in  terms  of  an  arbitrary  input.  Here, 
a recursion  equation  for  evaluation  of  Eq.  9^  cannot  be  written  until  a 
specific  R is  given.  That  is,  the  recursion  equation  for  must  be  found 
on  a case-by-case  basic.  Nevertheless,  the  basic  method  described  in  this 


i 

) 1 

■ 
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subsection  will  be  valuable  later  for  analysis  of  multi-rate  closed- loo£ 
systems . 

An  example  illustrates  the  computational  method.  Let 


G 


a 

s + a 


M 


, -sT/N 
1 - e ' 

s 


R = sin  bt  (91') 


First,  compute  (GM)T/N  z-transform.  This  is  straightforward,  since  the 
sampling  operation  on  GM  and  the  data  hold  are  for  the  same  sampling  period, 
T/N.  Then  obtain  RT/N(z)  from  the  z- transform  table,  and  form  the  product 

[(GM)T/NrT/N]T/N. 


(gm)t/n 


(1 


s ( s + a , 


T/N 


— e 


-aT/N 


z - e 


■aT/N 


z f esT/N 


(96) 


The  re  fore , 


(gm)T/nrt/n 


I _ e“flT'/N  z sin  bT/N 

z _ e“^  ^ ’ (z  - cos  bT/N)?  + (sin  bT/N):' 


(97 ) 


Expand  in  partial  fractions: 

ft/n  {gm)t/nrt/n 

Az Bz  sin  bT/N  + Cz(z  - cos  bT/N) 

z - (z  - cos  bT/N)1'  + (sin  bT/N)-' 

Compute  f(t)  such  that  f1  /N  *=>  KT/N 

+ 13  sin  bt  + C cos  bt 


sT  N 


(96) 


(oo) 


f(t)  = Ac 


Compute  f'T 


fT  = f(GM)T/NRT/N)T 


Notice  the  clear  resemblance  of  Eq.  98  to  Eq.  100.  The  only  difference 
Is  that  T/N  in  Eq.  98  is  replaced  by  T in  Eq.  100.  This  emphasises  that  the 
intermediate  step  of  computing  the  generating  time  function  is  merely  a con- 
venient conceptual  step  which  may  usually  be  skipped  in  the  analysis  of 
single-rate  systems.  However,  in  a later  section  we  shall  see  that  execu- 
tion of  this  step  can  be  quite  useful  in  sorting  out  delayed  time  functions 
arising  in  the  analysis  of  closed-loop  multi-rate  systems. 

To  complete  the  example,  it  may  he  verified  that: 

A = -€  A0  sin  bT/N 

. /*W 

B = -A0(e  - cos  bT/N) 

where 

, — aT 

a 1 - e 

Ao  — ZTjpw - * ■( 

(e  - cos  bT  /N ) ‘ + (sin  bTH)‘ 

H.  VECTOR  BLOCK  DIAGRAMS 

As  noted  in  the  introduction,  one  objective  is  to  present  all  results 
in  a vector  notation  which  is  compatible  with  either  degree-of- freedom  or 
state  variable  problem  formulations.  For  example,  consider  the  block 
diagram  of  Fig.  9 wherein  the  dimension  of  the  G and  H matrices  are  com- 
patible with  the  dimensions  of  the  C and  E vectors.  One  can  now  develop 
equations  and  pulse  transfer  functions  for  this  closed-loop  system  In  an 
orderly  manner. 


( 1 0 1 ) 
( lO.'l 

( 10M 


As 


s — e 


-uT 


( Bs  sin  b'f  + Cz(z  — cos  tvf) 
(z  — cos  bT)’’  + (sin  bT)- 


esT 


( 100) 


At 


Figure  9.  Simple  Closed-Loop  Configuration 

E = R - HGET  ( 1 Ok 

so  that 

ET  = rt-(hg)tet  (105 

Thus, 

Et  = [I  + (HG)T]~1RT  (106 

Using  Eq.  106,  write  an  equation  for  the  continuous  output  vector  C: 

C = G[l  + (HG)T]“1  RT  (107 

One  can  now  conceptually  add  a "phantom  sampler"  at  the  output  C and  obtain 

CT/N  = gT/N[I  + (HG)T]_1  RT  (108 

by  using  the  results  of  Section  II-E.  Thus,  values  for  the  continuous  out- 
put C can  be  calculated  at  as  many  inter-sample  points  as  desired. 

Consider  another  example.  Suppose  the  open-loop  plant  is  formulated  in 
the  first-order  state  variable  form 


x 


Fx  + Gu 


(109 


where  x is  the  state  vector  and  y is  the  output  vector.  F,  G,  and  H are 
constant  matrices  of  appropriate  dimension.  Laplace  transforming,  we  obtain 

• X(s ) = (Is  - F)  ^Gl^s)  + x(0) ] = A(s)U(s)  + B(s)x(0)  (ill) 


Y(s)  = HX( s ) 


(11.) 


In  Eq.  Ill,  x(0)  represents  the  initial  condition  vector. 

The  open-loop  plant,  given  by  Eqs.  Ill  and  112,  might  then  be  embedded 
in  the  closed-loop  configuration  of  Fig.  10  wherein  and  represent 
compensation  matrices  which  may  be  selected  to  achieve  closed-loop  design 
objectives. 

Suppose  the  objective  is  to  find  an  expression  for  the  ( continuous  1 
state  vector  x.  First,  write 


Therefore, 


= W>T  - WV[HBX(0)  + HAMJ^ 


rprn  TT  rP  TT  rn  rn 

WjR  - W(WC(HB)  x(0)  - Wy(HA.M)  U, 


= [I  + (Wy)(HAM)Tr'  [W>T  - wy(HB)Tx(0)] 


(m) 


(11U) 


(II1') 


u,  / u; 


Figure  10.  A Representative  Closed-Loop  System 


T 

Thus,  if  is  a function  of  the  input  vector  R and  the  initial  condition 

vector  x(0), 

t*  -? 

To  find  the  continuous  state  and  output  variables,  x and  y, 

(Ho) 

(117) 

The  inter-sample  response  can  be  determined  to  any  desired  resolution: 


X Bx(0)  + AMli^ 

Y HRx(O)  + HAMlf[  - Hx 


XT/N 
XT  N 


dt  Nx(o)  -*  (am)t/nu* 

9 > 

PT  Nx(0)  ♦ (am)t  nJ(i  -*  (w^  w^)(  ham)t  f”'  [wV  - W^WT-(HB)Tx(0)]| 

(IIP) 


In  general , throughout  the  text  care  has  been  taken  to  develop  all 
equations  in  vector  notavion,  even  though  most  of  the  illustrative  examples 
are  scalar.  Note  that  this  transform  domain  notation  is  applied  just  as 
easily  to  state  vector  problem  formulations  as  it  is  to  degree-of-freedom 
formulations . 

I.  SWITCH  DECOMPOSITION  CONCEPT 

Tiie  pliant om  sampler  approach  was  introduced  in  Section  II-G  in  connec- 
tion with  the  fast  input /slow  output  sampling  format . Now  that,  vector 
notation  has  been  discussed  (Section  II-H),  we  are  able  to  present  an 
alternative  description  that  includes  the  phantom  sampler  as  well  as  the 
"T  N"  approach  treated  earlier  in  Section  II-K.  This  alternative  approach 
is  called  "switch  decomposition"  (Ref.  o). 

In  essence,  switch  decomposition  is  a procedure  wherein  systems  having 
multiple  sampling  operations  (occurring  at  fixed  but  unequal  sampling  inter- 
vals, but  with  a sampling  pattern  which  is  repeated  over  a fixed,  finite  time 
interval)  are  converted  into  an  equivalent  single  sample  rate  form.  As  ori- 
ginally introduced  by  Kruno . the  method  used  a summing  point  methodology 


wh:  -li  prove,!  to  l>e  extremely  cumbersome  when  the  rut i os  of  the  sumpling 
perl. vis  become  high.  For  this  reason  and  also  because  evolving  state  tran- 
sit ior  metluKls  wort’  tending  to  push  transform  methods  into  the  background, 

••  > method  fell  into  disuse.  However,  there  is  much  to  recommend  the  switch 
i • imposition  concept  for  use  in  both  time  domain  and  transform  domain  analy- 
In  th ■ subsections  that  follow  we  will  review  the  basic  concept  and 
r-mov  xne  • irlicr  restrictions  by  recasting  it  in  vector  form.  Tiie  vector 
''onr.  ••  implif'es  matrix  block  diagram  manipulation  for  multiloop,  multi-rate 
systems;  but,  as  we  shall  see,  the  dimensionality  of  problems  for- 
mulate! using  this  approach  can  present  a practical  limitation. 

J.  VICTOR  SWITCH  DECOMPOSITION 

Consider  the  continuous  signal  shown  in  Fig.  11a  to  be  sampled  at 
'T  samples  'second.  This  results  in  the  sample  sequence  shown  in  Fig.  lib. 
The  sampled  values  have  been  numbered  for  easy  reference.  Suppose  now  we 
sample  the  continuous  signal  with  a sampling  period,  T.  This  results  in 
the  sample  sequence  consisting  of  1 , h , 7,  10,  13,  ...  shown  in  Fig.  lie. 
Define  this  sample  sequence  to  be  R^. 

Next,  advance  the  continuous  signal  R by  T/3.  Then  sample  the  advanced 
signal  with  a sampling  period,  T.  This  results  in  a sample  sequence  con- 
sisting of  7,  8,  11,  lit,  ...  shown  in  Fig.  lid.  Define  this  sample 

sequence  to  be  (es^'R)^.  Finally,  advance  the  continuous  signal  R by 
T/%  and  sample  it  with  a sampling  period,  T.  This  results  in  the  sequence 
consisting  of  3,  6,  9»  1?,  15,  •••  shown  in  Fig.  lie.  Define  this  signal 
sequence  to  be  (e  'R)^  . 

The  significance  of  the  switch  decomposition  concept  resides  in  its 
ability  to  provide  an  alternative  expression  for  the  original  sequence, 
g/5.  This  alternative  expression  for  R^-/^1  consists  of  the  sum  of  R^  , 

(e''^  *'R)'5-  , and  (e-  ST//"'  r)T  each  delayed  by  a time  interval  corresponding 
to  the  advance. 


,T/3  dT  l , sT/3  ,T  -st/3  ^ , ,?sT/3  ,T  -FsT/3 
S'  = R + (e  • R)  e ' + (c  ' R)  e 


(119) 


E luation  119  has  a simple  factored  equivalent  which  is  the  product  of  two 
vectors  and  the  scalar,  R, 


Figure  11.  Decomposition  of  a Sample  Sequence 


T 


where 


and 


3t/5 


— sT/3  — PsT 

e ' e 


W = 


sT/3) 


2s T '? 


W(W,R)J 


i „-sT/?  — 2sT  /? 

1 e e 


esT/3 

2sT/3 


(120) 


(121) 


(122) 


(12?) 


This  result  can  be  generalised  for  any  major  sampling  period,  T, 
which  is  an  integer  multiple  of  the  minor  sampling  period.  That  is, 
whenever 

T = N(T/N) 


where  N is  a finite  positive  integer  and  (T/N)  is  the  minor  sampling 
period.  The  result  for  arbitrary  positive  N is 


,T  'N 


R = W(W  R ) 


J R ^ 


where 


and 


W W(s) 


1,  e 


-sT/N  -2sT/N 


W,  W,(s)  * W ' (— s ) 

In  Eq.  126,  the  prime  denotes  a transpose. 


ho 


— ( N- 1 )sT/N 


(126) 


The  operator  matrices  W and  W*  can  be  used  to  represent  multi-rate 
sampling  operations  in  terms  of  a single-rate  sampling  operation  in  vector 
block  diagrams.  This  is  illustrated  in  Fig.  12. 

Consider  an  example.  Let  R be  a vector  with  three  components.  Let 
the  first  component  be  sampled  with  period  T/6,  the  second  with  period  T/3, 
and  the  third  with  period  T/2.  That  is, 


The  objective  is  to  compute  W in  order  to  obtain  an  explicit  expression  for 
R*  via  Eqs.  129  and  1 2b  (which  is  equivalent  to  Fig.  12b).  For  this  example, 

P = 3 

T is  the  least  common  sampling  period 
T]  = T/6,  To  = T/3,  T5  = T/2 

Nt  - 6,  N2  = 3,  % = 2 

T/6  is  the  greatest  common  subinterval 
N = 6 

Therefore, 


1,  e-^/S,  e-2sl/6)  e— 3sT/6^  e~ Ust/6^  e-5sT/6 

0 

0 

0 

1,  e-23T/6,  e^31/6 

0 

0 

_ 

0 

1,  e~5sT/6 

(132) 

This  example  gives  us  some  insight  as  to  how  increased  dimensionality  can 
complicate  problems  in  practical  application.  Consider  the  vectors  R,  R* , 
and  (W,R)  . These  vectors  have  3,  3,  and  11  elements,  respectively.  The 
vector  (W#R)^  will  have  pXN^  elements  in  general;  whereas  R and  R*  will 
only  have  p elements  each.  This  is  significant  in  that  analyses  will  tend 
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Multi- Rate  Sampling 
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b)  Equivalent  Single  Sample  Rate 
Vector  Block  Diagram 


Figure  If.  Vector  Block  Diagrams  for  Multi-Rate 
Sampling  Operations 

to  Be  conducted  in  terms  of  vectors  like  (W»R)*  in  distinction  to  vectors 
like  R*.  Consequently,  the  potential  for  expanded  dimensionality  in  con- 
nection with  analyses  of  multi-rate  sampled  problems  is  great.  For  example 
consider  a problem  wherein  there  are  two  minor  sampling  periods,  ->9  ms  and 
•HO  ms.  It  is  easy  to  verify  that  the  dimensionality  expansion  factor,  N, 
is  156O. 

Fortunately,  minor  sampling  periods  which  are  so  little  different  are 
not  usually  of  practical  interest.  Furthermore,  it  is  important  to  stress 
the  fact  that  the  increase  in  dimensionality  affects  input  and  output  vec- 
tors, but  it  does  not  affect  the  dimension  of  the  system  state  vector. 

On  the  more  positive  side,  matrix  operations  are  routine.  Consider 
the  system  shown  in  Fig.  t '*> . Once  the  vector  multi-rate  sampling  operation 
in  Fig.  1 An  have  been  replaced  by  the  switch  decomposition  equivalent 
(Fig.  1?b),  analytical  manipulations  are  routine: 

E'J  (W,,R)T  - (W),1IW.,)T(W >GW,)TEf  (1^ 

H7; 


Therefore 


Figure  1%  A Multi-Rate  Closed-Loop  System 


I + (W  j #HW.  ,)T(Wp#GW , )T  1 (W^R)1 


X GW,E, 


Notice  that  the  dimension  of  the  inverse  in  Eq.  131*  is  determined  by  the 
column  dimension  of  Wj.  If  this  dimensionality  is  high,  and  if  the  column 
dimension  of  W ■>  is  lower,  then  we  can  develop  an  alternate  equation  having 
lower  dimension.  The  alternative  equation  is  in  terms  of  the  xj'  vector. 

The  dimension  in  this  case  is  determined  by  the  column  dimension  of  W->. 

X{  I + (Wp.GW,  )T(W|J[W?)T  (W2*GW1)T(W1<(R)T  (136] 


xT  N (gw,)t/n  (wur  )t  - (wunw.OTxf 


A more  concrete  scalar  example  is  shown  in  Fig.  I1*.  The  objective  is  to 
develop  the  pulse  trunsfer  function  relating  and  R^. 
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UawJ 


W = 


T/2  Sampler 


Figure  I1!.  A Scalar  Two-Rate  Example 


To  develop  this  pulse  transfer  function, 


for  E1  in  terms 


of  RT  and  ET. 


start  by  writing  an  equation 


I 


(1--9) 


Coefficients  of  like  terms  in  Eq.  1 *>9  are  collected.  Next,  the  various 
factors  are  expressed  in  terms  of  the  .-.-transforms  for  sampled  time  func- 
tions which  correspond  to  the  Laplace  transforms  inside  each  of  the  sampling 
operators,  (O1.  The  result  is  then  simplified  to  obtain 


1 


(t7’^ 


)T 


7 - O’ 


T 

TP  -t 


nT 


« sT 

' 0 


( 1 ’.O') 


Tiie  re  fore , 


( 1 !<  1 ) 


I 





Once  t:'1'  is  known  as  a 1X111011011  of  RT,  it  is  a simple  task  to  compute  the 
equations  for  the  two  continuous  states: 


and 


However,  note  that  it  was  necessary  to  consider  terms  such  as 


Ft(s) 


T 


0*2) 


(1^3) 


(1U4) 


in  order  to  obtain  the  2- transformed  functions  required  for  Eq.  1 -<0.  Thus, 
a method  for  obtaining  z- transforms  for  functions  of  this  form  is  needed  at 
this  point.  An  approach  using  the  so-called  advanced  ^-transform  is  pre- 
sented in  the  next  subsection. 

K.  ADVANCED  t -TRANSFORMS 


Matrix  switch  decomposition  requires  the  ".-transform  of  functions 
advanced  in  time  by  some  fraction.  A,  of  the  least  common  sampling  period, 

T.  For  example, 

f(t)  >AT'i  . Qi  (t+AT)  < (t+T)  (l^1) 


By  definition. 


y-»  — a(n+A)T  n 

n 0 


so  we  may  write 


t 1 !>o) 


F(^ 


y %-unT  n 
n O 


— fu\I' 


( i !*7 ) 


The  infinite  summation  is  recognized  as  a geometric  progress 
can  be  expressed  in  closed  form.  Therefore, 


ssion  which 


<2 


4=* 


-uAT 

ze 

^Tt 


for  0 < A < 1 


(148) 


Equation  148  can  be  used  t.o  generate  a 

, „ — i a 

same  manner  as  the  z-transtorms  ol  e 
ordinary  z-transforms  in  Section  II-C. 
Eq.  148  take  on  new  values  t.o  generate 
a -•»-  a + ,ib: 


table  of  advanced  z-transforms  in  the 
were  used  to  generate  a table  of 
Thus,  we  can  let  the  parameter  a in 
new  entries  for  the  table.  Let 


— (a+jb)AT 

ze 

— (a+jb  )T 

7 — o N 


f(,ns  bxr  _ -|  ajn  bAT) 


z — e 


<iT 


(1!»9) 


(cos  bT  - j sin  bT) 


Rationalize  Eq.  1 49  by  multiplying  numerator  and  denominator  by  the  complex 
conjugate  of  the  denominator  to  obtain  the  intermediate  result.  F,(z): 


F,(z) 


ze~!^Tfcos  bAT  - ,')  sin  bAT  1 f ( z - e^lT  cos  bT)  - je'"*1'1'  sin  bT 
rz_e~HT  cos  bT)  4 je^T  sin  bT][(z-e'*T  cos  bT)  - je^1  — ^ 


Multiplying  out  the  factor  and  equating  real  and  imaginary  parts  gives 
after  a small  amount  of  trigonometric  manipulation: 


/.  ze_°‘^'rcos  hATz  — e""*41,  cos  b(1-A)T] 

e-a(t4AT)  jcoa  b(t  + AT)  - J sin  b(t  + AT)  <=>  4 ^-a'f  'J„.  . _-AT  ^ 


_ 2c  cos  bTz  4 e 


3 bATz  F sin  b(  1 -A)T] 

— * 5* _.S»T 

- 2e  co s bTz  + e 


This  is  the  advanced  z-transform  for  the  exponentially  damped  cosine  and 
sine  waves. 
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w 

l / 


Next,  let  a = 0 in  Eq.  148  to  obtain  tbe  advanced  z-transform  for  a 
unit  step  function: 


- - 


(1*30) 


This  states  the  obvious  fact  that  a unit  step,  advanced  by  AT  seconds,  has 
the  same  z-transform  as  the  step  function  itself.  Next,  consider  the 
z-transform  of  the  time  function: 


(t  + ATle 


u(t^)  aAT  te~a^'  + ATe~at 


(150 


Since  we  already  know  the  z-transform  of  each  individual  term,  a straight- 
forward computation  results  in: 


.(2—  e~aT)2 


zAT 

^oT 


ze~QArrgflT  + T(  1 — A)e~*lT  ] 


(z  - e^T)2 


(Vt) 


Similarly, 


(t  + AT)2e^t+^) 


= (t?  + 2ATt  + (AT)2)e_aAr 


(153) 


leads  to  the  transform  of  the  remaining  terms  needed  in  order  to  augment 
Table  2 with  an  additional  column  of  advanced  z-transforms.  Table  1*  is  an 
abbreviated  table  that  is  sufficient  for  the  purposes  of  this  report. 

L.  A COMPARISON  OF  VECTOR  SWITCH  DECOMPOSITION 
AND  THE  PHANTOM  SAMPUER 


We  are  now  in  a position  to  compare  vector  switch  decomposition  and 
the  phantom  sampler.  As  already  noted,  the  case  of  slow-input/fast-output 
sampling  poses  no  conceptual  difficulties  (refer  to  Fig.  15),  since 


CT/N  = [gmrt]t'/n  = [gm]t/n[rt] 


(W 


TRANSFORM  TABLE 


■ 


cos  bT)  + (e  sin  bT)  z — 2e  cos  bTz 


. 

I 

I 


Figure  1%  Slow- Input /Fast -Output  Sampling 


Thus,  the  output  CT  N is  represented,  ns  a product  of  two  factors,  one  of 
which  describes  the  properties  of  the  system  and  the  other  describes  the 
input.  Separation  of  system  and  input  characteristics  is  crucial  as  far 
as  conceptual  closed- loop  manipulations  are  concerned.  For  example,  it  is 
necessary  lor  the  use  of  signal  flow  graphs,  block  diagram  algebra,  etc. 
The  state  of  affairs  is  somewhat  different  for  the  fast-input/ slow-output 
sampling  case  shown  in  Fig.  16.  The  output  in  Fig.  16  is  expressed  as 

CT  [GMRT/N]T  [(GM)T/NrT/N]T  (] 


Figure  lo.  Fast-Input/Slow-Output.  Sampling 


where  the  extreme  right-hand  side  is  obtained  by  using  the  phantom  sampler 
concept  (Eii.  9M.  No  simple  product  factorisation  for  Eq.  I'vs  into  a system 
pulse  transfer  function  and  a transform  of  the  input  signal  is  apparent. 
While  vector  switch  decomposition,  depicted  in  Fig.  17.  provides  the  sought - 
after  product  1 actor i oat ion , this  is  at  the  expense  of  an  increase  in  the 
dimensionality  of  the  problem  since 


Figure  17.  Vector  Switch  Decomposition  Model; 
Fast -Input /Slow-Output  Sampling 


Specifically,  the  row  dimension  of  fW^Rl^  and  the  column  dimension  of 
[GW  ]T  are  N for  the  fast- input /slow -output  sampling  case  in  Fig.  1 6. 

On  the  other  hand,  if  G(s)  and  R(s)  are  known  prescribed  functions, 
the  evaluation  of  Eq.  155  can  be  carried  to  completion  using  either  the 
vector  switch  decomposition  or  the  phantom  sampler  concept. 

Consider  use  of  vector  switch  decomposition  for  analysis  of  the  slow- 
input/fast-output  sampling  case  (Eq.  I5U  and  Fig.  15).  The  equivalent 
block  diagram  is  shown  in  Fig.  18. 

CT'/N  = WfW.GMRT]1  = W[W,GMlTRT  ( 1 57 ) 

Notice  that  system  properties  are  separable  from  the  input  as  is  always  the 
case  for  the  vector  switch  decomposition  formulation.  Also  notice  that 
from  an  input/output  viewpoint  [i.e.,  considering  (W[W#GM]T)]  there  is  no 


Figure  18.  Vector  Switch  Decomposition;  Slow-Input/ 
Fast-Output  Sampling 
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increase  in  the  dimension  of  the  pulse  transfer  function  matrix.  If  the 
C vector  has  m elements  and  the  R vector  has  p elements,  then  (W[W#GM]t) 
is  an  m x p matrix.  However,  one  should  also  notice  in  connection  with 
Eq.  1c'7  that  an  intermediate  vector,  p^,  has  been  defined  (conceptually, 
at  least)  which  has  N elements.  It  is  frequently  the  case  that  the  analyst 
iias  no  interest  in  PT.  When  this  is  so,  application  of  vector  switch  decorm 
position  does  not  result  in  increased  dimensions  for  slow-input  'fast-output 
sampling  elements  like  that  shown  in  Fig.  l3. 


We  close  with  a specific  numerical  example  to  demonstrate  t lie  utility 
of  advanced  o- transforms  in  conjunction  with  switch  decomposition.  Consider 
Fig.  Id,  set  up  in  a phantom  sampler  format. 


By  switch  decomposition 


Figure  19.  Illustrative  Example 

Note  the  manner  in  which  w,  written  in  terms  of  delay  factors  in  Eq.  1 59 
is  converted  to  advanced  factors  in  Eq.  160.  Equation  160  is  in  a form 
suitable  for  the  use  of  the  advanced  z- transform 


Equation  160  becomes 


z e_2bT/?  e~bT'/5 

-bT  * -bT  > -bT 

— e z — e z — e 


Therefore, 


-2aT/3 

ze 


r — PbT -nT/3  -W/?  -2aT/3- 

-,T  z[z  + e ' e 1 + e e 

, — aT  \ , — bT  v 

(z  - e )(s  - e ) 


Next,  repeat  the  problem  using  the  phantom  sampler  (Eq.  9M : 


-bT/?  > 


z --  esT  — 


r 

l 


„T/3  t/3 
= G ''  R 


Bz 

-bT  '*■ 


(Id1*) 


It  is  readily  determined  that 


-bT  <1 


-bT/3  — uT  /T  ' A — bT  T*  -nT  ^ 

e — e e — e 


(165) 


Equation  1 6‘-  implies  a continuous  "generating"  time  function  of  the  form 


C(t)  = Ae-^1  + Be“bt 


so  that  upon  resampling  in  the  T frame  time  we  obtain 


CT  = + - B" 

-nT  — bT  ’ 

2 — e 2 — e 


Substituting  for  A and  B and  then  clearing  gives 


(Ido) 


(167) 


— fbT/3  -0T/3  —bT/3  — 2aT /3  ■ 

? e + e ' e 


, -<iT  u„  -bT  v 
(2  - e )(=  - e ) 


( 1d8) 


Equation  I08  is  in  agreement  with  the  switch  decomposition  result  of 
Eq.  Id.'. 

One  may  conclude  that  the  phantom  sampler  approach  is  more  efficient 
than  the  vector  switch  decomposition  approach  when  both  approaches  apply. 

On  the  other  hand,  vector  switch  decomposition  always  provides  a means  for 
analysing  system  properties  independently  of  system  input  signals  which  are 
unspecified  or  when  various  input  signals  are  of  interest.  The  computa- 
tional approach  one  uses  for  any  given  problem  should  be  selected  on  the 
basis  of  these  considerations. 


We  may  also  note  that  the  vector  switch  decomposition  tool  can  be  used 
to  model  the  computational  throughput  delay  of  the  computer.  This  topic  is 
treated  in  Appendix  F- 

M.  SECTION  SUMMARY 

A table  of  abbreviated  transforms  has  been  developed  using  the  transform 
for  an  exponential  function  of  time  together  with  several  transform  domain 
properties.  This  table,  coupled  with  partial  fraction  expansion,  enables 
us  to  move  back  and  forth  between  the  t-,  s-,  and  z-domains  with  relative 
ease.  Basic  background  material  for  analysis  of  multi-rate  sampled  systems 
was  also  developed.  Included  were  discussions  of  the  "phantom  sampler," 
vector  switch  decomposition,  and  the  "T/N"  transform  approach.  Examples 
have  been  used  to  demonstrate  the  strengths  and  weaknesses  of  each  tool  and 
to  call  attention  to  points  which  will  be  developed  in  greater  depth  in  the 
text  which  follows. 
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SECTION  III 


ANALYSIS  IN  THE  V -DOMAIN 


A.  INTRODUCTION 

During  the  19r'0's  and  19<jO's  the  need  for  simulating  open-  and  closed- 
loop  aircraft  responses  furnished  an  appreciable  impetus  for  the  develop- 
ment. and  refinement  of  discrete  algorithms  for  the  simulation  of  continuous 
systems  using  digital  computers  (e.g.,  Refs.  1 1 — 1 *" ) . Thus,  the  Tustin 
transform  method  and  other  similar  techniques  for  the  approximate  discre- 
tization of  continuous  systems  made  it  feasible  to  replace  analog  computers 
with  digital  computers  for  the  simulation  phases  of  a design  effort.  Now, 
this  approach  is  also  commonly  applied  to  design  digital  control  laws  for 
fly-by-wire  aircraft.  First  a continuous  control  law  is  synthesized,  and 
then  it  is  adapted  for  digital  implementation  using  one  of  the  approximate 
discretization  methods  originally  developed  for  simulation  (e.g.,  Refs.  1 - 
!o).  This  procedure  is  called  emulation. 

This  application  of  emulation  for  the  design  of  a digital  controller  is 
motivated  both  by  a fundamental  reliance  on  system  design  criteria  developed 
for  analog  systems  and  the  justifiable  desire,  to  preserve  the  large  body  of 
design  experience  built  up  over  the  past  twenty-five  years.  Emulation  is  an 
approximate  procedure  when  the  Tustin  transform  or  other  "direct  substitu- 
tion" methods  are  used.  In  general,  emulation  procedures  fail  to  account 
for  multiplexer /data  bus  effects  and,  more  importantly,  require  use  of  high 
update  rates  (short  computation  frame  times)  in  order  that  live  inherent 
approximations  be  valid.  Thus  there  is  a need  for  a direct  digital  design 
procedure  which  is  exact  (in  distinction  to  approximate),  accounts  for  the 
effects  of  data  holds,  computational  delays,  etc.,  and  yet  preserves  the 
experience  and  physical  insight  developed  over  the  years  using  conventional 
design  proced\ires.  The  main  purpose  of  this  section  is  tu  highlight  the 
properties  of  a design  domain  wherein  these  objectives  arc  realized. 

Specifically,  attention  is  focused  on  the  w ' -domain  — a domain  related 
to  the  well-known  w-domain  (Ref.  8)  by  an  all-important  scale  factor.  It 


Is  a domain  wherein  the  non-minimum  phase  effects  of  the  sampling  and  data 
hold  operations  can  he  directly  accounted  for  using  conventional  frequency 
domain  design  tools  such  as  root  locus  and  Bode  plots.  These  conventional 
frequency  domain  design  tools  can  be  used  to  considerably  greater  advantage 
in  the  w' -domain  than  in  the  w-  or  z-domuins  because  several  more  powerful 
analogies  between  the  s-domain  and  the  w' -domain  exist.  These  analogies 
are,  in  a sense,  the  key  to  exploiting  the  w' -domain  for  design  purpo;  s, 
making  direct  design  in  the  w' -domain  an  attractive  alternative  to  design 
by  emulation. 

In  the  sections  which  follow,  we  first  review  basic  properties  which 
make  the  w '-domain  preferable  to  the  s-domain  or  w-domain.  Following  this, 
illustrative  examples  are  used  to  highlight  the  analogies  between  s and  w ' , 
Next  we  demonstrate  that,  satisfactory  designs  can  he  obtained  using  conven- 
tional design  approaches  even  for  low  data  rates  where  approximate  discre- 
tization techniques  are  so  seriously  in  error  as  to  be  invalid.  Herein  lie 
the  main  contribution:  recognition  that  tile  w'-domain  models  the  sampling 
and  data  hold  operations  exactly,  regardless  of  the  sampling  rate  employed, 
and  that  the  w '-variable  is  analogous  to  the  s-variable  in  the  sense  that 
all  familiar  frequency  domain  design  concepts,  procedures,  and  interpret a- 
' ions  can  be  carried  over  directly. 

B.  RELATIONSHIPS  BETWEEN  z,  w,  AND  V 

In  the  analysis  of  linear  sampled  data  systems,  use  of  z * exp(sT) 
results  in  transfer  functions  which  are  rational  polynomial  functions  of  z. 
This  is  in  distinction  to  the  correspond i tig  transfer  functions  in  s-domain 
terms  which  are  transcendental  functions  of  s.  It  is  also  the  case  that 
till'  left  half  of  the  s-plane  is  mapped  into  the  interior  of  a unit  circle 
in  t lie  .•-plane.  A drawback  of  t he  . -domain  is  that,  convent  ional  design 
criteria,  such  as  root,  locus  and  Bode  plots,  re  more  difficult  to  inter- 
pret. Moreover,  at  high  sampling  rates,  the  z-platie  poles  and  zeros  tend 
to  cluster  on  the  unit,  circle,  creating  numerical  problems  of  a subst.au- 
t ial  magnitude.  The  region  of  the  z-plane  correspond ing  to  stable  system 
behavior  is  interior  to  the  unit  circle  which  prevents  the  direct  applica- 
tion of  Rout.h's  stability  criteria.  Historically  (e.g.,  see  Ref.  8),  this 


fact  prompted  the  application  of  an  additional  bilinear  transformation 
which  maps  a function  of  z into  a domain  where  the  region  of  stability  is 

once  again  the  left  half  plane.  This  is  the  so-called  w transformation 
defined  by  the  equation: 


4 2 ~ 1 

z + 1 


069) 


One  may  use  root  locus  and  Bode  plot  methods  in  the  w-domain  with  more 
facility  and  insight  than  would  be  possible  in  the  z-domain,  even  though 
each  domain  contains  exactly  the  same  information. 

However,  the  w-domain  still  lacks  other  desirable  properties.  Most 
important  of  these  is  the  property  that  w approach  s as  the  sampling 
interval  approaches  zero.  This  property  is  not  provided  since 

w = i'-  ~ 1 _ eST  - 1 _ ST  + (sT)2/,?.'  + ... 

sT  ~ ■ ' - ( ] 

2 + 1 e +1  2 + sT+  (sT)72I  + ... 

and  in  the  limit,  as  T 0,  w approaches  zero  rather  than  s.  A simple 
scaling  of  the  w-plane  changes  this  situation  dramatically.  Define 


w ' 


T 


z - 1 

z + 1 


is  the  so-called  "w-prime"  transformation  wherein 
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In  the  limit  as  T 


0 in  Eq.  17;',  w'  approaches  s. 


This  property  is  significant  in  that  it 
conceptual  basis  for  defining  a quantity 
which  is  analogous  to  frequency  in  the  s 
more,  the  analogy  becomes  an  identity  in 


establishes  the 
in  the  w' -domain 
-domain.  Pur the r- 
the  limiting 


(171) 


(172) 


case . 


We  are  unable  to  cite  a readily  available  reference  for  the  transformation 
given  in  Eq.  172  even  though  the  relationship  is  well  known  to  many  prac- 
ticing control  engineers. 

One  may,  of  course,  use  root  locus  and  Bode  plot  methods  in  the  w' -domain 
as  well.  However,  the  relationship  between  angular  frequency,  cn,  and  the 
imaginary  part  of  w',  v,  is 

v = -jp-  tan  cjdT/2  (v  = u>  for  |oo|  < jt/2T ) (173) 


The  approximate  relationship  between  v and  is  significant  in  that  the 
designer /analyst  may  regard  v as  angular  frequency  (for  |v|  < 2/T)  during 
qualitative  phases  of  design  development.  Conversion  to  actual  angular 
frequency  units  is  almost  always  unnecessary.  Moreover,  (w')— ^ itself  is 
the  trapezoidal  integration  operator  analogous  to  s ^ for  continuous  systems. 
Finally,  the  unit  delay,  z \ when  expressed  in  the  w’ -domain,  has  break 
parameters  which  are  a function  of  sampling  period,  T,  and  has  the  form  of 
a first-order  Pade  approximation  for  a transport  delay  in  the  s -domain. 
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To  illustrate  another  basic  relationship  which  exists  between  the 
s-domain  and  the  w' -domain,  consider  the  z-transform  for  a continuous  low- 
pass  filter  section 


H(s) 


a 

s + a 


(175) 


obtained  assuming  the  input  signal  to  the  filter  has  been  reconstructed 
using  a zero-order  hold  (ZOH).  The  result  is 
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Applying  Eq.  171  to  Eq.  176  gives 


H(w’)  = 
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The  equal-order  over  equal-order  nature  of  Eq.  177  is,  of  course,  a direct 
consequence  of  the  use  of  Eq.  171.  At  first  one  may  feel  that  the  analogy 
between  the  s-  and  w' -domains  is  weak,  since  a proper  rational  function  of 
s will  always  map  into  a rational  equal-order  over  equal-order  function 
of  w'.  However,  this  is  not  the  case  if  zeros  at  infinity  in  the  s-domain 
are  considered,  for  then  the  s-domain  zeros  at  infinity  correspond  to  the 
"extra"  finite  zeros  in  the  w' -plane.  For  example,  Eq.  175  shows  that  H(s) 
has  a pole  at  s = — a and  a zero  at  infinity.  In  the  w' -plane  the  pole  at 
s = — a is  mapped  into 
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(w'  = -a  for  | a | < it/2T) 


(178) 


while  the  zero  at  s = «>  is  mapped  into  a zero  at  w'  = 2/T.  Obviously,  as 
T 0,  w'  -»■  s,  the  w'  -domain  pole  goes  to  —a  and  the  w'  -domain  zero 
approaches  infinity,  its  proper  s-plane  location.  This  is  a general  result; 
every  pole  and  zero  in  the  w' -plane  has  its  counterpart  in  the  s-plane,  as 
long  as  the  zeros  at  infinity  in  the  s-plane  are  counted. 

Notice  now  the  clear  resemblance  of  Eq.  177  to  its  s-domain  counterpart. 
This  is  in  distinction  to  the  z-domain  counterpart  (Eq.  1 76 ) which  has  a pole 
that  approaches  the  unit  circle  as  T -►  0. 

The  zero  at  w'  = 2/T  is  especially  significant.  This  zero  is  introduced 
by  the  ZOH  used  to  reconstruct  the  input  to  H(s).  It  provides  a nonminimum- 
phase  contribution  which  is  the  effect  of  the  data  hold  sampling  rate  para- 
meter. This  provides  another  major  advantage  of  the  w 'domain  in  comparison 
to  the  z-domain  since  the  effect  of  the  data  hold  and  sample  rate  will  be 
quite  apparent  in  root  locus  or  Bode  plots.  To  emphasize  this  point,  let 
the  input  to  the  continuous  filter  be  reconstructed  with  another  type  of 


data  hold  called  the  "slewer  data  hold."  The  slewer  data  hold  results  in 
constant  rate  output  between  sampling  instants  and  has  no  discontinuity 
at  the  sampling  instant.  Using  the  slewer,  one  computes 
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Or,  using  Eq.  171, 


H(w')  = 
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Observe  that  the  slewer  data  hold  has  both  a pole  and  zero  in  the  w' -plane 
model.  Thus  it  is  seen  that  the  use  of  a zero-order  hold  introduces  a non- 
minimum phase  zero  at  w'  - 2/T  (refer  to  Eq.  177),  whereas  the  slewer  intro- 
duces both  a pole  and  zero.  Finally,  note  that  Eq.  101  reduces  to 


lim 

T — 0 


H(w ' ) 
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+ a 


(10.') 


as  was  the  case  with  the  zero-order  hold. 

Table  u summarizes  the  relationships  between  the  various  planes  for  three 
different  filter  sections  assuming  the  use  of  a zero-order  hold.  One  may 
evaluate  these  transfer  functions,  using  T as  a parameter,  to  obtain  a feel 
for  the  relative  positioning  of  poles  and  zeros  in  the  three  domains. 

Thus  far,  three  important  properties  of  the  w’ -domain  have  been  enu- 
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• The  nonminimum-phase  effects  of  data  holds  become 
clearly  evident . 

• Conventional  methods  for  stability  analyses  continue 
to  be  applicable  since  the  left  half  of  the  w' -plane 
corresponds  to  the  left  half  of  the  s-plane. 

A most  important  property  of  the  w' -domain  which  remains  to  be  shown 

is: 

• Conventional  frequency  domain  synthesis  methods  (both 
scalar  and  vector)  continue  to  be  applicable,  even  when 
the  sampling  rate  is  so  low  as  to  cause  large  differ- 
ences between  s-plane  and  w' -plane  pole  and  zero  loca- 
tions for  a given  plant. 

This  property  is  the  topic  of  the  next  two  subsections. 

C.  ILLUSTRATIVE  DESIGN  PROBLEM  — 

SHORT-PERIOD  AIRCRAFT  MODEL 

Conventional  frequency  domain  design  methods  are  directly  applicable 
in  the  w' -domain.  This  will  be  demonstrated  in  two  stages.  In  this  sub- 
section we  first  consider  the  design  of  a stability  augmentation  system 
using  a representative  model  of  a short-period  aircraft.*  The  use  of  the 
second-order  model  will  allow  us  to  illustrate  clearly  the  manner  in  which 
the  non-minimum  phase  contribution  of  the  data  hold  affects  closed- loop 
system  properties  and  responses.  (Furthermore,  these  effects  are  different 
for  a disturbance  input  than  they  are  for  a command  input.)  The  short- 
period  model  will,  however,  tend  to  have  modal  frequencies  well  below  the 
folding  frequency  even  for  data  rates  as  low  as  ten  samples/second.  Thus, 
the  most  dramatic  effect  that  folding  (or  aliasing)  may  have  upon  the  plant 
w'  transfer  functions  will  not  be  evident.  Therefore,  in  the  next  subsec- 
tion we  consider  a fourth-order  aircraft  model  which  includes  a high  fre- 
quency, lightly  damped  bending  mode  for  which  the  effects  of  folding  are 
dramatic.  This  will  enable  us  to  demonstrate  that  the  w' -domain  continues 


*The  technical  details  of  the  models  used  in  this  subsection  and  the 
next  are  contained  in  Ref.  17. 

i 

A 


to  yield  satisfactory  closed-loop  designs  even  when  the  differences  between 
the  s-domain  and  w '-domain  modal  representations  are  large. 

Consider  first  the  simplified  aircraft  model  of  Eq.  18^: 
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In  Eq.  18%  6e  is  the  control  input,  ag  is  the  gust  disturbance  input,  and 
q (body  axis  pitch  rate)  and  a (inertial  angle-of-attack  component)  are  the 
states  of  the  system.  Suppose  the  continuous  controller  is  designed  using 
t tie  mult  iloop  auaLysis  technique  of  Refs.  l8  and  19.  The  block  diagram  of 
Fig.  CO  defines  the  closed-loop  configuration  with  postulated  compensation 
networks  H],  H. ■,  and  H*,.  In  Fig.  .’0,  Hi,  Hp,  and  are  to  be  determined 
(designed),  R is  the  command  input,  and  x-j  and  Xp  are  used  to  represent  the 
states  q and  a,  respectively.  Application  of  the  multiloop  analysis  method 


*1  (q  or  qg) 


xg(  a or  nz) 


Figure  .'0.  block  Diagram.  Illustrative  Example 
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yields  the  matrix  of  closed-loop  transfer  functions  directly.  The  theo- 
retical details  of  the  multiloop  analysis  method  are  documented  in  Chapter 
of  Ref.  1 8. 


X1 

V*5 

\ 

1 

1 

X.  X.Xo 

N ' + HoH-Ji  V' 

Clp  >-  J ugOc 

R 

Xo 

Nscil-s 
L oe  5 

~\ 

1 

1 

Xp  x^oc, 

N_  •+  H,H,-Nn  c,‘ 
op  i j agoe J 

A i 

(Wh) 


[ < 

' i 


( 

I 


The  various  numerators  of  Eq.  1 Sh  are  found  using  Cramer's  rule  applied  to 
the  Laplace  transform  of  Eq.  1 8^  (see  Eq.  1 8*0 ) . 
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N^;.  for  example,  is  found  by  substituting  the  be  column  on  the  right-hand 
side  of  Eq.  18b  into  the  Sn  column  on  the  left-hand  side  of  Eq.  1 85  and 
evaluating  the  determinant  of  the  array  that  results.  A is  the  character- 
istic open-loop  polynomial  and  is  the  determinant  of  the  left-hand  matrix 
in  Eq.  18%  To  illustrate 
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Figure  J1 . Digital  Control  Block  Diagram 
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by  a sampler.  All  the  remaining  elements  of  Fig.  FI  are  to  be  implemented, 
on  a digital  computer.  For  simplicity  of  presentation,  the  data  hold  is 
assumed  to  be  a zero-order  hold;  however,  the  synthesis  procedure  we  are 
about  to  explore  will  be  directly  applicable  for  other  types  of  data  holds 
as  well. 

Taking  the  z-transform  of  the  first-order  set  of  differential  equations 
given  in  Eq.  187  gives: 
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Equation  19' 
frequency  domain 


may  be  developed  from  Eq.  18^  using  either  time  domain  or 
approaches  (refer  to  Appendix  Fl • We  have  assumed  the  gust 
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input  is  sampled  at  1 /T  samples/second  and  held.  This  is  somewhat  at  vuri 
ance  with  the  physical  reality  of  the  problem  since  the  plant  is  actually 
excited  by  a continuous  disturbance . This  assumption  has  been  made  merely 
because  it  has  a simplifying  effect  upon  this  illustrative  problem.  It  is 
not  essential. 


Because  of  the  sampling  assumption  on  a.g,  one  can  proceed  conceptually 
in  terms  of  a completely  discretized  system;  the  state  vector  has  been 
sampled  and  fed  to  the  digital  computer,  as  has  the  scalar  input  R.  One 
may  think  of  Fig.  00  (and  Eq.  18U)  in  terns  of  the  z-  and  w'-domains  as 
well  as  in  terms  of  the  s-domain.  That  is,  in  Eq.  18U  consider  x^s)  to 

Xo  Xo 

be  replaced  by  x((z),  N=^(s)  by  N^(z),  etc.  Thus,  one  can  proceed  to 
generate  the  z-domain  equivalent  of  Eq.  190,  given  a sampling  rate  (assume 
10  samples/second),  and  upon  making  computations  of  the  characteristic  poly- 
nomial  and  numerators  of  the  first  and  second  kinds. 


-Domain 


The  closed-loop  z-domain  equation  is  not  easily  interpreted  by  methods  use- 
ful for  interpreting  the  closed-loop  s-domain  equations.  To  facilitate 
interpretation'  it  is  our  assertion  that  the  w'-transform  should  be  applied 
to  Eq.  191*,  i.e.,  z = [1  + (T/2)w' ]/[l  - (T/r)w' ].  The  result  is  given  in 
Eq.  19%  The  reader  should  compare  the  numerical  values  for  the  gains  and 
break  frequencies  in  Eq.  195  with  the  corresponding  quantities  in  Eq.  1 9?. 
Close  correspondence  for  many  of  these  numerical  values  should  be  noted. 
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It  is  our  contention  that  the  effects  of  the  sample/hold  operation,  the 
data  rate,  etc.,  are  readily  apparent  in  Eq.  195*  Comparing  Eqs.  195  and 
192  it  is  seen  that: 

• Numerators  of  the  first  kind  pick  up  additional  zeros  at 
2/T  in  such  a way  that  they  have  the  same  degree  as  the 
characteristic  polynomial. 

• Numerators  of  the  second  kind  pick  up  additional  zeros 

at  2/T  in  such  a way  that  they  are  equal  in  degree  to  the 
characteristic  polynomial. 

• Every  entry  in  the  s-domain  equation  can  be  considered  to 
be  equal-order  over  equal-order  if  the  zeros  at  infinity 
are  included. 

• Every  entry  in  the  w' -plane  equation  is  equal-order  over 
equal-order.  Zeros  at  infinity  in  the  s-plane  move  to 
either  2/T  or  some  other  location,  for  example: 

-1.25H3(s)  + l)(l  - ~)h5(w<)  (196) 

• The  numerical  values  of  gains  and  time  constants  in  the 
w' -plane  are  very  similar  to  their  s-plane  counterparts. 
(Although  not  demonstrated  in  this  example,  this  obser- 
vation holds  only  for  modes  having  an  s-plane  modal  fre- 
quency which  is  well  below  the  folding  frequency  (2Ws 
equals  Un/T'i ■ The  fact  that  this  observation  holds  true 
only  under  the  stated  conditions  does  not  limit  validity 
of  the  w' -plane  analysis  techniques  in  any  way  when  the 
stated  conditions  are  not  satisfied.  That  this  is  so  is 
demonstrated  in  the  next  subsection. ) 

Thus,  one  may  proceed  in  the  w' -domain  using  all  the  familiar  synthesis 
tools  of  the  s-domain.  The  direct  digital  design  in  the  w'-domain,  however, 
proceeds  with  additional  explicit  knowledge  of  the  nonminimum-phase  effects 
introduced  by  the  A/D  and  D/A  conversion  through  the  zeros  introduced  at  2/T 
(equal  to  20  rad/sec  for  this  example).  If  these  nonminimum -phase  effects 
become  significant  (the  zeros  move  closer  to  the  origin  as  the  sampling  rate 
is  decreased)  then  they  can  be  regarded  in  exactly  the  same  manner  and 
treated  using  the  same  techniques  as  are  used  in  the  s-domain. 

The  actual  synthesis  will  not  be  carried  out  since  the  prime  objective 
was  to  highlight  the  analogy  between  the  closed-loop  transfer  functions  in 
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the  s-,  z- , and  w'-domains.  A synthesis  will  be  carried  to  completion 
for  the  more  complex  aircraft  model  used  in  the  next  subsection. 

D.  SHORT-PERIOD  AIRCRAFT  WITH  BENDING  MODE 

It  remains  to  demonstrate  the  utility  of  the  w' -domain  in  the  face  of 
substantial  folding  effects.  This  can  be  done  by  modifying  the  example  of 
the  previous  section  to  present  a more  realistic  design  situation.  Speci- 
fically, the  aircraft  model  is  modified  to  include  a lightly  damped  bending 
mode  at  25  rad/sec  (this  will  be  close  to  the  10  samples /second  aliasing 
frequency  of  31.4  rad/sec)  and  to  include  accelerometer  and  rate  gyro  out- 
put equations.  Thus,  we  will  have  one  control  input,  one  disturbance  input, 
four  components  in  the  state  vector,  and  two  components  in  the  output  vector. 
The  details  of  the  model  are  given  in  Ref.  17-  For  our  present  purposes,  it 
is  sufficient  to  modify  our  interpretation  of  the  block  diagram  of  Fig.  20. 
Let  the  bending  mode  state  be  called  x^,  the  output  of  the  accelerometer  (n2) 
be  x-| , and  the  output  of  the  rate  gyro  (qg)  be  x0.  The  open-loop  transfer 
functions  are  given  in  Eqs.  197  and  198: 


-Domain 

-1 17.72[.  l4l , 35.2] [— .031'+,  26.357] 

x^s) 

nz 

^-3. 7975 [3"—  + 1]  [.0152,  24.346] 

x2(s) 

qg 

[-355,  6.675] [.01 • 25] 

w 1 -plane  (T  - 0.1) 


( w'  + 1 

W w’  + 

[ 27.53 

H 13020 

-3.7974  ( 3^22  + + l)  [.069,  107.22]  J 

[.3789,  6.864][.o4l7,  60.1371 

(198 
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In  these  equations  quadratic  factors  [s:'  + 2^8  + ag}  have  been  shown 
in  the  shorthand  form  [£,  o^].  Comparison  of  the  corresponding  numerator 
and  characteristic  polynomial  roots  in  Eqs.  197  and  1 98  leads  to  the 
following  observations: 

• The  short-period  quadratic  is  essentially  the  same  in 
either  the  s-  or  w’ -domains. 

• The  bending  mode  at  25  rad/sec  s-domain  has  been  shifted 
upward  to  60  rad/sec  in  the  w' -domain.  (A  still  lower 
sampling  rate  could  shift  it  downward.) 

• The  nz  numerator,  which  is  equal-order  over  equal-order 
in  the  s-domain  does  not  have  the  zero  at  w'  = 2/T. 

Clearly,  we  now  have  a design  problem  where  the  folding  effects  are  signi- 
ficant. However,  the  closed-loop  design  can  still  be  synthesized  using 
conventional  multiloop  frequency  domain  techniques.  The  following  compen- 
sation was  arrived  at  using  mainly  root  locus  techniques. 


HjH^  = 0.006 

(199) 

H0H3  = 0.008*1 

(200) 

This  translates  back  into  the  z-domain  (and  gives  the 
equation  for  the  discrete  control  law)  as; 

required  recursion 

h1h5  - 0.006 

( 801  ) 

^ 2 - .71878 

(80.') 

The  closed-loop  relationships  are  given  by  Eq.  184,  which  is  valid  for 
both  the  s-  or  w' -domains.  Rather  than  give  numerical  comparison  between 
s and  w'  using  Eq.  I'/’l  (as  was  done  in  the  previous  section),  it  is  our 
preference  to  show  Bode  plots  for  the  closed-loop  system  in  the  s-  and 
w' -domains.  The  Bode  plots  for  the  transfer  functions  of  interest  are 
plotted  in  Figs.  2%  and  24.  Notice  in  particular  the  effect  of  the 
nonminimum-phase  zero  on  the  magnitude  and  phase  plots. 

71 


At  this  point  we  come  face  to  face  with  the  central  question,  "Are 
s-domain  design  procedures  and  concepts  effective  in  the  w'-domain?"  We 
proceeded  with  the  design  in  the  w'-domain  on  the  hypothesis  that  these 
design  procedures  and  concepts  are  effective  and  with  the  knowledge  that 
the  w'-domain  transfer  functions  for  the  plant  properly  and  comprehen- 
sively incorporate  all  the  effects  of  the  data  hold  and  sampling  rate  for 
all  modes.  To  demonstrate  that  this  hypothesis  is  correct  we  must  check 
the  time  responses  for  the  discrete  controlled  system  to  see  if  they  are 
acceptable.  This  is  done  for  the  q variable  in  Fig.  25.  Notice  that 
there  is  no  basis  to  conclude  that  the  performance  of  the  digitally  con- 
trolled system  is  inferior  to  that  of  the  continuously  controlled  one  or 
vice  versa.  In  fact,  the  use  of  an  accelerometer  feedback  gain  of  0.006 
has  resulted  in  a lower  q/ctg  amplitude  ratio  at  lower  frequencies  for  the 
digitally  controlled  system. 


a)  tdeg  Step,  R = Sec,  H3  (s)  = t.O 


(deg/sec) 


Figure  25.  q Transient  Responses  of  Continuous  and 
Digitally  Controlled  Systems 
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E.  SECTION  SUM4ABY 


Analogies  between  system  form, nations  in  the  s-  and  domains  have 
been  drawn.  To  arrive  at  the  .'-plane  formulation  one  must  first  discre- 

or;r,P;rem  * °f  8 ''8lid  Umt**  if  effects 

the  d-  '°IT  “r°  t0  be  rePr,iSe"ted  T»ib  beads  to  a statement  of 

the  diseretised  problem  in  the  s-dcmain.  The  s-domain  statement  of  the 

problem  is  then  converted  to  a V. domain  statement  by  means  of  a bilinear 
algebraic  transformation. 

It  has  been  demonstrated  that  direct  digital  control  law  synthesis  in 
he  w -domain  is  a viable  and  practical  alternative  to  design  by  emulation 
C a .«  muons  system.  Key  properties  of  the  w '-domain  have  been  stated 
an  e Visibility”  of  data  hold  and  sampling  rate  non-minimum  phase 
Ciuct..  m the  w '-domain  has  been  demonstrated.  Most  importantly  it  has 
been  Pointed  out  that  conventional  frequency  domain  design  procedures,  TOh 
as  multiloop  analysis.  Bode  plots,  root  locus,  etc.,  are  valid  and  useful 
procedures  in  the  V -domain  even  in  the  presence  of  significant  aliasing. 

.no  ly,  there  is  the  convenience  resulting  from  the  fact  that  the  imaginary 
Part  o.  w , v,  approximates  angular  frequency,  for  |,„|  < ,/;T  or  , J „ 

* °f  tMS  l*  th8t  the  ««W  designer  can  now  synthesise  digital 

controllers  using  considerably  lower  sampling  rates  than  are  required  when 
an  emulation  design  approach  is  used,  furthermore,  the  direct  digital  con- 
trol law  synthesis  approach  presented  here  requires  no  new  analytical  tech- 
niques beyond  those  classical  frequency  domain  procedurea  required  by  the 
emulation  design  approach.  The  analytical  techniques  are  merely  applied 
and  interpreted  in  the  novel  manner  which  we  have  described  in  this  section 
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SECTION  IV 


MULTI-RATE  TRANSFORM  DOMAIN  APPROACH 


A.  INTRODUCTION 


The  existing  state  of  affairs  with  regard  to  the  multi-rate  problem 
(introduced  in  Section  II)  can  be  reviewed  with  the  aid  of  Fig.  26.  It 
may  well  be  that  our  major  interest  in  Fig.  26  is  the  response  of  the  con 
tinuous  output  C in  response  to  the  input  r(t).  This  can,  of  course,  be 
computed,  since 


A Basic  Multi-Rate  Configuration 


For  example,  if  G(s)  = l/( 


then  we  may  write,  setting 


Equation  20b,  while  cori'ect,  is  not  very  illuminating.  We  can  gain  more 
insight  and  increase  computational  efficiency  by  looking  at  the  inter-sample 


response  at  a finite  number  of  points.  This  is  indicated  by  the  use  of  a 
sampler  working  in  a T/N  frame  time.  The  basic  property  of  the  transform 
domain,  discussed  in  Section  II,  can  be  invoked  to  give 


c1*  = [arfrf/"  - gt-V 


(=06) 


Letting  z = esT/N  gives,  for  our  illustrative  example, 


.T/N 


— "T/N 


(f07) 


We  could  pursue  a cpuf'se  of  action  that  would  yield  a "closed  form"  answer 
by  finding  the  N + 1 poles  of  Eq.  :X>7,  expanding  in  partial  fractions  and 
inverting  the  result  back  into  the  time  domain.  However,  this  "increased 
dimension"  approach  is  computationally  burdensome,  especially  if  the  order 
of  N is  high. 

Suppose  we  elect  to  proceed  via  the  recursion  equation  route.  What 
additional  complexities  are  introduced  by  the  multi-rate  nature  of  the 
problem"  Suppose,  in  Eq.  f07,  that  N is  unity  so  that  the  recursion  equa- 
tion becomes 

CnT  = e Tc(n-1)T  + rnT 

In  practice,  one  usually  makes  a mental  note  of  the  frame  time  and  sup- 
presses the  use  of  it  in  the  recursion  equation: 


But,  i”  ^ merely  indicates  that  r(t)  is  being  sampled  every  T seconds  and 
therefore  non-sere  values  result  only  1'or  every  Nth  occurrence.  Therefore 


It  is  seen  that  the  memory  requirements  of  Eq.  .°10  are  basically  unchanged 
from  Eq.  .'09.  There  is,  of  course,  the  calculational  logic  needed  to  imple- 
ment Eq.  ’ll.  For  the  example, 


n/N  * integer 


n/N  integer 


While  this  pleasant  state  of  affairs  remains  essentially  unchanged  when  zero 
order  holds  are  used  as  couplers,  the  use  of  higher-order  data  holds  (such 
as  the  slewer)  does  significantly  increase  the  modeling  complexity.  These 
cases  are  discussed  in  succeeding  sections. 


B.  USE  OF  ZERO-ORDER  HOLDS  AS  COUPLERS 


A'  to  include  a zero-order  hold  (fOH)  gives  Fig.  As 

use  of  the  zero-order  hold  will  cause  no  increase  in  the 


hall 


dimensionality  or  storage  requirement 


igurc 


1 -e’*T 



G(s) 

s 

From  Fig.  27 


r , -sT v 

C = G(s)  ,e l rT 

s 

Evaluate  Eq.  213  at  a sampling  interval  of  T/N: 


(E13) 


CT/N  = | IT  - e~*T)rofs)] 

s 


I T/N 


RT 


(E14) 


Let  z = e^1 /N  so  that  Eq.  21 L becomes 


CT/N 


ZN  _ 
^N~ 


T/N 


or 


CT/N  „ j£^  f(|  - a~-T/NUafs1  ,p/M 

“ (2  -1)1  S | 


(::i5) 


In  Eq.  215  we  have  "inserted" 


1 = 


(2  -TT 


i£-rJ)  _ 


-sT/N 


) 


zr thc  bmcketed  tcm  t,,e  same  r°™  - - — — 

Equation  1r^  simplifies  to: 


CT/N 


fl  + z'1  + ...  z-(N~->) 


sT /N  ,T/N 

•Lr-°  T>g(*)1  rT 


(216) 


Let  the  power  series 


operate  on  the  input  RT: 


80 


I | 


-s-r  yVoC^liV 


« „-i  + • ••  + * 

rt  i + a 


{?M) 


— £>  JL  / »»  |HiV  O / j_l  t 

\ (1^ — u^m  1 „ r(s'),  on 

c'I''S  ' 3 s the  *-«»',st0rB  ling  iBter- 

is  recognise4  - " the  T/« 

. 21  f ^ UoX^-j  *■  rn  -result  ^ 

. + term  in  _ ■»er°'or^er  that  the  T 

We  C tlo„  of  smoothing  W ,„me  »»« ; , a partlcu- 

«*  — ?*?*  u co^4 *' 

val-  '"'“"’next,  observe  that  time  domain: 

cobv»“4-  ,nterpTet»fi»"  « (S'M 

lorly  »»“  ' T t * ...  +f„-!H-') 

-(R-0^  — Rn+Rn-1 

Tr  ^ + • ••  + • ndex  tafc.es  oh 

.oero  »1»  — ***Z  » times  ~<«* 

lc  only  has  “ "0"  the  T v»l“  18  " rore  « can  n°te 

T for  example,  Thus,  tl  Therefore, 

*“l  ! is  n multiple  «f  4-  „ce  of  Ed-  -'8- 

value  »hich  t>  ,f  ^ signif^n 

it  1=  ***•*>  1 , \ 

the  foU^in&-  « , R^-N  I 

z„  Rr/„/TjYLrJ  I 


*h.N 

] ^tCn/^TNT 


^ 1 

-1  + ..•  + z 

*-Vl'  + " 


•w „ , „ « — 

ing  index  diVldedmuiyUplying  tbl®  B8Bie  pr°ce‘ 

inEq* is.  ^ «*»**»' 

phe  value  0 ion  equation  seeKin6  the  1 vioull  eive 

,o  find  the  recurs44"  y ^ „cre  see  ^ looh-up  in  a 

die  as  in  partial  ^ 

0(s >/•  «=“14  teolonents,  -*  *£  “ of  the  *»  ^ “ 
the  ihdlvldua4  e *®  resolt  uould 

rational  format- 
an  mth  order  system)  • 


Rn-m,^ 


*f  • • • -f- 


°T  ■ f^O  -e-sT)jT 


- , „ m-1 

,T  _ aQz  + aj2 


z"  + + ...  ~ R(Z)  (22°) 


Z2 ,Eq;  220  the  reCUrSl°"  eWati°n  1S  "riUe"  "»<**  <**  nT,  the  time 
index,  be  represented  by  n): 


:n  - -blCn_,  - boCn_ 


bmC 


cn-m  + a0Rn  + alRn-1  + •••  + a,nRn.m  (221 ) 


The  recursion  equation  for  eft*  has  identically  the  same  form  as  *.  £2, 

„7  rC°re  th<!  a'S  Md  b'S  <the  -to,  . T/N  sampling' 

nterval.  Now  let  the  index  nT/N  be  n and  write: 


n = -b1Cn-l  ~ b2pn-2  - •••  - bfflCn.m  + 


aORn,N  + a1Rn- 1 ,N  + •’•  + a^_. 


m,N 


(222) 


prefer  toT^l  ^ “*  T/"  ““  iS  «"  — 

We  illustrate  with  several  examples. 

Example  1 


I 


Let 


G(s) 


a 

s~+  a 


(22?) 


The  T recursion  equation  is 


Cn  = e-^Cn.-,  + (1  - e^T)  Rn-l  , 


z = e 


sT 


whereas  the  T/H  recursion  equation 


is 


(22b) 


°»  ’ Ml-.-/*)*,,,  , s = esT/N 


-qT/N-v  „t 


(22^) 


8 


s + a 


so  that 

cT/B  = ? - [1  - BT[,  + ,-1  + ...  + *-<»-!)  (227) 

z — e ' 

The  T recursion  equation  is 

Cn  = + Rn  + - 1 “ ~ e^T)Rn_i  , z = esT  (228) 


whereas  the  T/N  recursion  equation  is 


(229) 

In  the  next  section,  these  basic  ideas  will  be  extended  in  a manner 
that  permits  the  analysis  of  multi -rate  closed-loop  systems  containing 
zero-order  holds.  Before  proceeding,  it  will  be  instructive  to  exercise 
Eq.  225  for  a cosine  wave  input.  This  is  done  in  Fig.  28  for  the  case 
where  T = 1,  thus  causing  a sampling  frequency  of  2n  rad/sec.  The  input 
frequency  in  Fig.  28a  is  jt/2  rad/sec,  which  bears  an  integer  relationship 
to  the  sampling  frequency.  On  the  other  hand,  the  input  frequency  in 
Fig.  28b  is  1-5  rad/sec,  giving  a ratio  of  the  sampling  frequency  to  the 
input  frequency  which  is  an  irrational  number.  The  effect  of  tnis  on  the 
"steady  state"  response  is  clear.  In  the  steady  state,  Fig.  28a  has 
acquired  the  additional  attribute  of  periodicity,  whereas  Fig.  28b  shows 
the  modulation  caused  by  the  non-integer  relationship  between  the  input 
frequency  and  the  sampling  frequency. 

The  concept  of  the  continuous  frequency  response  of  a discretely  excited 
system  will  be  developed  in  Section  V.  This  will  provide  the  tool  for  iden- 
tifying the  spectral  components  of  waveforms  such  as  those  shown  in  Fig.  28. 
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Figure  28.  Response  to  a Cosine  Wave  Input 

C.  MULTI-RATE  FORMULATION 

The  type  of  recursion  equation  developed  in  the  previous  section  can 
be  used  to  advantage  in  the  analysis  of  closed- loop  systems.  To  develop 
the  basis  ideas,  we  note  that  the  transform  approach  states  that  if 


C(s)  = F^sjT^s) 


(230) 


CT/N  = fT/NfT/M 


(231) 


given  N/M  is  an  integer  (see  also  Appendix  B).  The  impact  of  this  can  be 
appreciated  by  considering  Fig.  2$. 

The  output  equation  is,  letting  z = esT/N, 


cT/»  , [,(1  -e-61/*1;  a(s) 

S 


T/N  T/M 


(232) 


f 


Figure  P9.  An  Open-Loop  Multi-Rate  System 


Rearranging  gives 


Again,  let  the  power  series  in  z 


and  then  interpret  the  time  domain  implications 


R„  , for  example,  only  has  a non-zero  value  when  the  n index  takes  on  a 
value  which  is  a multiple  of  T/M.  Thus,  the  T/M  value  is  "held"  N/M  times 
before  it  is  updated.  Therefore,  we  can  write 


RT/Mt , + + . . . , a-fO»/M)- J ] 


r = Rf 

|N_/  n \ _T_  1/  n \ _T_ 

\w\mhm  n \\mu  m 


rt/m, 

Hn,N/M 


(237) 


In  Eq.  237,  T/M  is  the  sampling  interval  of  the  input,  n is  the  running 
index  of  the  recursion  equation,  N/M  (which  is  an  integer)  is  the  number 
of  indexing  periods  for  which  the  input  is  held  constant  before  updating. 

D.  A CLOSED-LOOP  MULTI-RATE  APPLICATION 

The  type  of  recursion  equation  developed  in  the  previous  section  is 
very  useful  in  the  analysis  of  closed- loop  systems.  To  illustrate,  con- 
sider the  problem  given  in  Fig.  30.  From  Fig.  30  we  can  write  the  trans- 
formed equations 


1 - e~ST/ ' „T /? 


1 — t /? 

Y - -sTsVrr x 


(238) 


Applying  the  results  of  the  previous  section: 

xT/»  . et/5  |,  , a->  + ...  t z-t(»/3)-0|  ( „ „ , (2W) 

1®/"  = 1 — 0 J-.  XT/g  |l  * z~'  + ...  + l|  j (21.0) 


Figure  30.  Multi -Rate  Example 
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Any  desired  integer  value  of  N can  be  used  which  is  commensurate  with 
3 and  2 (for  example,  6,  12,  18,  etc.).  The  difference  equations  can 

be  written  directly  from  Fig.  30.  Noting  that  z = es^/^  and  E = R — Y 
gives 


V _ Y -i;/5  , + X rT/3  , 

xn  ' xn-1  TYn-1,N/3  N Rn-1,N/3 


Yn  " e 


N/2 


(241) 


where,  for  example. 


(?h2) 


Note  that  these  two  recursion  equations  completely  define  the  system  no 
matter  what  value  of  N is  selected.  The  response  of  this  closed-loop 
system  to  a sine  wave  input  is  shown  in  Fig.  31.  Since  the  sampling  fre- 
quencies (T  = l)  and  the  input  frequency  of  >t  are  commensurate,  we  see  that 

the  "steady  state"  in  Fig.  31  takes  on  the  additional  property  of  periodi- 

* 

city.  This  would  not  be  the  case  if  the  periods  of  the  input  wave  were  not 
commensurate  with  the  sampling  frequencies. 


Figure  M.  Transient.  Response,  Cosine  Wave  Input. 


As  mentioned  previously,  means  for  identifying  the  spectral  content 
in  response  to  sine  and  cosine  wave  inputs  will  he  developed  in  Section  V, 


It  is  important  to  note  that  the  basic  simplicity  of  the  multi-rate 
results  is  critically  dependent  on  the  use  of  zero-order  holds.  An  indi 
cation  of  the  type  of  complications  which  result  from  the  use  of  more 
sophisticated  couplers  will  be  given  in  the  next  subsection,  where  we 
treat  the  slewer  data  hold. 


I.  THE  SLEWER  DATA  HOLD 


The  impulse  response  of  the  familiar  idealized  triangular  data  hold 
is  shown  in  Fig.  32.  The  transfer  function,  in  the  s-domain,  is 


Figure  .32.  Impulse  Response  of  Triangular  Data  Hold 


If  the  impulse  response  of  Fig.  32  is  delayed  in  time  by  T seconds,  the 
unrealizable  triangular  data  hold  becomes  physically  realizable  and  is 
called  the  slewer  data  hold,  a description  first  coined  by  Goff  (Ref.  10) 
Since  a time  advance  of  T seconds  is  introduced,  the  transfer  function 
becomes 


Thus,  the  transfer  function  of  the  slewer  data  hold  is  simply  the  square 
of  the  transfer  function  of  the  ZOH  divided  by  T. 

In  the  presence  of  no  additional  dynamics,  the  implications  of  Eq.  2^ 
are  clear  — one  simply  samples  the  input  waveform,  advances  this  value 
T seconds,  and  ramps  ("slews")  to  it  (refer  to  Fig.  33).  The  basic  property 
of  the  slewer  is  that  it  provides  ramp-like  waveforms,  with  no  jump  discon- 
tinuities at  the  sampling  instants.  An  inherent  lag  of  one  frametime  is 
introduced.  This  additional  effective  lag  of  one-half  a period,  when  com- 
pared with  the  ZOH,  must  be  ;aken  into  account  in  the  synthesis  of  closed- 
loop  control  systems. 

It  will  be  helpful  to  review  the  development  of  the  ZOH  recursion  equa- 
tions in  order  to  place  the  differences  and  similarities  encountered  in  the 
use  of  the  slewer  into  perspective.  The  basic  open-loop  equation 


= (OMf/V 


(2*5) 


reduces  to 


G R1 


(214-6) 


Figure  33.  Reconstruction  Via  the  Slewer 
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if  M is  permitted  to  be  a ZOH  and  N is  set  equal  to  unity.  More  specifi- 
cally, rewrite  Eq.  : 'U6  as 


C 


T 


) 


I 

I 


RT 


(2**7) 


where  1 — z * is  constrained  to  operate  within  the  "curly"  brackets.  This 
constraint  is  exercised  because  it  is  recognized  that  the  [G(s)/s]^  con- 
tains the  pole  (1  — z”').  Failure  to  recognize  this  by  implementing  the 
recursion  equation  via 


C 


T 


(1*8) 


adds  an  additional  state  which  is  (eventually)  cancelled  out  by  a "zero" 
operating  on  the  input. 

In  tiie  T/N  case  for  the  ZOH,  set  z - es^/"  and  write 


,T  'N 


IT  /N 


0 


- ) ! rT  iczi 


1 — Z 


-N 


(■*9) 


The  1 - z 1 inside  the  curly  bracket  takes  care  of  the  pole  of  [c(s  )/s  ]T/^ 
at  z 1 , and  the  term 


+ z“(N~'  > 


(-50) 


operates  on  tiie  input,  yielding  the  interpretation 


=£>  R 


T 

n.N 


(251 ) 


That  is,  the  zero-order  hold  results  in  the  "T"  sanpled  function  being  held 
at  a constant  value  over  the  entire  "T"  frame  time. 


Now  consider  the  slewer  where  a similar  situation  pertains.  Setting 
M = Ms  and  N = 1 gives 


- to' 

Tsc 


(1  - Z-1  )2  J 


(252) 


In  Eq.  252,  the  term  (1  — z—1 )‘  is  constrained  within  the  curly  brackets 
T 

because  [G(s)/Ts-  ] will  always  produce  a second-order  pole  at  z = 1 , and 
hence  should  be  cancelled  out  at  the  onset.  Letting  N be  non-unity  gives 


CT/N  _ \ 

' I 


£LlL 


(T/N)s^ 


1 n[i  - z-1  r 


z = esT/N 


(255) 


Now  the  term  in  curly  brackets  is  simply  the  "T"  transform  recomputed 
with  a T/N  frame  time.  That  is,  if  we  know  how  to  compute  the  pulse  transfer 
function  using  a slewer  and  a T frame  time,  then  a recomputation  for  a T/N 
frame  time  is  a straightforward  operation.  The  term  operating  on  the  input 
can  be  rewritten  as 


N 


[1 

[1  - z-1  f 


[l  + z ^ + . . 


N 


+ z-(N-i>r 


(25M 


The  first  [l+z  +...  +z operates  on  the  input  creating  the  con- 
stant value  over  the  T frametime.  Now  the  additional  [1  + z—1  + ...  + z (N— . ) ] 
factor  no  longer  operates  on  inter-sample  points  that  are  zero  but  rather 
sums  the  previous  constant  value  with  the  present  constant  value.  The  only 
complication  is  that  not  all  the  previous  N-1  constant  values  are  in  the 
same  T frame;  we  now  span  values  taken  from  two  adjacent  frames.  The  situa- 
tion is  depicted  in  Fig.  for  the  case  where  N = 5.  Setting  z = esT/5, 

T 

we  first  construct  Rn^  . Two  cases  will  demonstrate  the  additional  logic 
required.  In  Case  1 we  add  three  constant  values  from  the  frame  where 
ZT  i t < 3T  and  two  constant  values  from  the  frame  where  T < t i 2T.  In 
Case  2,  we  add  five  values  from  the  frame  4T  i t < 5T.  Thus  it  is  seen  that 
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Figure  5**.  Slower  Sampling  Logic 


the  additional  complication  introduced  is  the  need  for  an  indexing  register 
call  it  k,  which  keeps  track  of  numbers  summed  in  the  present  frame  and 
those  summed  from  the  past  frame.  We  adopt  the  notation 


where 


The  index  k varies  between  ) and  N and  increments  in  unison  with  n,  the 
running  index  of  the  recursion  equation.  It  is  reset  to  unity  whenever 
^ takes  on  a new  value.  That  is, 

k 1 if  (n/N)™™  = Integer 


(257) 


or  equivalently 


1 if  (n/N)™ ap  - 0 


(258) 


T T T 

In  Eq.  2Sb  R is  the  first  "back  value",  while  R - R „ „ represents 

n-N,N  * n,N  n-N,N  * 

the  first  back  difference. 

The  reader  should  bear  in  mind  that  the  pulse  transfer  obtained  using 
the  slower  will  always  contain  a free  z in  the  denominator  which  will  even- 
tually end  up  operating  on  the  input;  hence,  it  is  prudent  to  point  out  that 


-lnTr,  -t  — (N-  ? ) 

z R [ 1 + z + • • • + z ' 


F.  ILLUSTRATIVE  EXAMPLE 


^ Rn-(N+1),N  + N Rn-1,N  “ Rn-(N+1),N 


An  example  will  serve  to  bring  the  details  of  the  discussion  into  focus. 


G<s>  = FTT  * 


so  that 


2T/N  = | 


1-5-  s2(s  + o 


(i  - z-1  r j rt  Il+e  \+,w+,~ 


( 2t>0) 


\ 


Therefore, 


The  recursion  equation  is  written  directly  as 


-n 


a0  C , + a 
2 n-1  0 


C / T m 

+ T (Rn-1,N  ~ Rn- 


(HN),N 


+ a. 


ST  k2  / T T \ 

n-(2+N)  + N \Rn-2,N  ~ V(2+N),N/ 


In  Eq.  26^*,  the  indexes  k1  and  k , vary  between  1 and  N; 


'< 


\ 


k,  — 1 

when 

[(n  - 1 )/N]frac 

ko  1 

when 

[(n  - 2)/N]p^AC 

(265) 


An  actual  computation  of  the  terms  in  Bracket  B is  unnecessary,  since 
Bracket  B is  Bracket  A delayed  by  the  factor  e~ . If  N is  set  to  unity 
(recall  that  uq,  a^,  and  a - are  functions  of  N),  we  obtain  the  "T"  recur- 
sion equation  as  a check: 

Cn  = a2Cn_i  + a0Rn_i  + Rn_2  (266) 

Using  a sine  wave  input,  the  open-loop  characteristic  of  the  slewer  can 
be  contrasted  against  the  ZOH.  This  is  done  in  Fig.  55.  Note  the  added 
attenuation  and  phase  shift  introduced  by  the  slewer.  On  the  other  hand, 
the  output  is  considerably  smoother.  We  remark  at  this  point  that  the 
"strength"  of  the  slewer  lies  with  its  magnitude  plot,  whereas  its  weak- 
ness is  linked  with  the  phase  characteristic.  This  point  will  become 
clearer  in  Section  V. 


Figure  55.  Comparison  Between  Slewer  and  ZOH 
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G.  A MOTIVATING  EXAMPLE 

We  have  now  demonstrated  the  application  of  the  T/N  approach  both  to 
open-  and  closed-loop  systems  and  to  multi-rate  systems.  The  primary 
advantage  of  the  technique  is  that  the  inter-sample  response  of  multi- 
rate systems  is  obtained  without  a significant  increase  in  the  complexity 
of  the  defining  recursion  equations.  The  importance  of  being  able  to 
compute  the  inter-sample  points  will  be  demonstrated  with  the  use  of  a 
"motivating"  example.  We  use  this  particular  adjective  because,  even 
though  a step  response  is  discussed,  our  true  intent  is  to  interest  the 
reader  in  the  idea  of  finding  the  "continuous"  frequency  response  of  a 
discretely  controlled  system. 

Consider  the  block  diagram  of  the  single-rate  closed-loop  system  shown 
in  Fig.  36.  In  Fig.  3 6,  let 


G(s) 


10 

s + 10 


or  G(s) 


so  that  the  open-loop  plant  is  either  strongly  stable  or  unstable.  Let 
M be  either  a slewer  or  a ZOH.  As  one  of  our  design  objectives  we  will 


Figure  3 6.  A Closed-Loop  Single-Rate  System 
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require  a closed- loop  pole  in  the  z-plane  that  has  (s  + 0.5)  as  its  s-plane 
counterpart.  The  other  objective  will  be  to  force  a unity  steady-state 
1 step  response.  In  addition,  a 1 sec  frame  time  will  be  utilized,  forcing 

a sampling  frequency  of  Pit  rad/sec  and  a folding  frequency  of  n rad/sec. 
Hence,  the  folding  frequency  will  always  be  lower  than  the  modulus  of  the 
open-loop  characteristic  equation. 

■j  The  details  of  picking  the  compensation  networks  G^  and  Go  which 

achieve  these  objectives  are  discussed  in  Appendix  H.  The  actual  design 
can  be  easily  carried  through  in,  among  other  domains,  the  w' -plane.  The 
resultant  compensation  is  given  in  Table  6.  (Note  that  the  slewer  design 
was  carried  through  for  only  one  value  of  the  open-loop  root,  since  the 
inter-sample  response  is  quite  smooth  and  therefore  very  similar  for  either 
case . ) 


TABLE  6.  COMPENSATION  NETWORKS 


DATA 

HOLD 

OPEN-LOOP 

POLE 

Gi(z) 

Go(  z ) 

ZOH 

(s  - 4) 

-.00734 

137.2 

ZOH 

(s  + 10) 

•393 

- 1.541 

Slewer 

(s  + 10) 

.i+372z 

z + . 1 1 1 

- 1.541 

The  step  responses  are  shown  in  Fig.  ^7.  Note  that  at  the  sampling 
instants  all  three  closed-loop  designs  give  exactly  the  same  performance. 
However,  the  inter-sample  responses  present  quite  a different  picture. 
There  is  a very  noticeable  roughness  or  "ripple"  in  the  ZOH  designs, 
whereas  the  slewer  response  is  smooth  and  well  behaved.  In  addition, 
notice  that  it  was  possible  to  overcome  the  additional  phase  lag,  intro- 
duced by  the  slewer,  through  the  use  of  a lead  filter  for  Gj(z). 

Suppose  this  system  was  excited  with  a sinusoidal  input  signaL.  In 
the  next  section  we  will  show  that  the  classical  sampled  spectrum  of 
z-trnnsform  theory  would,  like  the  step  response  evaluated  at  the  sampling 


1 
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Stewer 


Figure  Closed-Loop  Step  Responses 


instants,  be  unable  to  differentiate  between  the  three  designs.  This  is 
so  because  the  sampled  spectrum  considers  only  the  single  sinusoid  that 
matches  the  output  at  the  sampling  instants.  Our  objective  in  Section  V 
will  be  to  develop  a tool  that  can  define  the  spectral  content  for  the 
continuous  states  of  the  system.  Such  a tool  would  have  a capability  to 
clearly  differentiate  between  the  three  closed-loop  systems  of  Fig.  ?'7. 


H.  SECTION  SUMMARY 


A classical  property  of  sampled  data  control  theory  has  been  exploited 
in  order  to  develop  a creditable  tool  for  the  analysis  of  the  inter-sample 
responses  of  multi-rate  systems  under  the  influence  of  discrete  controllers 
In  particular,  it  was  shown  that  the  use  of  Oil's  leads  to  recursion  equa- 
tions capable  of  defining  the  inter-sample  response  with  any  desired  degree 
of  fineness  — yet  the  order  of  the  recursion  equation  is  not  increased. 


The  characteristics  of  the  slower  data  hold  were  also  developed.  Here 
a noticeable  increase  in  the  logic  requirements  was  needed,  even  though  the 
order  of  the  recursion  equation  is  not  increased  (when  compared  io  the 
single-rate  recursion  equation). 


SECTION  V 


FREQUENCY  RESPONSE  OP  A DIGITALLY  CCNTROLLED 
CCNT  HWCUS  SYSTEM 


A.  INTRODUCTION 

When  a continuous  (stable)  linear  system  is  excited  by  a sine  wave, 
the  steady-state  waveform  is  comprised  of  a single  wave  at  the  same  fre- 
quency as  the  input.  It  differs  from  the  input  wave  only  by  a phase  angle 
and  a magnitude  factor.  Moreover,  it  is  unnecessary  to  compute  the  actual 
transient  response  of  the  system  when  the  behavior  for  large  values  of  time 
is  of  interest,  since  both  the  magnitude  factor  and  phase  angle  can  be  read 
from  a Bode  plot. 

As  we  shall  see,  a similar  but  more  complex  situation  pertains  to  a 
discretely  excited  system.  Given  that  the  system  is  stable,  the  continuous 
output  waveform  will  contain  more  than  just  a wave  at  the  fundamental  fre- 
quency; it  will  consist  of  the  fundamental  and  all  of  its  aliases.  Thus,  if 
the  system  is  forced  with  1 sin  bt,  0 < b < 2n/T;  the  output  will  contain 
terms  at  frequencies  b,  b+(2n/T),  b +(Un/T),  ....  The  relative  amplitudes 
and  phase  angles  will  depend  on  the  data  hold  employed  as  well  as  the  sys- 
tem transfer  function.  Nevertheless,  given  a data  hold  and  transfer  function, 
the  magnitude  and  phase  angle  of  each  and  every  component  can  be  read  from 
a particular  "Bode  plot."  Note  that  this  concept  of  frequency  response  is 
more  comprehensive  than  the  traditional  concept  of  the  "sampled  spectrum," 
which  is  limited  to  determining  the  single  sinusoid  that  fits  the  system 
output  samples  at  the  sampling  instants. 

In  the  subsections  to  follow  we  will  first  review  the  Bode  plot  con- 
cept for  a continuous  system  and  then  proceed  to  the  frequency  response  of 
a discretely  excited  open-loop  system.  We  will  conclude  with  demonstra- 
tions of  the  manner  in  which  closed-loop  systems  (both  single  and  multi- 
rate) are  analyzed. 
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B.  CCHTIRU0U8  SYSTEM  BODE  PLOTS 


It  will  be  helpful  to  first  review  the  Bode  plot  concept  for  continu- 
ous systems.  Let  R in  Fig.  58  be  a unit  input  sine  wave  with  frequency 
o^-,  rad/sec.  The  output,  in  the  frequency  domain,  is: 

C(s)  = G(s)  R(s)  = G(s)  (26 7) 


Figure  58 • Continuous  System 


Equation  267  can  be  expanded  in  partial  fractions  as: 


C(s) 


Bs 

2 2 + 2 2 

+ o£ 


Terms  associated 
with  characteristic 
polynomial  of  G(s) 


(268) 


Given  that  all  poles  of  G(s)  are  in  the  left  half  plane,  the  bracketed 
term  in  Eq.  268  represents  time  functions  that  vanish  as  t -•-<».  Thus, 
the  steady-state  behavior  is  completely  determined  by  the  partial  fraction 
coefficients  A and  B,  since  once  they  are  known  the  steady-state  time 
response  can  be  written  directly  as; 


c(t) 


lim 
t 00 


A sin  ob>t  + B cos  o^jt 


^/a?  + B^  sin  (ufct  + cp) 


where  9 = tan""1  (B/A). 


(&9) 
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The  details  of  solving  for  A and  B clearly  show  the  relationship  between 
the  Bode  plot  and  the  steady-state  waveform.  To  solve  for  A and  B,  one 
multiplies  Eq.  268  by  [s*  + o^']  and  evaluates  the  result  for  s = jo^: 


g(sHs=j 


Terms  associated 

(Aui©  + Bs)l  + with  characteristic  (s2+u£)| 

,s~Ja*>  polynomial  of  G(s)  I s J“o 


°(8>|«=W*  = *9*  + = 


A2  + B2  tan  B/A 


(270) 


(271) 


To  summarize,  we  see  that  a sinusoidal  input  at  frequency  o^j  produces  a 
steady-state  waveform  of  the  same  frequency.  It  differs  from  the  input 
only  by  a magnitude  factor  and  a phase  shift.  Both  the  magnitude  factor 
and  phase  shift  for  any  given  input  frequency,  o^,,  can  be  read  directly 
from  the  Bode  plot  for  G(jca).  That  is,  for  any  given  input  frequency  u*>, 


A + JB 


G(s>  s_. 


(272) 


In  the  subsections  that  follow  we  will  expand  this  "frequency  response" 
viewpoint  to  encompass  discretely  excited  systems. 

C.  MATHEMATICAL  PRELIMINARIES 

Let  R be  a sinuqoid  of  unit  amplitude  with  frequency  b rad/sec.  If  R 
is  sanpled  at  l/T  samples/sec  and  then  described  in  terms  of  a N/T  samples/ 
sec  model,  we  obtain 


zw  sin  bT 


z™  — (p  cos  bT)z^  + 1 


= sT/N 


(273) 


For  later  use,  it  is  necessary  to  find  the  N factors  of  the  denominator  of 
Eq.  273  In  a form  which  will  permit  a partial  fraction  expansion  containing 
terms  for  which  corresponding  time  functions  are  known.  For  example,  if 
f(t)  = sin  bt,  then 
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fT(t)  = [sin  bt]T  F(z) 


z sin  bT 

z2  - 2 cos  bTz  + 1 


07*0 


but  we  do  not  know  the  time  function  corresponding  to 


F(z) 


z sin  bT 

z2  + 2 cos  bTz  + 1 


(■£75) 


which  will  also  occur  among  the  N factors  of  the  denominator  of  Eq.  273. 
This  situation  is  remedied  by  adding  n to  the  urgument  so  that  F(z)  in 
Eq.  27*s  becomes  recognizable  in  the  time  domain.  For  example,  Eq.  273 
becomes : 


F(z) 


-z  sin  fb  + (it/T)  IT 
z-'  — 2 cos  [b  + (n/T)lTz  + 1 


=»  f(t) 


(276) 


The  reasons  for  adding  it,  when  obviously  kn  (k  is  an  integer)  would  also 
work,  will  become  apparent  at  a later  point.  Consider  the  denominator  of 
Eq.  273  for  the  special  case  of  N = 2: 


5^  — 2 cos  bTz2  + 1 


bT/2: 


z — 2 cos  bTz*-  + 


= |z*  - 2 cos  ££  z + 1 j |z‘ 


+ 2 cos 


bT 

— 2+1] 

ini) 

: jt  to  the 

argument 

( /,  2*\ 

| (b  * -) 

t|  ) 

078) 


Thus,  for  N = 2,  the  partial  fraction  expansion  of  the  output  will  require 
terms  which  account  for  not  only  the  sine  and  cosine  terms  of  the  input  fre- 
quency b but  also  for  the  first  alias,  b + (2n/T),  of  the  frequency  b as 
well. 

Now  repeat  this  exercise  once  more,  for  N = k.  After  this  the  result 
for  arbitrary  N will  be  apparent  by  induction. 


10:' 


■PPBM 


IWPP 


' 1 
I 


("79) 


* — 2 cos  bz^  + 1 = — 2 cos  — • z2  + ij^  z1*  + 2 cos  ^ z2  + ij 

= (z4-2  cos  ^ z2  + ^[zk-2  cos|^b+-|^  -|-j  z2  + lj  (?80) 


In  turn,  write 


k bT  ? 

z4  — 2 cos  ■—  z^  + 1 


and 


4 ~ bT  2 , 
z +2  cos  — z + 1 


|z2  — 2 cos  z + 1 j |z2  + 2 cos  ^jp  z + l| 

|z2-  2 cos  ^p  z + 1 j |z2-  2 cos  ||b  + -J-J  z \ 1 j ( 281 ) 


= z — 2 cos 


z2  + 1 


K>  * *)  t] 

(z2~2  cos  |(b  + ^pj  -Jp]z  +l)(z2  + 2 cos  |(b  + -J-|z  + l) 

(z2-2  cos||b+~j  x]Z  + 1)(z2“COS[(b+T)  T“]Z  + 1) 


Substituting  Eqs.  281  and  282  into  Eq.  ?80  gives: 


(282) 


z®  - 2 cos 


bTzU  + 1 = 1^  (z2-2  cos[(b+^  t]Z  + 1) 


(283) 


From  Eq.  283,  the  result  for  arbitrary  N is  apparent  by  induction: 
N-1 

n 

n=0 


z2^— 2 cos  bTzN  + 1 


l 

_(*—[(^¥)f]f+(*i»  [(»+¥)  fjr 


N-1 

n 

nO 


(284) 


(285) 
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Equation  Pah  is  true  for  arbitrary  N even  though  the  development  was 
in  te*»s  of  even  N. 

Thus,  the  b rad/sec  sine  wave,  which  is  actually  sailed  on  a T inter- 
val, can  be  m&then&t ically  described  in  terms  of  a function  having  sinu- 
soidal coaq>onents  at  frequency  b rad/sec  and  its  first  N— 1 aliases  (using 
a T/lf  interval). 


Consider  the  system  of  Fig.  59  where  G(s)  represents  an  arbitrary 
transfer  function  and  M represents  an  arbitrary  data  hold.  Suppose  R is 
a unit  amplitude  sine  wave  and  we  wish  to  sample  the  output  with  a T/N 

frame  time. 

CT^  = (GM)T/NRT  = (GM)T/N  si-S-bT , z = esT/N  ( 

ZPN  _ 2 cos  bTzN  + 1 

Expand  the  right-hand  side  of  Eq.  286  in  partial  fractions; 

T/N  tJ  Anz  3in  ^(T/N)  + Bnztz  “ cos  ‘^(T/N)] 

c = 2^  • ^ 1 ■ ..... 

n=0  z - 2 cos  a)n(T/N)z  + 1 

+ [Terms  due  to  modes  of  (GM)1^/^]  (; 

Assume  that  responses  in  the  modes  of  (GM)^/1*  approach  zero  as  t ♦ i 
that  all  modes  are  stable.  In  Eq.  287, 


Figure  59.  Open-Loop  Case 


F 


tUn  = b 4 , n = 0,  1,  2,  N-1 

For  the  present,  we  assume  that  b < 2it/T.  The  steady-state  waveform,  at 
the  sampling  instants,  can  be  written  as: 


C(t)  ^ - V(  (An  sin  ci^jt  + Bp  cos  ayjt) 

n=0 


(288) 


To  solve  for  An  and  Bn,  multiply  each  side  of  Eq.  287  by 


[z°  - 2 cos  (^(T/Njz  + 1 


Oiks  (N-  1 ) 


and  evaluate  for  z = lzko^(T/N).  Since  the  only  term  which  can  survive  on 
the  right-hand  side  occurs  when  n = k,  the  Z will,  of  course,  disappear. 
Then,  if  we  so  choose,  the  k notation  can  be  changed  back  to  n.  To  illus- 
trate, 


«jm)t/n  -m  a'nffw  V 2 cos  “\(T/N)*  + 1 J | 

cos  bTi"  4 1 ,, 


UWylT/H) 


. V C*n‘  [,2-j  co. 

n*0  *2-2  cos  tu^T/Jr)*  4 t s»»*iy(T/N) 


(289) 


For  any  n ft  k,  the  right-hand  side  of  Eq.  289  is  identically  zero  since 


z2  - 2 cos  cok( T/N)z  + 1 = [z  - cos  a^(T/N))2  4 [sin  a^(T/N)12  ( 290) 


vanishes  when 


z 1.xl^(T/N)  cos  a^(T/N)  4 j sin  co^T/N) 


(291) 


t 


■umBfiifa 


Specifically,  we  obtain 


(cos  o^(T/N)  + j sin  oj^T/N)  - cos  aij^T/N)]2  + (sin  gdj^T/N)]2  = 0 (292) 

For  n = k,  the  cancellation  of  the  common  factor  guarantees  the  survival 
of  an  n k term.  Factoring  out  a common  z gives 


z | Ak  sin  u^T/ft)  + Bk(cos  ^(T/N)  - cos  ^(T/N)  + j sin  o^(t/N)]J 


= (Ak  + jBk)z  sin  (^(T/N)! 


I z=1^o)k(T/N)  (29^) 


Therefore,  Eq.  289  becomes 


z sin  u>k(T/N)(Ak  + JBk)  = 


{0/M)T/W(zW  sin  bT)[z2-2  coa  uk(T/H)z  + J ] 


— 2 coa  bTzw  + 1 


|z=ltok(T/N) 


(?9M 


At  this  point,  let  k revert  to  n. 


An+JBn  = (OM)1 


zN~’  sin  bT  z ~2  cos  <Uh(T/N)z  + ) 


|z=14un(T/ff)  sin  u^TT/N)  z2^  — 2 cos  bTz^  + I 


Iz'lku^d/N) 


(‘35) 


The  last  term  on  the  right-hand  side  of  Eq.  295  is  indeterminate  (0/0)  when 
z 1.'So.Nn(T/N ) . Therefore,  apply  L'Hopital's  rule  once  and  obtain 


An  + JBn  = (CM)1 


1 • ZN'’  ‘ [Z~C°8  ai"(T/N)1  sin  bT 

Iz^l&i^T/N)  2NZ1*" ' ( z"  - cos  bT)  sin  uy.(T/N)  ( - 

z=l*uh(T/N) 


An+JBn  = T (0M)T/N 


z=1htun(T/N) 


[cos  qfr(T/N)  + J sin  »'n(T/N)-C03  a^,(T/W)]  sin  bT 

(cos  u>nT+J  r.ln  a^T-cos  bT)  sin  ciy^T/N) 


('97) 


A direct  substitution  for  -■  b + Pim/T  quickly  shows  that  the  last  part 
of  the  product  in  Eq.  097  is  unity.  Therefore, 


T/N 


A„  ♦ T (GM)  ' 


2?e3T/N 


z tVo^T/N) 


(298) 


The  superscript  notation  in  Eq.  298  is  for  the  purpose  of  calling  out  the 
definition  of  z being  used  in  the  evaluation. 

To  review  the  situation,  the  system  is  forced  with  the  sin  bt ; then  the 
steady-state  output  waveform,  sampled  with  a T/N  frametime,  has  the  form 


where 


cT/“(t) 


“'ll 


)T  'N 


b + 


Onu 


n 0,  1 , 0,  . . . , N- 1 


( ■'99) 


{*>00) 


The  coefficients  An  and  Bn  are  computed  using  Eq.  .x)8. 
For  example , let 


so  that 


M 


1 _ e 


G(s) 


T «-™>T  'N 


1 - 


N(r,  - 


■nT/N 
-aT  ;N 


i a 


-N> 


) ( 1 - z-’  ) 


OOi) 


CHV) 


it.  is  instructive  to  plot,  the  Bode  plot,  for  Eq.  ?02  using  N as  a para- 
meter. For  the  sake  of  clarity,  we  will  plot  versus  <0  rather  than  log  ,0 
and,  in  addition,  will  omit  the  phase  angle  plot.  Also,  for  reasons  of 
clarity,  the  ordinate  scales  will  be  displaced  for  the  different  values 
of  N (refer  to  Fig.  Ho).  Note  that.,  over  the  plotted  range  of  8m,  the 
N ) ease  repeats  itself  H times.  Tn  a like  manner,  the  N 
repeats  twice,  whereas  the  N U goes  through  one  cycle. 


. ' case 


Figure  1*0.  Magnitude  Flot-  for  N = 1 , ?,  4 


Imagine  that  a unit  amplitude  sine  wave  is  input  with  a frequency  of 
n/f.  In  the  N 1 case,  our  only  interest  is  matching  the  sampling  points 
with  a single  sine  wave.  The  magnitude  and  phase  angle  (not  shown  in 
Fig.  4o)  could  be  read  from  this  plot  at  a'  = it/?,  it/2  + 2n,  n/2  + l*it, 

*/?  + on.  . ..;  each  point  gives  the  correct  values.  Assume  next  that  the 
input  has  a frequency  b = {n'2  + Uir).  Clearly,  if  the  objective  is  to 
match  the  sample  points  with  a single  sinusoid,  the  frequency  of  the  out- 
put could  be  b plus  any  ,'ir/T  multiple;  the  sampler  cannot  tell  the  differ- 
ence. In  fact,  the  "sub"  aliases  at  b - 'n/T,  b-l*it/T  will  also  work.  These 
"sub"  aliases  are  the  difference  terms  so  prominent  in  modulation  theory. 

Figure  4l  depicts  this  situation  for  a steady-state  response  given  an 
input  frequency  of  b = n/C  rad^see  and  assuming  the  system  used  in  com- 
puting the  Bode  plot  of  Fig.  !k\  For  the  sake  of  clarity,  only  two  of  the 
many  waves  which  fit  the  sample  points  are  shown  — one  at  it/?  rad/sec.  the 
other  at  [(it  '.  ) - ( T)]  rad/sec. 


108 


O Sample  Points 

— — A sin[(r/2  )]t  ♦ 8 cos  [(r/2)]» 

A sin[(W2)-(2r/T)]t  + B cos[(ir/2)-(2»-/T)]t 


Figure  '‘I.  'Two  Continuous  Sine  Waves  Which 
Match  the  Sample  Points 

The  N 1 plot  in  Fig.  ••0  corresponds  to  the  "sampled  spectrum"  of 
sampled  data  control  theory.  Turn  now  to  the  N = 2 case  wherein  the  objec- 
tive is  to  match  one  inter-sample  point  as  well  as  the  sample  points.  Let 
the  input  frequency  be  jt/f  and  note  that  the  points  at  u>  = it/2,  it/2  + fit 
give  the  correct  answers,  as  would  the  points  jt/2  + Un,  it/ 2 + tore.  Suppose 
next  that  the  input  frequency  is  b it '2  + 2it.  Clearly,  the  second  required 
component  could  be  read  from  the  "first  alias"  at  b + Cit/T  or  the  first  sub- 
alias  at  b — . it 'T  (or,  for  that  matter,  a whole  host  of  other  frequencies). 

In  the  N •*  case,  four  sine  waves  are  required  to  fit  three  inter- 
sample points  as  well  as  the  sampiled  points.  If  the  input  frequency  were 
b it  ♦ i'it,  and  if  the  plot  of  Fig.  ;*0  with  its  limited  range  of  8n  were 


>9 


the  only  one  available,  clearly  it  would  be  to  our  advantage  to  use  the 
"difference”  frequency  points  at  o>  = b - 2rt/T,  b — 4jt/T,  and  b - 6jt/T  to 
establish  the  magnitude  and  relative  phase  of  the  three  remaining  sine 
waves . 

This  brief  discussion  serves  to  point  out  that  the  aliases  and  sub- 
aliases  can  be  associated  with  the  sum  and  difference  frequencies  of  modu- 
lation theory.  One  should  not,  however,  think  that  both  sum  and  difference 
components  must  be  simultaneously  present  in  the  output.  Clearly,  only 
N components  are  needed.  We  can  now  remove  the  earlier  constraint  that 
b < 2jt/T.  If  b is  less  than  2jt/T,  it  is  certainly  true  that 

= b + , n = 0,  1,  2,  ...,  N-1  (303 


However,  if  b > 2jt/T,  let  o>s  = 2n/T  and  use 


to  restate  Eq.  303  as 

= b + ^ > n=no>  no  + ^>  • ■ 0,  1,  2,  N—n0— 1 

For  our  previous  example  where  N = 4,  suppose  b = 6n  + ?t/2.  Then 


n 


o 


l6n  + (rt/2)\ 

\ 2lT  /INT 


-3 


Thus, 


“>n 


b 


n = -3,  -2,  -1,  0 


(304 


(305 


(306 


(307 


and  we  use  three  sub  aliases.  If  b = it/2,  then 


n = 0,  1,  2,  3 


and  we  use 


u>ri  = 


’nil 

T • 


the  "positive"  aliases. 


Keep  in  mind  that  all  this  represents  a convention  which  the  reader 
may  not  necessarily  elect  to  follow.  What  is  important  is  a clear  under- 
standing which  will  permit  one  to  pick  a consistent  set  of  N points  from 
the  Bode  plot. 


Of  interest  is  the  case  where  N is  extremely  large.  In  fact,  let 
N -»*  >«  after  evaluating  Eq.  ’'OS  at  z = ^^(T/N): 


1 _ e-»T/N  , _ „-N  (l  - e-nT/r* ) ( i - 14-o^T) 

N(z  - e_Kl’/N)  ) - z-1  N[l&q,(T/N)  _ e“*T/N][l  - 14-.«n(T/N)J 

z-l^(T/N) 

An  indeterminate  form  is  obtained.  Therefore,  use  L'Hopital's  rule  twice 
(substitute  14uy,(T/N)  = cos  uy^T/N)  + j sin  oy^T/N),  etc.)  and  obtain  (see 
Appendix  l) : 


(508) 

I 1 im 
|N-*-oo 


1 

TT 


T/N 

z = 1JVi^j(T/N  )| 


11m 
N—  * 


1 _ e-J<‘ViT  1 

JaJ  ' + J'On 


I 

8 + I 


(509) 


That  is,  as  N -»•  «>,  one  simply  divides  GM  by  T and  evaluates  the  coeffi- 
cients at  s jay,.  This  is  representative  of  the  general  result  discussed 
in  the  next  subsection. 

E.  OPEN  LOOP  FREQUENCY  RESPONSE  — 

CONTINUOUS  OUTPUT 


In  the  previous  section  it  was  shown  that 


i <«T/N 


z e 


sT/N 


1 1 


^n  *■  J^n 


a I (T/N) 


(^10) 


To  deduce  the  behavior  for  infinite  N,  rewrite  Eq.  310  as 


The  T/N's  cancel,  and  as  N gets  very  large  only  the  k = 0 term  contributes 
since  all  the  "aliased"  spectra  have  moved  to  infinite  frequency.  There- 
fore, 


where 


lim 

N -*-«> 


1 

T 


(GM)T/N 


z4esT/N 

z=14c%(T/N) 


‘lln 


b + 


Pmi 

~T~ 


n = n, 


o> 


no  + ^ > 


(310) 


The  finite  N example  of  the  previous  section  can  now  be  studied  for  the 
case  oi  ini inite  N — therefore  the  "continuous"  frequency  response.  Thus 


A„  + JBn 


1 

s + 1 


s=jccn 


(313) 


This  is  shown  in  Fig.  Up  where  the  components  for  input  frequencies  of 
h 0.1,  1.0,  U.O,  o.O  rad/sec  have  been  indicated  with  different  symbols. 

The  interpretation  of  Fig.  UP  is  as  follows.  Suppose  a unit  sine  wave 
at.  I rad /sec  is  input  to  the  sampler.  Then,  if  sine  waves  at  1 , i + Pn/T. 

I *■  Uir/T,  ...,  are  added  together,  the  resultant  waveform  will  be  an  exact 
match  of  the  actual  steady-state  output  waveform.  One  would  expect  this 
waveform  to  be  relatively  clean,  since  the  first  alias  is  down  by  something 
m the  order  of  <0  dB,  relative  to  the  input  component.  However,  the 


Figure  4?.  Frequency  Response  and  Spectral 
Components  of  Output 


transient  response  itself  does  not  bear  out  this  conjecture  as  can  be  seen 
in  Fig.  4 ■'>.  The  reason  is  that,  the  higher  terms  are  important.  They  do 
not  represent  "harmonic"  terms  but  are  rather  modulation  components  which 
must  add  together  properly  in  order  to  match  conditions  at  the  T transi- 
tion points.  As  mentioned  earlier  in  Section  IV,  it  can  be  seen  that  the 
"steady  state"  does  not  necessarily  take  on  the  additional  attribute  of 
periodicity.  This  occurs  only  when  the  input  frequency  and  the  S'impling 
frequency  bear  an  integer  relationship  with  respect  to  one  another.  As 
another  example,  the  "frequency  response"  for  a higher-order  system  is 
shown  in  Fig.  44.  In  this  example,  the  sampling  frequency  was  deliber- 
ately picked  to  suppress  the  bending  mode  responses. 


( 

i 


P.  INPUT  SIGNAL  WITH  PHASE  SHIFT 


If  the  input  signal  has  the  form 


r(t)  = k1  sin  bt  + kg  cos  bt 


(31*0 


the  results  of  the  previous  section  are  changed  only  by  a complex  constant. 
Following  exactly  the  procedures  of  Subsection  D,  except  for  using  the  more 
general  input  given  in  Eq.  311*,  gives  (see  Appendix  E): 


An  + j Bp 


z=esT/N 

z=14t%(T/N ) 


(k-|  + jko) 


(315) 


Limit  N — °° 


A-n  + 


(k]  + Jk2) 


(3 16) 


0.  SINGLE -RATE  CLOSED-LOOP  FREQUENCY  RESPONSE 

The  closed-loop  results  will  be  conf iguration  dependent.  However,  the 
mathematics  are  quite  tractable  and  can  be  followed  through  on  a case-by- 
case basis.  Thus  it  is  more  important,  for  the  present,  to  have  an  insight 
into  the  mathematical  structure  and  the  particular  simplifications  that  sur 
face  in  a closed-loop  analysis. 

Consider  the  (vector)  system  shown  in  Fig.  **5.  The  procedure  we  now 
follow  will  be  typical.  First,  solve  for  the  vector  component  at  the  input 
of  the  data  holds . 


ET  . <,y-GTo*(OK)V 

Therefore 

ET  = [I  + G^G^GM)1]”1  Gf  RT 

Next,  solve  for  C(s): 

C = (GM)[I  + G'jGT,(GM)Tr1  Gf  RT 
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(317) 

(318) 

(319) 


; 


Figure  4*?.  Illustrative  Vector  Closed-Loop  Configuration 


The  spectrum  of  C(s)  is  what  is  really  of  interest;  we  seek  it  by  finding 
first  the  spectrum  of  C and  going  to  the  limiting  case  of  N -»■<». 

Let  the  input  be  a sine  wave  at  frequency  b rad/sec  and  let  the  delay 
operator  be 


-1  -sT/N 

z = e ' 


(320) 


so  that 


z"  sin  bT 

■ 2 cos  bTz^  + 1 


(321) 


Therefore 


CT/N  = (GM)T/N  [I  + G^(GM)T]“1  G^  Ra 


(322) 


For  the  sake  of  brevity  write  Eq.  322  as 


CT/M  „ aT/»GTRT 


Expand  the  right-hand  side  of  Eq.  323  in  partial  fractions: 


(V3) 


T/N  _ pT/N  T 

~ GA  Gb  2N  , N , 

z - 2 cos  biz  + 1 


Anz  sin  <rn(T/N)  * Bnzfz  - cos  ^(T/N)] 
n^o  - 2 cos  a>^(  T/N)z  + 1 


t [Terms  due  to  modes  of  gT 

A D 


(324) 


Assume  that  responses  in  the  modes  of  G*/NG*  approach  zero  as  t 
that  all  modes  are  stable.  In  Eq. 


b + 2nn  ^ n = 0,  1,  2,  N-1 


[More  generally,  n = nQ,  nQ  + 1 , N-n0-1,  where  nG  = -(b/c*)^] • 

Note  that  Eq.  32^  is  exactly  the  same  as  Eq.  287,  except  (GM)  / has 
been  replaced  by  g£/NG*  . Hence  we  can  write  a crucial  result  using 

Eq.  298: 


1 rT/NpT 
An  + OBn  = -jf  Ga  GB 


g|  * gbC=") 


z^^T/N) 


(326) 


(327) 


Therefore,  using 


[l4mn(T/N)]N  = I^T 


= cos  <jonT  + j sin  c%T 

= cos  [b  + (2nn/T)]T 

+ j sin  [b  + (2im/T)]T 

= cos  bT  + j sin  bT 


we  obtain 


gZ  = G[l4o)n(T/N)f 

B l^d/N) 


3 G[l4<JJnT]  = G[UbT] 


(328) 


Thus,  we  can  replace  zR  in  g|  with  z and  evaluate  it  at  z = UbT.  At  this 
point  we  have 


I 


Equation  Vo  is  the  basic  result  for  the  finite  N ease.  To  reiterate,  to 
find  the  coefficients  of  the  N sine  waves  for  the  T/N  sampled  output  of  C, 
we  compute  the  normal  "T"  transfer  functions  for 

g'J  fl  + G? Gp(GM)Tr'  af  (*^0) 

and  evaluate  t.  for 


l.xbT 


( AM  i 


Next,  compute  the  normal  T 'N  pulsed  transfer  function  for  G^ ; evaluate  Lt 
at  :t  l.\ip(T  N)  where  up  b t (.  irn  T). 


Thus , 
it  suff i 
s imp  ! i fit 
the  ease 


T • 

Gp  ut  periodic 
es  to  use  bT. 

the  proeedure 
of  N — »■»,  the 


in  (.'mi  Ti  and  it  is  superfluous  to  use  uy,(T/N); 
T 'N 

Moreover,  only  the  G.  is  a function  of  N;  tiiis 
involved  in  the  limiting  case  tremendously.  For 
continuous  ease,  we  obtain: 


+ JBn 


M( s ) G ( s ) 

T 


e‘"n 
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Equation  A is  the  desired  result  for  the  given  closed-loop  configura- 
tion.* However,  the  mathematical  ideas  are  what  count;  one  can  follow  the 
details:  til  rough  for  other  configurations  quite  easily. 


* T i 

Accurate  numerical  determination  of  Gp  I ixbT  m;i-v  prove  difficult  at 

high  sampling  rates.  This  is  the  result  of  "small  differences  between  large 
numbers  which  occur  in  the  computations  as  poles  and  zeros  approach  the  unit 
circle.  In  this  event,  one  is  well  advised  to  carry  out  equivalent  compu- 
tations in  a domain  where  numerical  conditioning  is  much  improved  (e.g.,  in 
terms  of  w'  or  wi. 


K.  THE  ''MOTIVATING"  EXAMPLE  REVISITED 


Recall  the  single-rate  design  problem  of  Section  IV  wherein  a closed- 
loop  system  was  designed  to  have  exactly  the  same  step  transient  response 
at  the  sampling  instants  even  though  the  open-loop  parameter  took  on  two 
different  values.  For  convenience,  the  step  response  is  repeated  in 
Fig.  ho. 

Using  the  result  of  the  previous  subsection,  we  iire  now  in  a position 
to  examine  the  frequency  response.  The  magnitude  plot  is  given  in  Fig.  h'j 
(the  magnitude  is  plotted  against  u>  rather  than  log  a'). 


Figure  -o.  Closed-Loop  Step  Response 

Note  that  the  "sampled  spectrum,"  the  N 1 case,  is  unable  to  differ- 
entiate as  to  whether  the  open-loop  parameter  is  either  a — or  a - 1C. 
Moreover,  the  difference  in  the  frequency  content  of  the  steady  state  when 
the  ZOH  is  used  cannot  be  distinguished  from  that  when  the  slower  is  used. 
Certainly  (at  least  for  this  example)  the  classical  "sampled  spectrum"  is 
of  limited  vise.  On  the  other  hand,  the  continuous  spectrum  (N  case) 

clearly  displays  the  differences  in  the  closed-loop  designs.  Note,  in 
particular,  the  added  width  of  the  notches  when  the  slower  is  used;  the 
COH  notches  are  very  sharp  by  comparison.  Thus,  one  would  oxpcet  the 
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Before  leaving  this  example  ’et  ms  talte  the  "rportunitv  to  solidify 
the  meaning  of  the  frequency  re*  vmse  for  :.•»  fli.it>-  j|  cuse.  Consiaer 
•he  a = 10  case  and  use  the  recurs  lot  «n;n  .t  . to  nen  t-  th,  cont  l:.uou 
transient  response  for  a uni‘  anpl  • • nde  i-  wav,  ..  t.h  a frequency  of  n 
This  is  shown  in  Fig.  U8.  According  *o  theory  wc  should  he  able  to  vt 
N 1 and  from  the  Bode  plot  read  the  magnitude  anl  phase  of  the  single 
sinusoid  that  fits  the  sample  points  at  th ' 'amp's  *nrhant.  This  indeed 
♦'roves  to  he  the  ease  and  is  shown  in  F!,-.  ho.  pi^  !.o  :x  section  f 
*he  transient  response  has  been  "eep  ed"  uni  . vert  .i  with  the  sine  wave 
whien  results  from  the  N 1 amputation,  namely 


(Clt)l 

\T 

/ 

\ 

t~) 

;•  in  — t 

f bo  o 

' ~t) 

Equation  indicates  that  the  value  of  the  wave  at  the  sampling  instants 

is  of  interest.  However,  for  expository  purposes,  we  have  shown  a complete 
cycle  in  the  figure. 

Next,  consider  the  (N  ■»  P)  case  wherein  the  desire  is  to  fit  not  only 
the  sample  point  hut  one  inter-sample  point  as  well.  This  case  is  shown 
in  Fig.  [-o.  The  T/P  response  equation  is 


c<1’  ss 


a0  sin  ~ + b0  ('os  — t 


+ a,  sin 


(-r-  + -rrl  t + b,  cc 


(^  + t) 


J 


Figure  '0.  Ste&dy- state  Sinusoidal  Components  n 


*HF  X 


Again,  for  expository  reasons  we  have  shown  the  continuous  waveform  which 

results  from  sines  and  cosines  at  co  = rt/2  and  its  first  alias  at 

a)  = ( it/2 ) + (?jt/T).  A half  period  for  the  N = 10  case  is  shown  in  Fig.  1. 

The  reader  is  reminded  that  the  steady-state  wave  of  this  example  is 
periodic  and  free  of  modulation  effects  simply  because  the  selected  input 
frequency  bears  an  integer  relationship  to  the  sampling  frequency.  In 
general,  the  expansion  of  the  steady  state  in  terms  of  aliases  is  a non- 
orthogonal  series.  For  those  special  cases  where  the  steady  state  is 
periodic,  it  is  readily  apparent  that  the  N = °»  case  degenerates  to  a 
Fourier  series. 


Figure  ‘ 1 . Steady-State  Sinusoidal  Components,  N = 10,  Half  Period 
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a0  0.002998683 
b0  -0.302606086 


a-|  -0.01+6197029 
b-,  -0.01+354801 8 


a0  -0.174468021  a0  -0.048579760 
b0  -0.287649137  b0  -O.3O6381976 

a1  -0.125888261 
b!  0.018732839 


a2  -0.051578443 
b2  -0.003775889 

a^  -0.079691232 
bj  0.062280857 


I.  A PARTICULAR  MULTI-RATE  CONFIGURATION 

As  with  the  analysis  of  closed-loop  single-rate  systems,  the  analysis 
of  the  multi-rate  closed-loop  case  is  configuration  dependent.  However, 
the  mathematics  remain  quite  tractable,  as  we  shall  attempt  to  demonstrate 
in  this  and  the  succeeding  subsections.  Here  we  consider  a particularly 
simple  multi -rate  configuration;  the  next  subsection  will  treat  a more 
complex  configuration. 

Consider  the  two-rate  system  shown  in  Fig.  52.  In  Fig.  52,  Wi  and  Wp 
are  compensation  networks,  M is  a data  hold,  and  G represents  the  open- 
loop  system  dynamics.  One  may  consider  these  to  be  matrices  of  the  proper 
dimensions . 

The  objective  is  to  find  the  "frequency  response"  for  the  output 
vector  C.  As  in  the  single-rate  case,  we  shall  assume  that  C undergoes 

124 

f 1 


TABLE  7-  COMPONENT  COEFFICIENTS  N = 1 , 2,  4 


It  is  also  important  to  bear  in  mind  that  we  are  "matching"  the  sample 
points;  thus  the  a 0 term  for  N = 1 will  be  unequal  to  the  aQ  term  for 
N = 2.  This  is  demonstrated  in  Table  7 for  N = 1,  2,  and  4 (let  b = jt/2) . 


Figure  A Specific  Two-Rate  Closed-Loop  Con  figuration 


a phantom  T/N  sampling  operation  and  then  seek  the  limit  as  N -»<w.  Fr 
Fig.  ‘d'  it  is  seen  that 


or 


c gmet  m 


CT*  (OM)T/NET/M  , N/M  an  integei 


Thus,  the  first  task  is  to  solve  for  ET/M  . This  is  a non-trivial  task 
the  details  must  be  followed  with  care. 


E W,RT  - W]W?,(GMET/M)T 

Therefore, 

ET/M  w'f/M  RT  - w'[/Mw'r.(C.MET/M)T 


Pre-multiply  Eq.  M8  by  .GM  and  sample  at  a T interval. 

(GMEt/m)T  (GwT/m)TRt  - (GMWT/M)TWT(GMET/M)T 

' c* 


Solve  Eq.  ^9  for  (GMEt/m)T: 

(GMEt/m)T  [I  + (0NWT/M)TWT]“'  (GNW^/Mf  RT 


Substitute  Eq.  ^tO  into  Eq . V^8  and  elear  through.  The  result  is; 


ET/M  WT/Mjl  - WT[I  + ( GMVT/m)T wT ]— ' (GMW^/m)T!  Rt 


( *'1»  I ) 


For  brevi ty , let 


et/m  „T/M0TrT 


('ll.) 


The  evaluation  of  oj[  is  non-tHvial;  for  example,  the  (GNW’jj  1 element  of 
G^  will  have  to  be  eomputed  using  either  switch  decomposition  or  the  phantom 
sampler.  Via  the  phantom  sampler. 


(GWT/M)T  [(GM)T/MWT/Mf 


(^M 


So  far,  the  two-rate  example  yields: 


ct/h  (0„)T/»wt/mqtrt 


Oik) 


and  we  see  that  the  only  new  element  added  over  the  single-rate  case  is  the 
addition  of  a term  sampled  on  a T/M  interval  and  the  additional  constraint 
that.  N/M  be  an  integer. 


(A!,-,) 


so  that  a unit  amplitude  sinusoidal  input  at  b rad/sec  has  the  transform 


c.N  sin  bT 

z->N  - ? cos  bTzN  \ I 


(M«>) 


As  in  the  single-rate  case,  substitute  Eq.  J>h6  into  Eq.  Jhh: 


CT/N  = (GM)T/NWT/MGT  '^sin  bT 

z™  — 2 cos  bTz^  + 1 


(W 


Expand  Eq.  3^7  into  partial  fractions*: 


CTA« 


n=0 


Anz  sin  ayi(T/N)  +Btlz[z-cos  ( T /N  ) ] 

~2  

z -2  cos  ap(T/N)z  + 1 


Terms  due  to 


(318) 


where 


a\. 


( J'!,9) 


(Keep  in  mind  that  we  can  no  longer  pick  N arbitrarily.  Thus,  if  M 3,  we 
cannot  look  at  CT  within  the  present  problem  structure;  however,  we  can  look 
at  CT/\  CT/‘’,  CT/-,  etc.,  since  N/M  is  an  integer.  Thus,  as  we  let  N — - 
we  would  have  to  do  so  in  multiples  of  3.) 

Assume  that  the  terms  due  to  (GM)T/N  w'f/M  gJ  vanish  as  t — ».  Thus, 
the  steady-state  response,  at  the  sampling  instants,  is  given  by: 

N-1 

C(t)|t— " 2Zo  An  sin  aht  + cos  a^t  (3^0) 

To  find  An  and  Bn,  multiply  Eqs.  ?1»7  and  *<h8  by  [z:"  - p Cos  cck(T/N)z  + l ] 
(note  the  dummy  index). 


"At.  this  point  we  may  fall  back  on  the  results  of  previous  subsections 
and  proceed  directly  to  Eq.  358.  However,  it  is  a worthwhile  exercise  to 
complete  the  proof  in  order  to  gain  familiarity  with  the  notation. 


1 7T 


T zN  sin  bT[z"'_  2 cos  <%(T/N)z  + 1 ] 


z — 2 cos  bTz  + 1 


r?  „ /m/„N  ,i  v'1  Anz  sin  o^T/N)  + Bnz[z-cos  a*(T/N)] 

[z^-2  cos  u^T/Njz  + 1 J 2^  ^ 

n=0  zc  —2  cos  (^(T/Njz  + 1 


+ [z2-cos  <jdjJT/N)z  + 1 ][(GM)T'/Nw’J'//MgJ]  (350 


Evaluate  both  sides  for  z = cos  cd^T/N)  + j sin  oj^(T/N).  Thus,  the  only 
term  which  survives  in  the  sunroation  on  the  right-hand  side  is  when  n = k. 
We  may  as  well  drop  the  k notation  and  retain  n (the  choice  is  ours). 

So  far, 


Anz  sin  (Ufi(T/N)  + jBnZ  sin  (^(T/N) 

- ft-M)TA».,TVMcT  =”  Sln  bTtg‘’  ~ COS  ^(T/")--  « 'I 
1 ’ 1 * =''»  - . cos  OTs"  * 1 


z = 14ay,(T/N)  (352) 


/^w\T/N,,T/NUT  2 sin  bT[z‘'-2  cos  ^(T/N)z  + 1 ] 
An  4 (GM)  W.  G.  ™ 

1 A sin  a)n(T/N)(z  -2  cos  bTz.  + 1 ) 


z = 1.^avn(T,/N) 


Note: 


(553) 


1zk0n(T/N)  = cos  a^,(T/N)  + j sin  a^(T/N) 

The  right-hand  side  of  Eq.  353  is  indeterminate  when  z - 1:Smn(T/N). 
Therefore,  use  L'Hopital's  rule, 


L 


■».  . com )***/*; 


2z-2  cos  oyj(T/N) 


2N(z'°N_1 -cos  bTz1*"1) 


(351*) 


140^  (T/N) 


- (GM)^^ 


sin  bT 
sin  cnnC^/W) 


J^y[2 -cos  u)n(T/N)] 
( zN  — cos  bT ) 


^(T/N) 


(355) 


An  + jBn 


(gm)T/nwt/mgt 


^(T/N) 


sin  bT[cos^Bfj(T7N)  + j sin  a)n(T/N)  - cos_i*)nfT7N)  ] 
[sin  ay^T/N J JN(  cos  ay,T  + j sin  uyT-cos  bT) 


Therefore, 


An  + J®n 


i (GMJ^W^G l 


14a^(T/N) 


(sin  bT)[j  jjjw^(f/N)] 

X [sinjiijjiX/ft)  }(cos  oyjT+j  sin  oyjT  — cos  bTj 


(356) 


(357) 


but, 

[l4cun(T/N)]N-cos  bT  •=  14a)nT-cos  bT 

= cos  a)nT  + j sin  a^T  - cos  bT 

= cos  [b  + ( 2?tn/T) ]t  + j sin  [b  + (2m/T)]T-cos  bT 
= cos  bT+j  sin  bT  — cos  bT 
- j sin  bT 

Therefore, 


An  + J®n 


1 

IT 


(GH)T'\T/X 


14a)n(T/N) 


(358) 
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Next , 


ga(zN) 


Therefore, 


ga(-n) 


'^(T/N) 


= ga  [l^j(T/N)  ]N 

Ga(!4(%T)  = GA(l4bT) 


(359) 


That  it,  take  the  "T"  s-transform  of  Ga,  and  evaluate  « 
the  "new"  element,  : 


at  z = UbT.  Now. 


W^/M  t W (2n/m} 


(360) 


Therefore , 


^|l^(T/N)  - W1 fl4ufc(T/N)]N/M 

= V’‘WT/M)] 


(361) 

(362) 


That  is,  take  the  "T/M"  z-transfo™  of  W,  and  evaluate  i 
At  this  point,  only  GM  depends  on  N and 


N 


it  at  z = ^(T/M). 
we  can  go  to  the  limit  of 


An  + JBn  = ± 


(GM)T/n 


z4esT/N 


: = 14u>n(T/N) 


WT/M 


z^esT/M 


z = 14«%(T/M) 


z^esT 


z=14bT 


(363) 


This  is  the  desired  result  for  finite  X (remember  N . M,  a,,  ete.).  As 
’ tl"'  C”mclents  of  ‘he  continuous  spectrum  are  given  by 


An  f JBn  = j-  (GM) 


WT/M 


ztesT/M 


UA 


sT 


z=e 


z 14bT 


(36h) 


no 


An  example  which  illustrates  the  use  of  Eq.  36k  is  given  in  Fig.  55. 

This  frequency  response  (which  is  plotted  on  a linear  scale  rather  than 
log  cc)  can  be  contrasted  with  the  single-rate  results  shown  in  Fig.  3k.  In 
each  case,  the  design  objective  was  to  force  unity  DC  gain  and  a closed-loop 
short  period  at  A = (s  + 2)2 + ( 2)2  when  the  open  loop  was  D = (s+2)2  + (6)2. 

Note  that  the  two-rate  configuration  requires  a significant  increase 
in  the  feedback  gain.  Also,  the  first  "notch"  does  not  occur  until 
cn  = 100  rad/sec. 

J.  A FAST  INNER-LOOP,  SLOW  OUTER-LOOP  PROBLEM 

The  vector  configuration  of  Fig.  55  is  typical  of  the  "fast  sampled" 
inner-loop,  "slow  sampled"  outer-loop  format.  This  type  of  configuration 
is  particularly  easy  to  analyze;  closed-form  analytical  solutions  can  usually 
be  obtained.  This  will  be  demonstrated  using  the  configuration  of  Fig.  55. 


c,  = 

G^E11  - HG^rC^'/M 

(565) 

Therefore 

or 

cf/M  = 

(G1M1  )T//MET-  (HG2M2)T//MCyM 

(566) 

Is 

•11 

[I  + (HG2M2)T/Mr1(G1M17'/M  ET 

(567) 

Let 

rT/M 

ga 

[I  + (HG2M2)T/Mr1(G1M1)T'/M 

(-68) 

so  that 

11 

£ 

tHj- 

A 

gI/met 

(569) 

Next,  solve  for  ET  in  terms  of  RT:  . 

E = R - C 

= R - G^loCf^  = R - G^nG^//MET  (370) 


1?1 


Therefore 

ET  = 

R*  - 

f 

i ! 

j 

and  we  can 

T 

solve  for  E-1 

ll 

1 

eT  = 

[I  + (G^T/M)1]-1^ 

(371)  | 

Since 

c = 

G^C^/W  = Et 

(572) 

we  can  use 

continuous 

Eq.  571  to  obtain  a closed-form  analytical  relationship  for  the 

output  vector  C.  The  result  is 

i 

j 

C = 

GoM^gJ/M[I  + (GoMdGj/M)TI  1RT 

(373)  J 

| 

where  G^/M  is  defined  by  Eq.  568.  Notice  that  either  vector  switch  decom-  j 

position  or  the  phantom  sampler  concept  must  be  used  to  evaluate  . 

J 

. 

1 

Let 

o 

Cbt-3 

II 

[i  + (g^m2gT/M)t]  1 

(37*+) 

t 

so  that 

1 

c = 

Gr>UlJ//MGgRT 

(375) 
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Next  we  use  the  T^N  concept  to  obtain 

cT/N  (g.^Ot/ngT/mgTrt 


(576) 


which  requires  than  N/M  be  an  integer.  For  finite  N the  spectral  coeffi- 
cients  are  givren  by 


An  + J^n 


<1^  <*** 


A B 


lo(apT/N) 


where 


and 


b + 


. mi 


aYi  “ ' t 

r(t)  sin  bt 


(*77) 

(578) 

(579) 


In  the  limit,  as  N -»■  co,  we  obtain  the  continuous  signal  spectrum  coeffi- 
cients 


An  + J Bn 


G ; -M  ' 

rT/M 

°A 

T 

s jap 

z^esT/M 


B 


A'3'1’ 


2= 14bT 


(-80) 


As  usual,  the  notation 


T /M 


z4esT/M 


z 1A(uy,T/M) 

means  compute  the  z-trans  form  of 


s-l.’tohd/M) 

with  respect  to  a T/M  second  frame  time  and  evaluate  it  at  z cos  up(T/M) 
+ j sin  ap(T,/M).  In  a similar  fashion  the  z- transform  of  Gy  is  computed 
with  respect  to  a T frame  time  and  evaluated  at  z cos  bT  + ,i  sin  bT 

To  this  point  in  the  development,  all  results  apply  regardless  of 
whether  Fig.  *'1'  is  a scalar  or  a vector  block  diagram. 

It  is  worthwhile  to  go  through  a simple  scalar  illustrative  example 
in  order  to  bring  the  mathematical  details  sharply  into  focus. 

Example 

Let 


S -f  1 


H = -1 


Mp 


M, 


— e 


-sT  /M 


1 — e 


— sT 


(581) 

(*8f) 


At 


S f 2 


Reviewing  the  equations  it  is  seen  that  we  must  compute: 


0 (hgjOt/m 
O (g,m,)t/m 

')  [I  + (HC.  >MOt  Mf1 

4)  G J M [I  + (HGoMp^/Mf1  (G,M,  )T /M 

s)  (G  M ,gT/m)T 

6)  G'J  [I  t (G  ,M  ,GT/M)Tr' 

Let  us  now  compute  Items  1-6  noting  as  we  go  the  "tricky  points"  that  one 
may  have  a tendency  to  stumble  over. 

Item  1 


( HG  'W  0 


T M 


T 'M 


1 — e 
s ( s + 1 ) 


T M 


This  is  a trivial  computation  since  the  frame  time  of  the  sampling  opera- 
tion and  the  frame  time  of  the  data  hold  involved  are  the  same.  Let 


,sT/M. 


(hg  j*Ot  m 


1 - 

^Fm 


( ^ ) 


Item 


(G,Mf) 


T ;M 


2(  1 - e"67  ) 
s(s  + 2) 


T'M 


We  proceed  a little  more  carefully  with  this  computation  since  the  frame 
time  of  the  sampling  operation  is  shorter  than  the  frame  time  of  the  data 

% t 

hold.  Again  we  use  z e 1 


iv, 


we  obtain 


in  a more  illuminating  form 


The  incentive  for  rewriting  Eq.  387  in  the  form  of  Eq.  388  is  simply  that 
we  recognize  that  this  system  only  has  two  poles;  the  other  poles,  zM_1, 
have  been  forced  on  us  by  the  difference  in  the  frame  time  of  the  sampling 
operation  and  the  frame  time  of  the  computer. 


Next,  recall  that  a function  of  z has  a correspo 
function  which  we  can  visualize  as  being  the  "genera 
function.  For  example,  we  can  say  that  the  function 


was  generated  by  sampling  f(t)  at  intervals  T/M  where 


is  a unit  step  function.  This  insures  that  f(t)  = 0 for  t < 0. 
nd  by  taking  advantage  of  this  observation;  simply 

continuous  time  "generator"  for  (G2M2)t/mgJ/m  and  resample  it  using 
ing  interval!  Thus  we  can  rewrite  Eq.  388: 


I 111'  i\'||(  i lull'll 


• run;:  form. 


mu  i’o  mi 


whore  • is  understood  to  l>o  o"  . If  wo  koo]>  0° 

OttOO,  till'll 


ns  >u  r relVr- 


[(G  ,M1)T/:\;T  .'|T 


21U_ 


1)(  •• 


(I 


-t 

0 


e*-1  ) 


( ?'9(') 


For  frequency  response  purpose's , wo  can  use  Kq.  evaluated  at 

l.vi^T  or  Kq.  '>><■  evaluate  I at  c i.\n,,(T/r);  the  answers  will,  of  course, 
l>o  tlio  same. 

Return  now  to  the  basic  problem  at  hand.  Instead  of  setting  M lot 

it  bo  arbitrary  and  ro-evniuato  Eq.  ''l)  as 


[(G  ,M  ,Vr  MqT  M1 


( I 


( 597) 


To  derive  Kq.  '97 » simply  write  the  generic  term  for  t lie  delayed 
recognizing  that  the  sum  resulting  from  (z  1 t •••  > ^ ) can  be 
closed  form  since  it.  is  a geometric  progression.  This  concludes 
tion  of  Item  ■<;  it  was  the  most  difficult  step. 


.•,-traiis  form 
placed  in 
tin'  oval  un- 


it.em  (> 


fl  i (G.,M  ,C.|  M)  | 


r ,-i 


( *>98) 


HVW 


I 


We  are  now  ready  to  evaluate  Eq.  ?8o. 


An  + jBn  = 


i - 1 

U - 

ZM-1 

1 - e-®/"! 

s*r  s + 1 

s=J(Un 

(z  “ I)2 

•75TT- 

z - 

(z  - l)(z  - e^) 


+ a,z  + Uq 


z=1^((nnT/M) 

(UOI) 


z MbT 


where  b + (2m/1 T). 


One  must  remember  that  the  terms  in  Eq.  401  do  not  cancel.  For  example , 
the  z - 1 in  the  numerator  of  the  last  term  does  not  cancel  a z - 1 in  the 
denominator  of  the  second  term,  since  each  bracketed  term  is  being  evaluated 
at  a different  value  of  z. 

For  comparison  purposes,  we  can  define  a reference  "analog"  system  by 
removing  the  sample  and  holds  from  Fig.  The  transfer  function  would 

then  be 


C G2[I  + HG2r  G,[l  + G2(I  4 HG2r'G,  f'R 


(1*02) 


For  this  scalar  illustrative  example  the  reference  transfer  function  becomes 


C_ 

R 


s^  + 2s  + 


(1*0-) 


The  Bode  plot  obtained  using  Eq.  h 03!  can  be  compared  against  those  obtained 
using  Eq.  *401 . Some  comparative  results  are  shown  in  Fig.  for  the  case 
where  M = 2.  Note  that  for  rates  as  low  as  10  samples/second  (T  - 0. 1 ) 
that  there  is  no  discernible  difference  in  the  magnitude  plot  — - at  least 
in  the  portion  of  the  first  fold  that  is  shown  in  Fig.  56.  However,  the 
difference  in  the  phase  plots  is  noticeable.  Of  course,  the  difference 
becomes  quite  large  when  the  frame  time  is  set  at  1 second  (T  - 1). 


1 ho 


Amplitude  (dB) 


Phase (deg) 


K.  GENERAL  RESULTS 


It  is  worthwhile  at  this  point  to  puii  together  in  a compact  form  the 
general  results  of  Section  V.  The  reader  may  be  a little  apprehensive 
over  the  prospect  of  practically  applying  the  detailed  techniques  outlined 
in  this  section  for  obtaining  the  frequency  response  coefficients  (An  + Bnj). 
It  may  appear  that  each  time  a new  system  configuration  is  encountered,  a 
detailed  algebraic  and  limiting  process  using  partial  fraction  expansion 
must  he  c.'irried  out  to  formulate  the  frequency  response  expression  for  a 
discretely  controlled  system.  Fortunately,  some  generalities  exist  that 
allow  us  to  bypass  much  of  this  detailed  derivation  process  and  replace  it 
with  a simple  substitution  procedure. 


The  first  observation  we  make  is  that  the  vector  algebraic  expression 
for  the  output.  (see  Figs.  39.  5?»  and  55)  is  configuration  depen- 

dent. but  can  t iike  on  three  forms: 

Open -Loop 

CT/N  (GM)T/NRT  (i*04) 


Closed-Loop  0 Ingi e-Rate 
CT/N  : (GM)T/N  gTrT 

Closed-Loop  Multi-Rate 

CT/“  . (OM)T/"gJ  ioj*!1 


(405i 


( 4061 


N/Mk  = Integer  ; k 0,  1,  2,  3,  ... 

T/Vk 

The  number  of  Gk.  terms  in  the  product  expression  equals  the  number  of 
different  sampling  rates  in  the  system  (excluding  the  "T"  sampling  rate). 

For  example,  if  a particular  system  has  samplers  operating  at  T,  T/iO,  T/PO, 
and  T/‘)0  (N  100) 


^Z'00  <gmit/'00  gIgT/'0  RT 


(*07) 


i hr* 


With  the  CT/N  expression  formulated  using  block  diagram  algebra,  the 
following  identities  allow  the  frequency  response  expression  (An+jBn)  to 
be  written  directly  from  inspection  of  each  term  in  the  CT/N  vector  equa- 


tion. 
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The  notation  used  in  these  identities  indicates  the  definition  of  "z" 
to  be  used  in  calculating  each  z-transform  expression  and  the  subsequent 
evaluation  performed  to  calculate  coefficients  for  the  fundamental  and  alias 
terms  in  the  output  waveform  (An+ JBn).  That  is, 

Take  the  T/N  z-transform  of  (GM)  and  evaluate  at 

z = l^u^T/N). 

Take  the  T z-transform  of  GA  and  evaluate  at 
z = l^bT. 

Take  the  T/'Mk  z-transform  of  each  Gk  term  and 
evaluate  at  z - Liii^T/M^) . 

Then,  for  the  (GM)  term,  directly  calculate  the  normal  z-transform  expres- 
sion via  the  partial  fraction  expansion/table  lookup  approach  of  Section  II 
and  then  replace  each  "T"  in  this  expression  with  "T/N."  Follow  the  same 
procedure  for  G^  and  Gk  with  "T"  replaced  with  "T/Mk"  in  the  resulting 
z-transform  expression  for  the  Gfc  terms.  Now  with  (GM),  G^,  and  Gk  in 
z-transform  form,  each  term  is  individually  evaluated  at  the  appropriate 
value  of  "z"  as  indicated  by  the  notation. 

Even  though  the  tables  of  Section  II  essentially  assume  z = es^,  these 
tables  can  be  used  for  any  sampling  period  T/N,  T/Mk,  etc.,  by  simply  replac- 
ing "T"  in  these  tables  with  T/N,  T/M^,  etc.  The  definition  of  "z"  is  then 
completely  arbitrary  (as  far  as  sampling  rate  is  concerned)  and  must  be 
explicitly  indicated  with  each  pulse  transfer  function,  i.e., 
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Thus  it  is  not  necessary  to  carry  out  the  detailed  algebraic  and 
Limiting  process  each  time  a new  system  configuration  is  encountered.  Once 
the  CT/N 

expression  has  been  determined  for  the  system,  simple  substitution 
produces  the  exact  frequency  response  expression  for  either  the  finite  N 
case  or  the  Limiting  case  of  N -»-«>.  Using  this  expression,  the  coeffi- 
cients (An  + jBn)  for  the  fundamental  and  its  aliases  present  in  the  output 
waveform  can  then  be  determined.  A summary  of  the  steps  required  is  as 
follows : 


l4r« 


| I 


• Obtain  the  CT/N 

expression  using  block  diagram/ 
signal  flow  algebra. 

| 

• Apply  the  identities  outlined  in  this  subsection 
to  each  term  in  the  C^’/N  expression. 

• Calculate  the  individual  z-transform  expressions 
using  the  appropriate  definition  of  "z." 

t • Evaluate  each  z-transform  expression  at  the  appro- 

priate values  of  "z." 

As  a final  remark,  it  should  be  noted  that  these  simplifications  still 
do  not  relieve  us  from  the  non-trivial  nature  of  obtaining  z-transform 
expressions  such  as  ( GMwJ^') T . This  type  of  expression  will  require  the 
use  of  either  switch  decomposition  or  the  phantom  sampler  as  illustrated 
in  preceding  subsections. 

L.  SECTION  SUMMARY 

The  "sampled  spectrum"  concept  of  sampled  data  control  theory  is  con- 
cerned with  determining  the  simple  sinusoid  that  fits  the  output  samples 
of  a single-rate  system  at  the  sampling  instants.  In  this  section,  the 
frequency  concept  has  been  extended  to  encompass  the  continuous  spectrum 
of  the  continuous  variables  in  a discretely  controlled  system.  Moreover, 
the  theory  considers  the  finite  N case  wherein  one  is  concerned  with  the 
group  of  N sinusoids  that  matches  the  data  not  only  at  the  sample  points 
but  at  N— 1 inter-sample  points  as  well.  This  is  an  important  aspect  since 
bench  validation  of  digital  hardware  is  often  specified  in  terms  of  an  end- 
to-end  "frequency  response."  Since  output  data  are  taken  at  a finite  number 
of  points,  it  will  be  important  to  compute  finite  N results;  the  coeffi- 
cients will  differ  significantly  from  the  continuous  (N  *♦«.>)  values. 

The  results  for  the  closed-loop  cases  have  been  configuration  depen-  j. 

dent;  however,  the  basic  technique  is  relatively  clear.  One  starts  at  the 
continuous  state  vector  and  writes  the  system  equations  back  to  the  first 


input  point.  The  next  fundamental  step  is  to  convert  that  input  into  an 
equation  which  contains  the  sinusoidal  input  as  the  basic  forcing  function. 
Then  one  may  always  invoke  the  basic  equations  which  apply  to  the  open-loop 
case.  The  equations  which  summarise  the  cases  considered  in  this  section 
are  given  in  Tables  8-11. 
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TABLE  8.  AN  OPEN-LOOP  CONFIGURATION 
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TABIE  9.  A SINGLE-RATE  CLOSED-LOOP  CONFIGURATION 
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SECTION  VI 


FLYING  QUALITIES  APPLICATIONS 


A.  INTRODUCTION 


As  noted  in  Section  I,  increasingly  stringent  operational  demands  on 
military  aircraft  have  forced  the  use  of  more  complex  stability  and  control 
augmentation  systems.  An  augmented  aircraft  may  have,  for  example,  short- 
period  characteristics  that  apparently  satisfy  the  MIL-F-8785B  specification 
and  yet,  because  of  the  additional  modes  introduced  by  the  flight  control 
system,  the  aircraft  actually  responds  as  if  it  had  short-period  character- 
istics which  do  not  meet  the  specification.  Examples  of  these  anomalies 
are  given  in  Ref.  U,  one  of  which  is  used  here. 


The  purpose  of  the  (specially  constructed)  examples  in  Ref.  1 was  to 
point  out  possible  loopholes  in  interpretation  of  the  short-period  frequency 
requirements  of  MIL-F-BtS'tB  (Para.  3.2.2. 1.1).  The  problem  stems  from  the 
fact  that,  although  the  specification  refers  to  "short  period  response  in 
angle  of  attack, " the  requirement  was  based  on  configurations  for  which  this 
response  was  defined  by  a natural  short-period  mode.  The  specification  does 
not  consider  additional  modes  which  might  be  introduced  by  the  FCS.  Since 
the  pilot  senses  the  total  response  (sum  of  natural  airframe  modes  plus  FCS 
modes),  specification  of  the  short-period  roots  alone  does  not  necessarily 
define  the  physical  situation  or  insure  acceptable  flying  qualities.  Now, 
the  trend  towards  the  use  of  digital  flight  control  systems  gives  the  flexi- 
bility which  tends  to  yield  additional  FCS  modes,  but  also  introduces  addi- 
tional effects  due  to  the  digital  implementation. 


One  method  proposed  (e.g.,  Refs.  1 and  1)  to  treat  such  situations  is 
to  develop  lower-order  models  of  higher-order  systems  that  can  be  compared 
against  the  MIL-F-878hB  requirements  as  they  are  presently  stated.  To  see 
how  this  method  works,  we  will  review  an  example  from  Ref.  1 which  deals 
with  the  effects  of  FCS  modes  on  the  MIL--F-878r--B  requirements.  This  con- 
trived example  considers  a tactical  fighter  which  requires  some  type  of 
flight  control  system  (FCS)  to  maintain  the  short-period  frequency  within 


acceptable  limits  at  low  speeds.  Pitch  rate  feedback  is  intentionally  used 
to  illustrate  the  possible  pitfalls  in  designing  a control  system  to  meet 
MIL-F-8785B  requirements  using  an  equivalent  low -order  system.  It  will 
become  obvious  that,  pitch  rate  feedback  alone  can  not  increase  the  effec- 
tive unaugmented  short-period  frequency. 

Tiie  subsequent  digital  implementation  of  this  example  from  Ref.  It  will 
not  attempt  to  cure  the  obviously  inadequate  analog  design,  but  will  attempt 
to  bring  into  sharp  focus  the  manner  in  which  the  equivalent,  low-order  model 
should  be  modified  in  order  to  account  for  the  additional  artifacts  intro- 
duced when  the  analog  controller  is  replaced  with  a digital  implementation, 
it.  is  our  intention  to  illustrate  how  the  techniques  developed  in  this  report 
can  be  used  to  include  the  effects  of  digital  Implementation  in  the  parameters 
of  a lower-order  model  of  airplane  dynamics  and  not  digitally  design  for  good 
flying  qualities.  The  same  general  flying  qualities  design  techniques  and 
philosophies  used  for  analog  controller  design  apply  equally  as  well  in  a 
valid  discrete  domain  synthesis  such  as  in  the  w'-domain  (Section  III'!. 

It.  is  readily  concluded  Unit,  to  use  the  equivalent  model  concept,  an 
exact.  Bode  plot  for  the  closed-loop  digital  system  which  includes  all  the 
effects  Of  the  digital  artifacts  (e.g.,  sampling,  data  holds,  and  computa- 
tional delay)  is  needed.  This  can  be  obtained  for  any  single-rate  or  multi- 
rnt.e  system  using  the  digital  frequency  response  methods  of  Section  V.  These 
methods  are  in  contrast,  to  ! lie  tradit  ional  concept  of  "sampled  spect  rum"  of 
sampled  data  theory  which  is  limited  to  determining  the  amplitude  and  phase 
of  the  single  sinusoid  that  fits  the  output  samples  of  a single-rate  system 
at  the  sampling  instants. 

The  example  used  in  this  section  wilL  be  closed  loop  with  a washout  net- 
work in  the  pitch  rate  feedback  loop.  The  Tustin  transform  is  used  to  obtain 
a digital  implementation  of  the  analog  washout  filter  to  illustrate  the  dis- 
advantage of  this  emulation  procedure  (i.o.,  the  Tustin  transform  does  not 
account  for  the  phase  lag  introduced  by  the  P/A  conversion).  The  effects  of 
throughput  delay  will  also  be  incorporated  into  the  analysis.  This  is  done 
in  order  to  indicate  a procedure  for  modeling  computational  delay  via  vector 
switch  decomposition  (Section  111  and  to  ilso  indicate  the  effects  of  the 
delay  on  the  closed-loop  Rode  plot. 


r ’ 


Even  though  the  washout  filter  is  emulated  l’or  illustration  purposes, 
it  should  be  noted  that  direct  digital  synthesis  in  the  w' -domain  (Sec- 
tion 111)  would  produce  lead  compensation  to  onset  the  throughput  delay 
and  nonminimum  phase  A/D,  D/a  effects.  This  is  assuming  that  u match  of 
the  already  established  inadequate  analog  design  is  wanted.  As  covered  in 
detail  in  Section  111,  direct  digital  design  in  the  w ' -domain  is  straight- 
forward following  t.he  some  procedures  as  an  analog  s-domain  design.  It 
should  be  noted  that  the  power  of  a w’ -domain  direct  digital  design  is  that 
all  effects  of  digit. al  implementation  (e.g.,  sampling,  data  holds,  and  com- 
putational delay)  is  directly  and  naturally  accounted  for  in  t.he  resulting 
w' -domain  control  law.  That  is,  the  resulting  design  in  the  w' -domain  is 
exact  in  the  same  sense  as  the  analog  control  law  is  in  t.he  continuous 
8 -domain  synthesis.  Then,  from  a practical  and  realistic  standpoint,  a 
complete  redesign  of  the  control  law  for  t.he  example  In  this  section  using 
different  feedback  schemes  (e.g.,  angle -of -attack  feedback!  should  be  done 
either  in  the  s-domain  for  the  analog  system  or  in  the  w1 -domain  for  the 
digital  system. 

B.  TACTICAL  FIGHTER  — ANALOG  CONTROLLER 

Current  trends  in  fighter  design  frequently  require  some  type  of  KCP 
to  maintain  t.he  short -period  frequency  within  acceptable  l imits  at  low 
speed.  The  airpLane  characteristics  used  in  this  example  are  typical  of 
modern  tighter  aircraft  in  the  power  approach  configuration.  The  attitude 
transfer  function  is  given  as: 


JL  -m._(o.i,Ko.(Q  , • (U.S) 

be  [0.0b,  0.1 1(0.7.  0.b|  (sl[0>7,  O.b]  v 

where 

[<,  to]  =£>  S-'  4 P^inS  » <!»'  ; ( 1 /T 1 s * |/T 

Come  type  of  compensat. Ion  is  desired  because  the  low  short -period  frequency 
violates  Mll..-K-tV(8*SH  and  produces  a sluggish  pitch  response.  As  we  will 
see,  feedback  of  lagged  pitch  rate  results  in  an  apparent,  increase  in  short- 
period  frequency  (and  damping).  Tills  signal  would  be  washed  out  to  prevent 
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the  pitch  damper  from  receiving  low-frequency  signals  which  might  saturate 
the  system.  The  example  FCS  is  shown  in  Fig.  ‘'7  along  with  a system  survey 
to  indicate  the  effect  of  gain  (Kg)  on  the  system.  Since  the  phugoid  fre- 
quency is  only  20  percent  of  the  short  period,  the  two-degree-of-freedom 
constant-speed  approximation  is  used  for  the  sake  of  clarity.  Closing  the 
loop  at  Kg  = -0.63  results  in  short-period  poles  at  1.36  rad/sec.  The  aug- 
mented attitude  dynamics  are  then  given  as: 

_L  - -1 ■ 6(0.6) (0.?) (2.0)  ( . 

6p  s(0.68)(0.2)[0.86,  1.3?] 

From  Fig.  r'7  it  appears  that  the  FCS  increases  the  short-period  undamped 
natural  frequency  and,  apparently,  specification  compliance  has  been  achieved. 
However,  an  equivalent  lower-order  system  which  gives  approximately  the  same 
frequency  response  is; 

(h20) 

This  effective  second-order  system  fitted  to  the  actual  system  has  a natural 
frequency  of  only  0 . ? rad/sec  and  is  critically  damped.  Plotting  the  appar- 
ent augmented  ,oSp  and  the  "effective  <oSp"  on  the  MIL-F-878rd3  requirement 
illustrates  the  danger  in  considering  only  the  short-period  pole  for  speci- 
fication compliance  when  a FCS  is  employed  (see  Fig.  ^8).  Note  that  even 
though  the  FCS  apparently  increased  mSp  from  0.r^  to  1.3?  rad/sec,  the  effec- 
tive response  indicates  that  pitch  rate  feedback  only  increased  d;unping. 

The  sluggish  characteristics  of  the  augmented  airframe  are  indicated  in  the 
pitch  rate  time  history  to  a step  input  of  &p  as  shown  in  Fig.  ?9-  Note 
that  the  time  responses  for  the  actual  and  effective  transfer  functions  are 
in  close  agreement. 

The  above  augmentation  scheme  is  obviously  unsatisfactory,  making  a 
sluggish  response  even  worse.  It  still  serves  our  purposes  to  consider  the 
digital  mechanization  of  this  system. 
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c.  TACTICAL  FIGHTER  — DIGITAL  CCIfTROLIZR 

Consider  a possible  model  for  one  implementation  of  a digital  controller 

for  the  tactical  fighter  shown  in  Fig.  60.  The  associated  timing  diagram 
is  shown  in  Fig.  6l . 

In  Fig.  6l  the  timing  events  are: 


0 5p  is  sampled  and  stored  in  buffer  register;  q is 
sampled  and  stored  for  computation  of  the  feedback 
portion  of  the  control  law. 

(2)  Computation  completed;  control  algorithm  output  to 
control  actuator  through  a zero-order  hold. 

In  Fig.  60  advantage  has  been  taken  of  the  fact  that  a sampling  operation 
which  is  delayed  with  respect  to  the  sampling  operation  of  a reference  sampler 
by  T0  seconds  can  be  modeled  by  switch  decomposition  (refer  to  Appendix  F) 
as  shown  in  Fig.  62.  From  Fig.  60  the  equations  which  pertain  can  be  written- 
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and  the  table  of  advanced  z-transforms  (Table  4,  Section  II)  shows 
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given  0 i Tq  < T. 


Figure  60.  Digital  Implementation  of  FCS 
for  Tactical  Fighter 


Figure  62.  Model  of  a Non-Synchronous 
Sampling  Operation 
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Equation  ■■  1 r'  simplifies  to 
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The  output  is  then  given  by 
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The  "continuous"  spectrum,  given  that  bp  is  a sine  wave  input,  can  be 
obtained  (using  the  tools  of  Section  V)  by  a limiting  procedure  on  q'/N. 
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The  spectral  coefficients  for  the  ' response,  0 = q(si.  can  then  be  com- 
puted using: 
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The  procedure  is  as  follows : 
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2)  Compute  HT  via  Tustin  transform. 
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U ) Evaluate 


GA  * 


An  + JBn  = y 


/GMc-i;To\ 

V“H  a-*,, G 


z= 14bT 


Step  1 


e~sTo  _ -(l-A)sT  -sT  AsT 

“e  =e  e ;Osiil 


(GMe"^To)T 


j -bt  ( ---  e~~bT\  H.6s  - .o61e^ 
\e  \ s f s(s-  + ,7s  + 


AsT  1T 


L:'  (“1.6s  - .9 6)e*1 

[s-(s-  + ,7s  + .25) 


— 1 (a12~  + S2g  + a^)z 

z':  [(-  “ 1 )(z2  + bpz  + bT,) 


(GMe~6T°)T  = -a1Z"  + a22  + a~* 
z(z 2 + b2z  + b^) 


(1*291 


Values  for  the  coefficients  of  Eq.  629  are  tabulated  as  a function  of  A in 

Table  12  [A  = (T  - T0)/Tj. 


TABLE  1; 


COEFFICIENTS  OF  (GMe_'sT°)T 


I a,  0.000000000  -0.006398718  -O.O3.9676, 8 -0.0^736233  -0.057^098  "-0 .0638689<i3 
*2  ■C-06W^U  -0.00069S.95  0.026,606,0  0.0S9W9*  0.0625, 6508 
«,|  0.0605,6,09 1 0.0*, 08192  0.05„S7,68  0.0,8630978  0.006226,07  O.OOCOOOOOO 

Notes  b2  =.  -,.971993928;  b5  = O.972388367 
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The  Tustin  transform  of  the  washout  network  is  a direct  substitution 
technique  wherein  one  substitutes 

2 z — 1 
8 T z + 1 

into  the  s-domain  transfer  function,  H(s).  The  result  for  T = 0.0^  sec 
(2^  samples/sec)  is; 


H(z)  = 


ap( z2  ~ 1) 

2 + b0z  + h} 


(**30) 


aQ  --  -0.011995 
b0  = -1.903275 
b,  = 0.901*798 


The  computation  of  could  be  carried  out  on  a point-by -point  basis 

since  the  problem  being  considered  is  a scalar.  However,  it  is  illuminating, 

T 

even  if  somewhat  disconcerting,  to  examine  the  numerical  expression  for 
wherein  polynomial  coefficients  are  rounded  to  three-place  accuracy: 

z-  - 3.875Z1*  + 5.632^  - 3.735Z"  + o.880z 


z5  _ 3.875s4  + 5.63z5  _ 3.6362'-  + O.88O2 

Thus,  even  when  using  a sampling  rate  as  low  as  25  samples/second,  it  is 
impossible  to  distinguish  between  the  characteristics  of  the  open-  and 
closed- loop  systems  with  only  three-place  accuracy  for  the  algorithm  coef- 
ficients. This  example  serves  to  warn  one  of  numerical  roundoff  problems 
that  may  be  encountered  in  digital  systems  even  at  moderate  sampling  rates, 
and  indicates  the  need  for  at  least  double  precision  computation  at  higher 
sampling  rates. 


Step  U 


The  frequency  response  is  computed  with  respect  to  0 rather  than  q. 

Since  0 = q/s,  the  computation  is  carried  through  using  ( 1 /T)  (GMe~sT°/s') 
rather  than  ( 1 /T)GMe~sT°.  The  result  is  shown  in  Fig.  63.  Notice  at  the 
relatively  slow  sample  rate  of  25  samples/sec  that  the  difference  between 
the  magnitude  plot  for  the  "continuous"  washout  network  and  the  discre- 
tized version  only  becomes  significant  at  frequencies  above  ^0  rad/sec. 

There  is.  however,  a substantial  difference  between  the  phase  angles  of 
the  continuous  washout  and  the  discretized  version  at  frequencies  above 
2 rad/sec.  This  result  is  a typical  one  when  the  Tustin  transform  approach 
is  used.  The  only  recourse  for  decreasing  the  difference  between  the  phase 
responses  when  the  Tustin  transform  approach  is  used  is  to  decrease  the 
sampling  period.  This  is  the  case  because  the  Tustin  transform  does  not 
account  for  the  phase  lag  introduced  by  the  A/D.  D/A  conversion  process 
(i.e.,  the  data  holds ^ . 

The  phase  for  the  discretely  controlled  system  is  shown  for  three  dif- 
ferent values:  Tc  = 0 (A  = 1 ) , T0  = 0.5T  (A  = 0.7'),  and  To  = T (A  = 0) . 
Clearly,  the  closed-loop  phase  characteristic  is  sensitive  to  the  through- 
put delay  of  T0  seconds.  The  accuracy  with  which  the  lower-order  system 
models  the  higher-order  system  would  be  improved  by  inclusion  of  a pure  time 
delay  term.  However,  specification  of  acceptable  delay  is  not  presently  part 
of  MIL-F-878^B.  Therefore,  further  adjustments  of  the  effective  short-period 
damping  and  natural  frequency  in  the  lower-order  model  are  required  if  the 
"digital"  phase  angle  artifacts  are  to  be  accounted  for  in  the  format  of 
MIL-F-8783B.  The  result  will  be  an  increase  in  the  effective  damping  and 
a decrease  in  the  effective  short-period  frequency.  Either  in  this  form 
or  considering  the  lag  shown  in  Fig.  63  the  result  is  further  degradation 
of  flying  qualities  due  to  the  digital  mechanization. 

As  discussed  in  the  introduction  to  this  section,  an  exact  digital 
implementation  would  require  an  independent  direct  digital  design  in  a 
valid  discrete  domain  such  as  the  w' -domain  (Section  III).  The  resulting 
digital  control  law  would  in  general  differ  from  the  simple  Tustin  trans- 
form of  the  analog  control  law,  because  the  w' -domain  synthesis  would 
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produce  additional  lead  compensation  to  offset  the  throughput  delay  and 
nonminimum  phase  A/D,  D/A  effects.  The  approximation  procedure  using  the 
Tustin  transform  cannot  account  for  these  digital  artifacts. 

0.  LOWER -ORDER  MODELS 

No  effort  has  been  made  herein  to  define  a rationale  for  modifying  the 
parameters  of  the  lower-order  model  in  order  to  account  for  the  digital 
artifacts.  The  reason  for  this  is  that  several  modeling  procedures  for  the 
generation  of  lower-order  models  already  exist  (see,  for  example,  Refs.  1 
and  4).  In  general,  the  frequency  domain  modeling  procedures  define  a fixed- 
form  model  and  then  use  machine -implemented  "fitting"  algorithms  to  numeri- 
cally determine  the  parameters  of  the  model.  The  input  for  the  optimization 
program  is  obtained  by  reading  magnitude  and  phase  angle  data  from  an  appro- 
priate Bode  plot.  Typically,  ten  frequencies  between  0.1  and  10  rad/sec  are 
selected.  This  range  is,  of  course,  well  within  the  first  fold  of  any  digi- 
tally controlled  system  using  sampling  rates  on  the  order  of  25  samples/sec. 

It  is  apparent  that  these  existing  modeling  procedures  do  not  have  to 
be  drastically  changed.  It  is  only  necessary  to  read  the  input  data  from 
Bode  plots  which  already  include  the  digital  artifacts.  These  Bode  plots 
can  be  generated  using  the  tools  of  Section  V. 

E.  CONTROL  ROUGHNESS 

No  assessment  of  control  roughness  (ripple)  has  been  made  for  the 
illustrative  example.  The  reason  for  this  is  that  the  short -period  fre- 
quency of  the  aircraft  model  is  considerably  less  than  the  sample  frequency 
of  157-1  rad/sec,  insuring  that  the  ripple  will  be  small.  However,  if  the 
methods  in  the  report  had  been  used  to  their  extent  to  minimize  sample  rate, 
it  may  have  been  necessary  to  consider  control  roughness.  This  can  be  done 
with  the  approach  of  Section  V. 

F.  SECTION  SUMMARY 

An  illustrative  example  was  used  to  demonstrate  the  application  of  the 
digital  tools  developed  in  this  report  towards  lower-order  system  modeling 
of  aircraft  and  control  system  dynamics.  This  example  highlights  several 
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salient  features  of  the  digital  analysis/synthesis  techniques  dealt  with 
in  this  report: 

• The  Bode  plots  generated  using  the  frequency  response 
tools  of  Section  V contain  all  the  information  required 
to  assess  the  effective  time  delays  and  control  rough- 
ness resulting  from  digital  implementation.  The  phase 
plot  quantifies  the  effective  time  delay  and  the  magni- 
tude plot  the  control  roughness. 

• For  assessing  MIL-F-8785B  requirements,  the  effects  of 
the  digital  artifacts  on  the  equivalent  low-order  model 
can  be  accounted  for  exactly  using  the  digital  frequency 
response  techniques  of  Section  V. 

• Although  a single -rate  example  was  used,  the  frequency 
response  methods  of  Section  V are  also  applicable  to. 
multi-rate  systems.  This  is  in  contrast  to  the  tradi- 
tional concept  of  "sampled  spectrum”  of  sampled  data 
theory  which  is  limited  to  single-rate  systems. 

• A direct  digital  design  is  needed  if  the  phase  lag. intro- 
duced by  the  throughput  delay  and  a/D,  D/A  conversion 
process  is  to  be  offset  through  proper  control  law  design. 
The  w' -domain  synthesis  methods  of  Section  III  accomplish 
this,  whereas  the  emulation  procedure  using  the  Tustin 
transform  neglects  these  digital  phase  angle  artifacts. 

• The  procedure  for  modeling  computational  delay  is  easily 
handled  using  vector  switch  decomposition  coupled  with 
the  advanced  z-transform  methods  of  Section  II. 
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SECTION  VII 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A.  SUMMARY  AND  CONCLUSIONS 

The  criteria  of  MIL-F-8785  presume  an  analog  augmentation  system. 

There  is  now  an  additional  concern  for  response  characteristics  which  are 
unique  to  digital  augmentation  systems.  The  purpose  of  this  study  was  to 
identify  those  characteristics  which  are  unique  to  the  digitally  controlled 
system,  and  to  review  the  quantitative  tools  available  which  permit  an 
assessment  of  the  effects  of  these  characteristics  relative  to  the  MIL-F- 
8785  requirements. 

Two  important  characteristics  are  introduced  by  digital  control  laws: 

• The  effective  delay  introduced  by  the  A/D  and  D/A 
process  and  the  delay  introduced  by  the  digital 
algorithms  and  computational  frame  time. 

• The  "control  roughness"  or  inter-sample  ripple 
introduced  when  the  digital  computer  is  coupled 
to  the  control  actuators  using  data  holds. 

The  first  characteristic  is  of  great  concern  to  the  flying  qualities 
community  since  even  relatively  small  delays  are  potentially  important  in 
certain  closed-loop  piloting  tasks  involving  motion  cues  and/or  pilot - 
induced  oscillations.  The  second  characteristic  is  of  concern  since  air- 
craft response  resulting  from  control  roughness  is  in  effect  an  additional 
disturbance  source. 

Measures  of  the  effective  time  delay  introduced  by  the  A/D  computation, 

D/A  conversion  process  can  be  computed  using  a variety  of  analytical  tech- 
niques that  apply  at  the  sampling  instants.  For  example,  w’-,  w-,  or  z-domain 
analyses  and  discrete  frequency  response  techniques  can  be  applied  effectively 
to  obtain  both  a quantitative  and  qualitative  evaluation  of  the  delay  intro- 
duced . 

The  second  characteristic  is  more  difficult  to  assess  since  it  involves 
the  response  of  the  continuous  system  during  the  inter -sample  interval. 
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Traditional  analysis  tools  such  as  the  "sampled  spectrum"  (discrete  fre- 
quency response)  are  of  little  value  for  this  assessment  because  they  gener- 
ally ignore  continuous  system  response  between  samples. 

This  leads  to  the  other  main  objective  of  this  study  — to  encourage 
the  practicing  engineer  to  gain  working  familiarity  with  three  analytical 
techniques  well  suited  for  the  analysis  of  digitally  controlled  systems: 

• Analysis  (and  synthesis)  in  the  w'-domain. 

• Multi -rate  transform  domain  analysis. 

• fvaluat ion  of  continuous  system  frequency  response 
to  discrete  excitation. 

The  w'-domain  (Section  III)  is  related  to  the  well-known  w-domain  by  a 
scaler  transformation  and  to  the  z-domain  by  a bilinear  algebraic  transfor- 
mation. We  believe  that  those  engineers  skilled  in  frequency  domain  design 
procedures  will  immediately  feel  at  home  with  analysis  in  the  w'-domain. 
since  all  analog  control  system  design  technology  transfers  completely  for 
digital  control  system  design.  In  the  w'-domain,  the  "analog"  control 
system  designer  is  transformed  instantly  into  a "digital"  control  system 
designer. 

The  second  item  deals  with  a multi-rate  transform  domain  approach  that 
we  developed  (Section  IV)  into  an  effective  tool  for  analyzing  not  only 
multi-rate  sampled  closed-loop  systems  (Section  IV),  but  also  the  transient 
inter-sample  response  of  discretely  excited  systems  (Section  V).  That  is. 
it  yields  recursion  equations  describing  the  inter-sample  performance  to 
;uiy  degree  of  fineness  desired  without  increasing  computer  storage  require- 
ments . 

The  latter  use  is  a departure  point  (Section  V)  to  extend  frequency 
response  concepts  to  include  the  continuous  frequency  response  of  a dis- 
cretely controlled  system.  The  discrete  frequency  response  concept  has 
not  been  particularly  productive  for  the  analysis  of  discretely  controlled 
systems,  since  it  is  limited  to  determining  the  amplitude  and  phase  of  the 
single  sinusoid  that  fits  the  output.  samples  of  a single-rate  system  at  the 
sampling  instants.  Development  proceeded  by  first,  removing  this  restriction 
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for  open-loop  systems  and  then  extending  the  results  to  single-rate  closed- 
loop  systems.  Finally,  the  solution  for  the  multi-rate  closed-loop  case 
was  developed  using  the  results  of  Section  IV. 

B.  RECOMMENDATIONS 

There  are  two  specific  recommendations.  First,  there  is  a need  for 
portable,  general  purpose  software  to  support  the  frequency  response  analy- 
sis of  Section  V.  The  basic  requirements  for  one  implementation  are: 

1)  Compute  the  matrix  of  open-loop  transfer  functions, 
in  the  s-domain,  given  equations  in  either  a state 
variable  or  degree-of-freedom  format. 

2)  Partial  fraction  expand  each  eleinent(s)  of  this 
matrix,  compute  the  advanced  z-transform  of  each 
term,  and  then  combine  the  individual  terms  to 
obtain  a polynomial  in  z.  The  result  is  a matrix 
of  transfer  functions  in  z for  the  discretized 
system. 

5)  Given  the  matrix  of  transfer  functions  in  z,  com- 
pute the  w’ -matrix  of  transfer  functions. 

An  alternate  implementation  might  replace  the  transform  domain  computations 
in  Items  1 and  2 with  state  transition  computations  which  produce  equiva- 
lent results.  These  requirements  are  easily  satisfied.  It  is  primarily  a 
matter  of  integrating  existing  codes  in  a manner  which  will  permit  easy 
interaction  between  the  user  and  program. 

The  second  recommendation  deals  with  the  application  of  switch  decom- 
position to  model  data  skewness  and  throughput  delay  effects  in  digital  com- 
putations. This  approach  to  modeling  multi-rate,  multiple  order,  throughput 
delay  and  data  skewness  effects  has  vast  potential  for  the  effective  analy- 
sis of  complex  systems  exhibiting  these  phenomena.  Only  a hint  of  the  power 
switch  decomposition  teamed  with  the  continuous  frequency  response  concept 
is  given  by  the  simple  illustrative  example  of  Section  VI.  Specifically, 
the  switch  decomposition/frequency  response  techniques  in  combination  can 
(and  should)  be  developed  into  an  analytical  approach  for  pinpointing  error 
sources  that  arise  because  of  digital  effects  in  the  control  and  simulation 
of  continuous  systems. 
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APPENDIX  a 

PROPERTIES  OF  THE  LAPLACE  TRANSFORM 


INTRODUCTION 

Some  of  the  more  Important  properties  of  the  Laplace  transform  that 
„ere  not  discussed  in  Section  II  will  be  reviewed  here.  Of  particu  ar 
interest  are  the  properties  of  ..translation  in  time"  and  "periodic  func 

tions. " 

LAPLACE  TRANSFORM  OF  THE  DERIVATIVE 
OF  A TIME  FUNCTION 

The  definition  of  the  Laplace  transform  of  a time  function  f(t)  is 

<£  [f(t)]  £ F(s)  = f f(t)e'St  dt 

J O 


(A-1) 


Given  F(s).  the  transform 
Let 


of  df(t)/dt  is  found  using  integration  by  parts. 


u = e 


-st 


df  a + 
dv  = gj  dt 


(A-2) 


so  that 


J u dv  = mv  - / v du 


f(t)e  St  + s 
o 


f°°  f(t)e  St  dt  (A-5) 

J n 


Therefore 


Xj||  = sF(s)  - f(0) 


(A-U) 


giving  the  pair 


df 

dt 


sF(s)  - f(0) 


(A-?) 


INITIAL  AND  FINAL  VALUE  THEOREM 


The  pair  given  in  Eq.  A-5  can  be  now  used  to 


prove  the  equalities 


Lim 


t^w  = f(oj  - 8F(b) 


(a-6) 


Lim 


f(t)  - f(oo)  = s_^o  sF(s)  , poles  of  sF(s) 

in  LHP 


Since 


(A-7) 


/: 


df  -st 

dt  e dt  = sF(s)  - f (o) 


(A-8) 


we  may  write 


sF(s)  - f (o)  = Lim  C 

t^~  Jo  du 


e— su  du 


(A-9^ 


so  that  letting  s -»-  » immediately  gives 


Lim 

S -fc-oo 


sF(s)  - f(o)  = 0 


which  proves  Eq.  A-6. 

Next,  let  s — o in  Eq.  A-9,  giving 


Llm  sF(s)  — f(0)  = tLim  rt 


s -»-o 


f - 

du 


du 


Lim 


(A- 10) 


i-J  -*■  Ul  r ..  . \ 

t U(t)  - f(0)]  (A-n) 


proving  Eq.  A-7. 

The  exigence  of  the  limit  as  . - o require,  the  restriction  on  the 
I’O  M o.  »F(s).  The  limiting  procedure  of  Eq.  A-7  may  give  a "nonsense" 
answer  if  this  restriction  is  not  observed.  For  example, 


1 T1 


so  that 


sin  bt  <$=> 


b 


s^  + bc 


f(oo) 


Lim 
s -»-o 


sF(s) 


Lim  sb 
s -^-o  s2  + b^* 


(A-12) 


(A-13) 


Since  the  poles  of  sF(s)  are  not  in  the  left-half  plane  (they  are  on  the 
imaginary  axis),  the  restriction  is  not  satisfied,  and  the  result  of  Eq.  A-13 
is  meaningless.  This  is  as  it  should  be,  since  no  one  knows  the  value  that 
the  sine  wave,  although  bounded  by  ± 1 , takes  on  when  t = ». 

HIGHER  DERIVATIVES 

The  results  of  the  second  article  in  this  appendix,  page  171  are  suffi- 
cient for  finding  the  transform  of  higher  derivatives.  Consider 


X 


where  g(t)  is  set  equal  to  df/dt.  Therefore. 

<j=o  sG(s)  - g(0) 


(A-lM 


(A-?1') 


which  gives,  upon  substitution  for  g, 


or 


dg 

dt. 


«=>  s[sF(s)  - f(0)  ] - 
(It/ 


1=0 


(A-16) 


s-  F^ s) 


sf(0) 


f ' (0) 


(A-17) 


I '3 


\ 


For  example, 


X - f"(0)  = s5F(s)  - S2f(0)  - sf’(O)  - f"(0)  (A- 

dt5  (dt^ 

and 

x = sX  ~ - f”’(0)  = s!lF(s1  - s^f(O)  - s2f’(0) 

dt'1  dt5 

- sf "(0)  - f"’(0)  (A' 

As  an  application,  consider  the  time  domain  differential  equation 

x’  + fix  + 1 lx  + 6x  = 0 l A 

where 

x(0)  = 2 . x(O)  = -6  . x(0)  = 20 

This  transforms  into  an  algebraic  function  of  s 
[s^X(s)  — s2x(0)  — sx(0)  — x^O1)] 

+ 6[sr’x(s>  - sx(0)  - x(0)]  + 1l[sX(s)  — x(0)  ] + 6x(s)  = 0 (A 

or 

[.3  + 6s2  + 11s  + (']X(s)  = s2x(0)  + s[x(o)  + (oc(0)  ] 

+ [x(0)  + (ix(0)  + 11x(0)]  (A 


Substituting  the  initial  conditions  and  dividing  through  by  the  character 
istic  equation  gives 


Xts1! 


2s2  + ns  ■*  h 

(s  4 1 ) ( S 4 ?)(s  4 J>) 


2s2  4 hs  4 (> 

S ^ 4 6s2  + 1 1 S 4 6 


(A 


Partial  fraction  expansion  gives 


(A-2h ) 


so  that  the  time  solution,  via  table  lookup,  is 


TRANSLATION  IN  TIME 


Given  f(t)  u(t.),  the  transform  of  the  translated  function  can  be  easily 


found 


x + a 


Thus 


Given  that  a -i  0,  the  application  of  Eq.  A-?9  is  neat  and  straightfor- 
ward for  finding  the  transform  of  a delayed  time  function.  However,  given 
a 2 0.  the  transform  of  an  advanced  time  function  must  be  treated  with 


greater  care.  Using  the  same  procedures  as  above  one  finds 


Let  t + a = x,  so  that  dx  = dt  and  t = x-a.  Substituting  into  Eq.  A-30 
gives : 

£ [ f(t + a)u(t) ] = f f(x)e  S^X  dx 

•'a 

= eas  j f(x)  u(x— a)e-SX  dx 


* eas  £[f(t)u(t-a)] 


(A-31) 


Comparing  Eq.  A-31  with  Eq.  A-?9>  we  see  that  the  advantage  of  using  the 
translation  theorem  has  disappeared.  Namely,  in  Eq.  A-29  we  need  not  recom- 
pute the  transform  of  the  delayed  time  function.  It  is  sufficient  to  know 
the  transform  of  the  unshifted  function  since  a simple  multiplication  by 
e aS  gives  the  desired  result.  However,  in  Eq.  A-31,  a premultiplication 
by  eas  does  not  remove  the  necessity  for  recomputing  the  transform.  Indeed, 
we  may  as  well  use  Eq.  A-30  directly  when  computing  the  transforms  of  time 
advanced  functions,  since  no  advantage  accrues  from  the  use  of  Eq.  A-3l • 

PERIODIC  FUNCTIONS 

By  definition,  a function  satisfying 


i 


f(t)  = f(t  + T) 


(A -32) 


is  periodic  with  period  T.  This  definition  is  adequate  even  if  f(t)  =0 
for  t < 0.  Begin  by  defining  a function  equal  to  f(t)  in  its  first  period 
and  zero  elsewhere  (see  Eig.  A-l). 


f^t) 


f(t)  0 < t < T 


0 Elsewhere 


(A-33) 


Figure  A-1.  A Periodic  Function 


From  Fig.  A-1  we  see  that  f(t).  Fig.  A- la,  can  be  described  in  terms  of 
f^t).  Fig.  A-lb, 


f(t)  = f,(t)  u(t)  + f ^t-T)  u(t-T)  + ... 


Taking  the  transform  of  Eq.  A-3!i  gives 


F(s) 


F,(sl  + F^sle  sT  + F^sle^"''  + ... 

,,  / \ r 1 sT  — s2T 

F 1 ( s ) [ 1 + e + e +...] 


FIs') 


1 


Fi(s) 

-sT 

— e 


(A 

(A 

(A 

(A 


In  using  Eq.  A-37,  it  is  important  to  develop  the  transform  of  F)(sl 
over  its  entire  first  period,  taking  careful  note  of  the  "0  elsewhere" 
restriction  in  Eq.  A-33.  We  defer  making  an  application  of  Fq.  A-37  unti 
the  Laplace  of  the  impulse  function  has  been  introduced. 


* 


' 


We  observe  the  convention  that  t>( t ) = Lim  s(t  — a1!  and  is  therefore  included 

a-»-o 

within  the  limits  of  integration.  Specifically, 


S(t) 


fA-44> 


An  exajnple  is  given  in  conjunction  with  the  previous  subsection.  We  wish 
to  find  the  transform  of  the  "periodic"  string  of  5 functions  shown  in 
Fig.  A -2.  Since  the  Laplace  transform  of  a periodic  function  is 


F(s) 


Ft  (s) 


1 - e 


-sT 


(A-45) 


it  is  only  necessary  to  find  the  transform  over  one  period.  With  the  nota- 
tion 

f(t)  = 5+(t)  + 5+(t  — T)  + 5+(t-2T)  + •••  (A-46) 

to  indicate  a limiting  process  towards  the  left  in  time. 

Fi(s)  = 1 


j 


* 


and  we  have 

f(s)  = (A-I47) 

1 — e 

as  the  transform  of  a "comb"  of  5 functions. 


Figure  A-2.  The  Dirac  Delta  Comb  Function 
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r 


H 


THE  TRANSFORM  OF  TIME  PRODUCTS 


The  inversion  integral  for  Laplace  transforms  is 


f c+ j» 

f(t)  = X-1  [F(s)]  = 2^  J F(s)e8t  ds 

c — j ® 


(A-U8) 


To  find  the  Laplace  transform  of  the  product  of  two  time  functions  proceed 
as  follows: 

Xtf^tjfaft)]  = J f2(t)e"fit  j^ij  J [Fl(p)ept  dpjjdt  (a-49) 

Interchange  the  order  of  integration  and  adjust  the  contour 


£Cf1f2l  = 577  f Fi(p)  dp  f f2(t)e  (S-P)t  dt  (A-50) 

J C-j»  •'n 


Therefore 


, rc+j~ 

& [f (t)f2(t)]  = J _ F1(p)F2(s-p)  dp 


(A-51) 


For  example,  let  f-|(t)  = e_at  and  let  f2(t)  be  the  Dirac  S comb  function 
of  the  previous  section,  so  that  the  product  of  f-jf2  has  the  general  appear- 
ance shown  in  Fig.  A-3.  The  periodic  property  gives 


and  since 


Eq.  A-51  gives 


F2(s)  = 


Fi(s)  = 


1 — e 


s + a 


-sT 


(A-52) 


(A-53) 


rc+joo  F (p^  1 /• 

Xt f lf2]  = 2jtj  Jc_Joo  1 _ e-T(s-p)  dp  = 2^j  Jc 


1 c+ > dp 

^c-j~  (p  + a)[l  -e~Jr<'S-?h 

(A-5*0 


l8o 


I 


Figure  a-3-  Product  of  e 


with  Comb  Function 


Closing  the  contour  to  the  left  gives  (see  next  subsection) 


[ f 1^*2^ 


(A-?5) 


and  we  see  that  the  use  of  the  Laplace  transform  of  a product,  taken  in  con- 
junction with  the  Laplace  transform  of  a particular  periodic  function,  is 
a sufficient  tool  for  finding  the  transform  of  an  impulse  modulated  waveform. 
The  next  subsection  reviews  some  of  the  fundamentals  for  evaluating  Eq.  A-fl 
when  one  of  the  time  functions  in  the  product  is  the  Dirac  delta  comb  function. 


REVIEW  COMMENTS  ON  CONTOUR  INTEGRATION 


The  properties  of  sampled  systems  have  been  developed  directly  in  terras 
of  the  delay  operator  z 1 in  Section  II.  It  has  now  been  shown  that  the 
Laplace  transform  properties  associated  with  periodic  functions  and  with 
the  product  of  two  time  functions  can  be  used  to  develop  an  alternative 
frequency  domain  description  (by  letting  one  of  the  product  functions  be 
the  Dirac  delta  comb).  Representation  of  the  product  of  a continuous  time 
function  with  a Dirac  delta  comb  is  cast  as  an  equivalent  "impulse  modula- 


tion" model  in  Fig.  A-U.  It  has  already  been  shown  that 


r (t) 

rT(  t) 

R(s) 

T 

RT(s) 

Figure  A-1*.  Representation  of  Impulse  Modulation 

If  the  poles  of  R(s)  are  finite,  one  can  say  that  RT(s)  is  the  summation 
of  the  residues  of  R(p)  by  closing  the  contour  of  integration  to  the  left. 

On  the  other  hand,  if  the  contour  is  closed  to  the  right,  enclosing  the 
poles  of  the  Dirac  delta  comb,  one  obtains 

KT<«>  « E (A-57) 

n=  -a*  ' 7 


Equations  A -‘'6  and  A-57  are  useful  in  investigating  the  properties  of 
multi-rate  systems,  and  it  is  thus  worthwhile  to  gain  some  insight  into 
the  "mechanics"  of  evaluating  Eq.  A-56. 

First,  some  review  points.  If  a pole  is  enclosed  with  a counterclock- 
wise contour  giving  a residue  k-|,  then  a clockwise  enclosure  gives  the 
residue  -k^ . That  is,  a counterclockwise  contour  gives  a value  equal  to 
2nj  T.  Res  » whereas  a clockwise  contour  gives  a value  of  -2nj  T,  Res  . Also, 
if  the  integrand  is  represented  as  the  ratio  of  two  polynomials,  N(p)/D(p), 
then  the  residue  for  any  simple  factor  (p-p0)  is  given  by 


Res 


P=Po 


N(p0) 


(A-58) 


where  the  prime  denotes  differentiation  with  respect  to  p.  In  the  interest 
of  brevity,  the  discussion  will  be  restricted  to  systems  containing  simple 
poles  (no  multiple  poles).  Frequent  use  will  be  made  of  the  identity 


£ 


n=-« 


1 

X — ?nnj 


1 


18? 


(A-‘-9) 


Given  the  integral  vanishes  on  the  .nfinite  contour  (see  Fig.  A-'?) 
enclosing  the  poles  of  R(p),  we  may  say  that  RT(s)  is  equal  to  the  sum  of 
the  residues  of  R(p).  If  the  path  is  closed  to  the  right,  the  infinite 


\ C2 

X 

\ 

\ 


Figure  A-5.  The  Complex  p Plane 

set  of  poles  at  p = s-2jtnj/T,  n • 0,  ±1,  ±2,  of  1 - are 

enclosed.  First,  complete  the  contour  to  the  right,  enclosing  the  singu 
larities  at  p ■ s-2*nj/T.  Let 


Ni(p) 


- Di(p)[1 


where  it  is  assumed  that  R(p)  is  s proper  rational  polynomial  in  p. 


8 


1 _ e-T(s-p)  = 1 _ e-T(2«nj/T) 

p = s-(?nnj/T) 


1 - cos  ?nn  + j sin  2wn  = 0 (A-62) 


so  that 


ResL  . -M.I 


|p  = s-(2jtnj/T) 


n = 0,  ±1, 


Since  the  contour  was  closed  to  the  right,  -2nJ  £Res  gives 


"T<->  * T t *('-*-¥■) 

n^-«  \ / 


which  waa  to  be  shown.  We  therefore  have  the  identities 


Bl(s)  - 53  X 


c-j«  1 - e 


-T(s-p) 


T £ 

n=-«  ' ' 


For  example,  let  R(s)  = l/(s  + a)  and  close  the  contour  Ci 


RT(«) 


_i_  r+J* dp_ 

•'c-jco  (p  + a)ri-e“ 


■J"  " P * -a 


1 — e 


-T(a-P) 


P*  -a 


1 - e 


-T(s+a) 


(A-65) 


(A-64) 


(A-6?) 


Note  that  the  substitution  z = eSi  then  gives 


Equation  A-66  was  arrived  at  using  conventional  residue  theory.  The 
N(p)/D'(p)  approach  also  works,  since  N(s)  = 1 and  (again  closing  the  con 
tour  on  Ci)  gives 


] + (p  + a)  (-T)e 


This  result  can  also  be  obtained  with  the  C2  contour 


] + ( — T ) (p  + a)e 


p = s-(2jtnj/T) 


Therefore,  taking  due  account  of  the  clockwise  contour,  we  obtain 


where  Eq.  A-59  has  been  used  to  place  the  infinite  summation  in  closed  form 


THE  PULSE  TRANSFER  FUNCTION 


INTRODUCTION 


The  identities 


Figure  B-1 . An  Open-Loop  Two-Rate  System 


sxnaxi -rah  pulse  transhr  function 


In  Fig.  B-1,  let  Ti  ■ Tg  = T and  write 


4 


Now,  let  us  look  at  R^s  - (2jtnj/T)],  via  Eq.  B-1, 


'(•-m  - £ f 


c-joo  1 - e 


R(p)  dp 
r[s-p-(2jtnj 


(B-3) 


- -1-  r 

2nj  J 


f RU>)  dP 

c->  1 - e-^8-^  e2*^ 


(B-U) 


But,  e0,Tnj  = (cos  2jm  - j sin  2jm)  = 1 for  n,  an  integer.  Therefore, 


nT(s-Sfi)  , „t(s) 


(B-?) 


Substituting  Eq.  B-f  into  the  right-hand  side  of  Eq.  B-2  gives 


i £_  s(.  - ^)sT(.  - ^fi)  - [-fEo(.-^)pU)  - st(s)bt<») 


(B-6) 


Since  R^(s)  factors  out  of  the  summation,  we  see  that  the  remaining  term 
is  the  transform  of  the  g(t)  "impulse"  response.  Therefore,  it  has  been 
proven  that 


CT(s)  = [G(s)RT(s)]T  = GT(s)RT(s) 


(B-7) 


At  this  point,  it  is  convenient  to  discuss  (using  Eq.  B-1)  two  other 
points  occurring  in  Sections  IV  and  V which  the  reader  may  find  trouble- 
some. The  first  concerns  the  use  of  the  notation  RT(s)  rather  than  R(z). 

SUPERSCRIPT  ROTATION 

The  superscript  notation  has  been  adopted  as  a matter  of  convenience 
since  at  times  it  will  be  expedient  to  redefine  the  delay  operator  z 1 . 
For  example,  it  will  often  be  to  our  advantage  to  define 


-sT/N  -1 
e ' - z 


(B-8) 


so  that 


e-sT  = e-(sT/N)N 


- te-^T/")N  , 


(B-9) 


Equation  B-l  offers  a rationale  for  the  superscript  notation,  given  that 
we  desire  to  keep  our  options  open  on  the  definition  of  z.  To  demonstrate, 
we  seek  the  pulse  transfer  function  of  l/(s  + a)  using  Eq.  B-1  (closing  the 
contour  to  the  left): 


c-j°°  (p  + a) [ 1 -e  T^s"p  ] 


1 - e 


*r(  s-p) 


(B-10) 


Note  that  Eq.  B-10  contains  no  "z's";  the  only  "sampling  parameter"  is  T 
itself.  Only  with  the  subjective  decision  to  define  z = es-’  can  Eq.  B-10 


be  written  as: 


G(z  ■ = 


If  z had  been  defined  according  t.o 


z = e 


z = esip 


(B-l 1 ) 


(B- 12) 


— T s+ai 


N „-aT 
z — e 


(B-13) 


Next,  suppose  one  seeks  the  pulse  transfer  function  of  1 '(s  * a'  with 
respect  to  aT  N frametime.  Then. 


I 


1 lT/"  > 


-aT/N 


z = e8T/N 


(B-15) 


In  the  text  the  prevailing  definition  of  z will  be  omitted  when  it  is  felt 
that  the  meaning  is  clear. 

The  implications  of  Eqs.  B-11,  B-13,  and  B-15  are  as  follows.  Proceed- 
ing according  to  the  partial  fraction/table  lookup  approach  of  Section  II 
essentially  assumes 


z = e' 


(B-l6) 


since  the  table  was  written  with  that  definition  in  mind.  However,  use  of 
Eq.  B-l  gives  answers  in  terms  of  esT  or  esT/^  (whichever  is  appropriate), 
so  that  one  can  assign  z at  his  convenience.  If  we  wish  to  continue  to  use 
the  tables  of  Section  II,  and  yet  keep  the  definition  of  z an  option,  then 
one  proceeds  as  follows.  Suppose  that  it  is  convenient  to  define  z = esT/N. 
Then  is  computed  directly  from  the  tables  using  partial  fraction  expan- 

sion. Simply  replace  (mentally)  T in  the  tables  with  T/N.  On  the  other 
nand.  G^  is  computed,  by  table  lookup/partial  fraction  expansion  in  the 
usual  manner,  retaining  the  T in  the  transform  table.  When  the  overall 
pulse  transfer  function  has  beer,  computed,  one  simply  changes  z to  zN. 

THE  MULTI -RATE  CASE 

A logical  point  has  now  been  reached  to  derive  a particular'  relationship, 
pertaining  to  multi-rate  systems,  that  will  be  of  value  in  Sections  II,  IV, 
and  V.  Consider  the  two-rate  system  of  Fig.  B-1.  Write 


C(s)  = G(s)  R ’(s) 


(B-17) 


I 


1 


so  that 


CT2(s)  = [G(s)  RTi(s)]T’ 


(B-18) 


Using  Fq.  B-l,  one  may  derive  a relationship  between  Ti  and  Tp  which, 
when  satisfied,  permits  the  simplification 

To 

CT2(s)  = [G(s)  RTi(s)]  ‘ = GT2(s)RTi(s)  (B- 19) 

Applying  Eq.  B-l  to  Eq.  B-l8  gives 


, , ,_Ti,  ,-,T2  1 fC+J“  G(p)Ri1(p) 

tG(s'R  '<-'i  ■ ™ Jc_,„  , . >r£pi 


* £ £ o(.-2gi)BT'(s-2fi) 

n=  -a>  ' o-  ' ' c ' 


(B-20) 


Now,  take  a close  look,  via  Eq.  B-l,  at  R ^(s  — 2Knj/Tp).  By  definition, 


(8-z«i i\Ti  = -L  r 

\ To  )\  2nj  Jc_ 


R(p)  dp 


, , -Ti[s-(2itnj/T2)-p] 

300  i — e 


L 


a2«njT 1 /Tp  _ 


Let  T1/T2  be  an  integer, 


c+^°°  R(p)  dp 

_joo  1 _ g-^lls-p)  J-eitnjT,^ 


if  T^/Tp  is  an  integer 


( B— 2 1 ) 


(B-22) 


(B-23) 


RTl|s_2gij  = rT1(s) 


(B-21*) 
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Equation  B-20  becomes 


J2 


J_ 

T2 


£G(s-affli)p.(8)  , gt2(s)rti(s) 


(B-25) 


It  has  thus  been  proven  that 


C [G(s)fi  '(s)]  = G^(s)R1(s)  , T-|/Tp  = integer  (B-26) 

giving  vhe  ruie  that  the  outer  sampling  operator  moves  through  the  inner 
one  when  the  ratio  of  the  inner  to  the  outer  is  an  integer. 

For  example,  let 


R(s)  = 


s + a 


G(s)  = 


s + b 


(B-27) 


When  Ti/Tj)  is  an  integer,  we  expect 


.T2  _ 


(j  _ e^2(s+bj][l  _ e^Ti(s+a)] 


(B-28) 


as  can  be  verified  using  Eq.  B-1. 


-T2 


j_  rc+ 
2*J  I* 


C+Joo 


dp 


c-j«  (p  + b)[  1 - e-111"1  ^p+a^][i  _e-T2(s-P)] 


(B-29) 


Therefore, 


_N_ 

O' 


[,_e^>(P^^[,_e-jr2(B-p)]+Tl(p  + b)e-jr,(iHa)[1_e^r?(s-p^  + (p  + b)[i_eJr,(p+a)](_jr2)eJr2(s.p) 


(B-30) 
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Closing  the  contour  to  the  right,  enclosing  the  singularities  at 
p - s - (2jtnj/T2),  gives 


Ts,-*2(2«nj/rS>[1  -e'Jr'[st*-(2”nJ/T2))][..b-(2»n3A2)] 


(b-31) 


n a 0,  i- 1,  ±2,  . . . 


s^T2(2nnj/T2)  = 1 ? ^ (2*nj/T2)  = y 


so  that  we  are  left  with 


(B-32) 


1 f ] 

(■,  _ e-Ti(s+a)j  T2(s  +b)  - 2«nj 


(B-33) 


Using  Eq.  A-*9  gives 


,T2 


[ 1 - e-T1(B+a>][l  - e-^2<8+b)] 


(B-34) 


If  Ti/T2  is  not  an  integer,  we  are  stopped  at 


T2  _ y.  ' 

rt.  [1  _e-ri(»*a)s2«njT,/T2]  [T?(atb)  _ 8im))] 


(B-55) 


If  T2 /T ■(  is  an  integer,  an  answer  can  be  obtained  by  closing  the  contour 
to  the  left.  Doing  this  encloses  the  singularities  in  Eq.  B-29  at  p = -b 
and  p = -a  - (2nnj/Ti)  [poles  of  1 - e ].  At  p = -b,  using  Eq.  B-30 

gives 


P=-b  [1  -e_Tl^a'b^][l  -e_T?^S+b^] 


(B-36) 


The  Res  at  p = -a  - (2«nJ/Ti)  is 


,1  1 

'|p.-a-(2mJA,)  * T^(b  — a - (annjAi )[  l — ®+ft+^*njyTi ) j 


n = 0,  ±1,  ±2,  ... 


Since 


, _ e-T2[s+a+(2imJAl)]  = , _ e“T2(a+a)  e~2«njT2Al  = , _ e-T2(s+a) 


if  Tp/T^  is  an  integer, 


-T2(s+a)  T 

-e  c n=  -® 


1 (b  - a)  — 2«nj  ] _e-T2vs+a)  1 _e-Ti(b-a) 


The  final  result  is 


i^2  _ 


l/[1-e-T'(a-b>]  ,/[l-e-T'(b-“)] 

1 - e-T2(s*t>!  * , _ 


Equation  B-^0  highlights  the  fact  that  it  is  not  really  necessary  to 
use  the  phantom  sampler  approach  when  T2/T1  is  an  integer;  a direct  evalua 
tion  via  contour  integration  is  feasible. 
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APPENDIX  C 

ILLUSTRATIVE  EXAMPLE  — METHODS  FOR 
FINDING  THE  CONTINUOUS  RESPONSES 
OF  DISCRETELY  EXCITED  SYSTEMS 


INTRODUCTION 

An  example  is  used  to  demonstrate  various  techniques  for  finding  both 
the  discrete  and  continuous  response  of  a discretely  excited  system.  The 
example  is  a scalar  one  mid  hence  will  not  make  the  dimensionality  con- 
siderations associated  with,  for  example,  the  state  transition  approach 
overly  clear.  However,  our  objective  is  to  expose  the  reader  (or  remind 
him)  of  the  large  number  of  methods  available  for  computing  transient 
responses.  All  the  methods  demonstrated  can  be  applied  equally  as  well 
when  ti\e  system  description  is  in  a matrix  format. 

ILLUSTRATIVE  EXAMPLE  MODEL 

The  model  for  the  illustrative  example  is  given  in  Fig.  C-i.  For 
numerical  convenience  let 


r(t) 


-bt 


e bT ' = 0.^  , T = 0.^ 


and 


so  that 


-aT 

a 

a 

b 


= 0.2* 


?.77Pr'0fl7  PP 
1. 386091*  36 1 


Although  the  tables  used  in  the  "lookup"  process  will  be  in  terms  of  e 
etc.,  for  convenience  of  presentation,  we  ask.  the  reader  to  keep  in  mind 
that  e bnT  = (e~b^)n  = 0.  • 


■bnT 


*7 


■■■HI 


I I 
- . 11 


R 

l-e"*T 

0 

X 

T 

s 

s + a 

Figure  C-l.  Illustrative  Example 


The  fundamental  equations  describing  the  responses  are 


First  we  enumerate  six  methods  for  finding  the  discrete  time  response,  X(nT). 

DISCRETE  TIME  RESPONSE 

Method  I,  Partial  Fraction  Expansion 

This  is  an  "analytical"  method,  in  that  the  response  is  given  in  closed 
form  for  all  values  of  n.  The  technique  is  basically  partial  fraction  expan- 
sion coupled  with  table  lookup. 

Expand  Eq.  C-<?  in  partial  fractions: 

. T'z 

(z  - .2*)l z ~ .5) 

A = 

B = 

X(z)  = 


Az  Bz 

+ 


z - .5  z - .2? 

•Is 


z - .25 


• 75 


z=.5 


(C-3) 


z — .5 


-3 


z=.25 

3z  3z 


z - .5  z - .2^ 


IQ'- 


Via  table  lookup: 


X(nT)  = 3(.5)n  - 3(.25)n  = 3[-5n-.25n] 


(C-4) 


Method  II,  The  Inversion  Integral 

The  inversion  integral  approach  is  essentially  a residue  theory  method 
which  eliminates  the  need  for  table  lookup.  For  all  intents  and  purposes, 
this  method  is  identical  to  Method  1 in  that  the  solution  is  defined  for 
arbitrary  values  of  n.  Using  the  relationship: 


f(nT)  = 


1 

2nd 


J" z11”1  f(z)  dz  = yi  Residues 


(C-5) 


we  obtain 


i-  f 7 — dz  = + 

nj  J (z  - .25) ( z - -5)  .5  - -25  .2*-  - .5 


or 


X(nT)  = 3r.5n  - .25n] 


(C-6) 


Method  III,  Continued  Fraction  Expansion 

In  this  approach,  write  X(z)  as  a ratio  of  polynomials  in  z and  divide 
the  denominator  into  the  numerator,  interpreting  the  negative  powers  of  z 
as  suitably  shifted  delta  functions  in  time.  To  demonstrate: 


x(2)  ■ r--".#*- . 


25) 


z-  - . 75z  + .125 


= . 75z  1 + .5625z  ? 4 .328125Z  5 4 . 17578l25z-Jt  4 


Therefore 


X ( t ) = .758(t-T)  ^ . 562^ t-2T)  4 . 328l255(t -3T)  4 


(C-7) 
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Mithod  IV,  Recursion  [R(s)  Given  Analytically] 

Suppose  a closed  form  is  given  for  R(z),  as  in  Eq.  C-2.  Then  we  may 
write 


(z2  - .75*  + .125)X1(z)  « .75a 

or 

X(nT)  = .75X[n  - 1 )T]  - .125X[n  - 2)T]  + -756(n  - l) 

Mithod  V,  Recursion  [R(s)  in  Syabolic  Fora] 

Suppose  we  retain  the  input  as  R(z)  without  taking  e(-t)t,)T  . Then, 
again  from  Eq.  C-2,  one  may  write 

or 

X[n  + 1)T]  = .25X(nT)  + -75R(nT) 

Method  VI,  State  Transition 


(C-8) 


(C-9) 


(C-10) 


State  transition  can  be  used  to  find  the  total  response.  Using  state 
transition  concepts  we  may  write  (refer  to  Fig.  C-2): 


-a  a If  x 


o o J L X2 


+ R(nT)  - Ax  + BR 


Figure  C-2.  Illustrative  Example  Relabeled 


1Q7 


-ax  , -ax 
e 1 -e 


so  that 


<p(  x)  * <K  '[  IS  - A]  = 


i(nT  + x)  = <p(  x)  x(nT)  + f 

Jo 


d?K(nT) 


(C-ll) 


gives,  tor  the  X)  component.  (Note:  x.-(nT)  » R(nT)  ) 


. (nT  + x)  - (1  - e aT)  R(nT)  + e nTx,(nT)  , 0 < x < T 


Let  x = T and  obtain  the  difference  equation 


Cl[(n  + 1 )T]  = (1  - e aT)  R(nT)  + e aTx1(nT) 


x,[(nt1)T]  - . 7‘'R(nT)  + .2*  xj(nT) 

Equation  C-15  is,  of  course,  in  agreement  with  Eq.  C-10. 

SUMMARY  OF  DISCRETE  RESPONSE  METHODS 


Method  I,  Partial  Fraction  Expansion 


(C-12) 


(C-l?) 


. 7r'X 

’ z - . ?V)(z  - 


5z ^z 

z - .*>  z - .25 


Therefore 


X(nT)  = 3[(.^)n  - (.2‘S)"] 


Mathod  II,  Invaraion  Integral 


x(ut)  . ~ 

^n.1  J l z - - 


TT  di 


■ 3( ,M( .M1'-1  - .?•')( ! 


^ ,*.n  _ ,">*n] 


Method  III,  Continuad  Fraction  Expanaion 


\(z) 


.7*' 

'z  - , 


Mat hod  IV.  Recuraion  fR(*^  01  van  Analytically] 

(■&'-  . !'z  ♦ . t?*  'Xu'  • . ••  u 


There  fore 


X(nT'»  - . -‘Xf(n-1)T]  - . PWf  (n -?)  ] + .7*^(n-n 


Mathod  V.  Racuraion  [R(a)  in  Symbolic  Fora] 


x(*'  - -T— j-",  X K(a' 

There  Tore 

Xf  (n  ♦ 1 )T]  - . 7*'R ( nT ) ♦ ..^X(nT) 

Mathod  VI,  Stata  Tranaition 

Kvnlunte  the  state  transition  equation  at  t T to  get  dlserete  te*vusi 
t ion  equal  Ion: 

X(nT  t)  - (l  - a aT)  R(nT)  ♦ e aTX(nT)  , 0 i t i T 

• ( 1 - e~f*T ) R( nT ) ♦ e" aT  X(nT') 

- . T*  R ( nT ) ♦ .P‘-X(nT> 
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xC(n  ♦ 1 Vr] 


Table  of  Tims  Response 


Using  any  of  the  six  approaches  we  arrive  at  the  summary  table  given 
in  Table  C-1. 


TABLE  C-1 . DISCRETE  TRANSIENT  RESPONSE 


n 

X(nT) 

n 

X(nT) 

0 

0 

6 

.046 142578 

1 

• 75 

7 

.025254395 

2 

.*>625 

8 

.011672974 

3 

.328125 

9 

.00*847931 

4 

. 17*7812* 

10 

.002926826 

c. 

.090820313 

COHTXlfUOUS  RESPONSE 

We  will  illustrate  two  basic  methods  for  finding  the  continuous 
response  of  a discretely  excited  system: 

I.  State  Transition 

II.  Continued  Fraction  Expansion 

A third  method,  the  T/N  approach,  is  treated  in  the  text. 

Msthod  I,  Stats  Transition 

The  state  transition  equation  is 

X(nT  4 t)  « (1 -e“aT)  R(nT)  4 e_aTX(nT)  , 0 £ t * T (C-l4) 

This  is,  in  a strict,  sense,  a recursion  equation.  That  is,  we  must 
know  R(nT)  and  X(nT)  in  order  to  compute  X(nT  4 t)  where  0 < r < T.  How- 
ever. IT  R(nT)  is  given  analytically  (as  in  this  example),  then  one  may 
use  the  "closed  form"  discrete  results  in  conjunction  with  Eq.  C-14  to 


Jump  ahead"  in  time  and  compute  the  continuous  response  for  nTSt< 
To  illustrate,  a "normal"  application  of  Eq.  C-14,  starting  at  t - 0 
give  the  results  shown  in  Table  c-2. 

TABI£  C-2.  CONTINUOUS  RESPONSE 


t 

X(nT + t) 

0 

0 

. 1 

.24214177 

.2 

• 4256?>0823 

.3 

• !>64 724718 

.4 

.670123022 

.5 

.6 

.689464*71 

.7 

. 64  3*87294 

.8 

.608818821 

• 9 

. *8246924* 

1.0 

.*625 

Now,  if  desired,  we  could  skip  ahead  to,  for  example,  t.  * 4,5  since 


t ■ nT  = 9(*5)  * 4.5  sec 

R(9T)  - (,^)9 

X(9T)  =>  3[.?9-.2^] 

Therefore 


X (9T  ♦ t)  - 0-e~aT)(.5)9  + e_aTr3(.5)9  - (,25)9] 

= (1 -e~aT)(. 0019^31??)  + e“aT(.  00^84  7930 


This  gives  Table  C-3. 


(n  1 )T. 
, would 


(C-1f>) 


001 


TABLE  C-3 

CONTINUOUS  RESPONSE  — "JUMPING  AHEAD" 


t 

X(nT + t) 

1* . ^ 

.00‘'8it793l 

4.6 

.OOitdOl*83(' 

^ • 7 

.00it19010it 

4.8 

.0036481438 

it. 9 

.0030379 ’ip 

■'.0 

. 00 ?9?68p6 

Method  II,  Continued  Fraction  Expanaion 

The  continued  fraction  expansion  method  is  easily  extended  to  cover 
the  continuous  response  of  a discretely  excited  system.  Consider  Eq.  C- 1 : 


X = 


s + a s 


-sT 

e RT 


which  has  the  form 


X = F1(s)Fp(z) 


(C-16) 


( C - 1 7 > 


Expand  Fp(z)  in  a continued  fraction  expansion; 


“1  — ,o 

X = F (s)[aQ  + + a^z  ' + •••] 


(C-l8) 


There  fore 


x(t)  = aQf^t)  u(t)  + a 1 f 1 ( t “ T ) u(  t - T ) + apf^t-OT)  u(t.-PT)  ...  (C-19) 
For  our  example 


X = 


— .''Z  - .2*'Z  ‘ 


-(.^)nz-n  - 


POP 


There  fore 

x(t)  = (1  - e_at)  u(t)  - - e_a(t_T)]  u(t  - T) 

- .25[1  - e“a(t-?T)]  u(t  - 2T)  - ••• 


(C-20) 


Equation  C-20  gives  the  same  numerical  values  shown  in  Tables  C-2  and  C-3- 

SUMMARY 

There  are  a variety  of  methods  for  finding  the  continuous  response  of  a 
discretely  excited  linear  time  invariant  system.  They  can  be  roughly  cata- 
loged as: 

• State  transition. 

• Continued  fraction  expansion. 

• Higher  rate  sampling  in  conjunction  with  either: 

— Partial  fraction  expansions. 

— Recursion  equations. 
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APPENDIX  D 


PARTIAL  FRACTION  COEFFICIENTS  FOR 
SECOND -ORDER  REAL  P0I28 


Express. ons  for  finding  the  partial  fraction  coefficients  for  real  second- 
order  poles,  in  terms  of  the  numerator  and  denominator  polynomials  (and  their 
higher  derivatives),  are  derived.  If 


F(s) 


N_  = A B 

P (s  + a)^  (s  + ai 


(D-l) 


the  A coefficient  is  found  by  multiplying  each  side  of  the  equation  by  (s  + a)? 
and  applying  L'Hopital's  rule  twice.  Thus 


A 


-a 


0 

0 


(r-?) 


is  indeterminate  if  (s  -t- a)^  is  not  explicitly  factored  out  of  the  denomina- 
tor and  canceled  with  the  numerator  factor  of  (s  + a)?.  Applying  L'Hopital's 
rule  once  gives 


P(  s + a)N  > 
D' 


0 


(D-3) 


Applying  the  rule  only  once  leaves  an  indeterminate  form,  since  each  numera- 
tor term  contains  free  (s  + a)  factors  and  the  denominator,  D'.  now  contains 
only  a simple  pole  at  s = —a. 

Using  L'Hopital's  rule  a second  times  gives 

A = P(s  ♦ a)N’  + PN  ♦ (s  + a)rN"  + g(s  + a)N' 

D" 


Therefore 


will  give  the  correct  value  for  the  A coefficient,  since  D"  is  now  free 
of  the  double  singularity  at  s = —a.  This  verifies  the  result  stated  in 
Section  II. 

To  find  B,  subtract  A/(s  + a)^  from  each  side  of  Eq.  U-1  and  use 
L'Hopital's  rule  three  times: 

_N A _ N(s  + a)g  — AB 

^ (s  + a)2  D(  s + a)2 

Multiply  by  (s  + a)  and  evaluate  the  result  for 


s + a 


+ • • 


s * -a: 


(D-5) 


N(s  + a)2  - AD  _ 0_ 

D(s + a)  ~ 0 


(D-6) 


A 

Since  D conteins  (s  + a)  , both  numerator  terms  are  zero  (as  is  the  denomi- 
nator) when  s = -a.  Using  L'Hopital's  rule  for  the  first  time  gives 

2(s  + a)N  + N’(s  + a)2  - AD' 

D' (s  + a)  + D 


8 = -a 


_0_ 

0 


(D-7) 


since  every  term  in  the  numerator  still  contains  a free  (s  + a),  as  does  each 
term  of  the  denominator.  Use  the  rule  for  the  second  time: 


U(s  + a)N'  + (s+a)gN"  + (2N  — AD") 
DM(s  + a)  + 2D’ 


s = -a 


_0_ 

0 


(D-8) 


Again,  each  term  of  both  the  numerator  and  denominator  is  zero  when  s = -a 
by  virtue  of  containing  a free  (s+a)  or  the  definition  of  A itself.  That 
is,  the  last  term  in  the  numerator  of  Eq.  D-8  is  identically  zero  at  s = — a 
because : 


2N  - AD" 


Is  = -a 


PN 

2N  - P D" 


= 0 


(D-9) 
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APPENDIX  E 
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h 

* 

H 


j 


SINE  WAVE  INPUT  OF  ARBITRARY  PHASE 


INTRODUCTION 


The  frequency  response  of  a discretely  excited  system  was  developed, 
in  Section  V,  tor  a "zero  phase"  sine  wave  input.  However,  it  is  a 
straightforward  exercise  to  extend  the  development  to  encompass  a sine 
wave  input  of  arbitrary  phase.  This  will  be  done  first  for  the  finite  N 

case  followed  by  the  limiting  process  (N  — «)  which  gives  the  continuous 
results. 

FINITE  N CASE 


In  Fig.  E-l,  let.  the  input  be 

r(t)  = k?  sin  bt.  + kp  cos  bt 
with  the  s-plane  representation 


R(s) 


k^b  kps 

s-’  + b2  s?  + b^ 


Using  z - eSI/N  gives  the  z-plane  representation 

rT  _ k ]ZN  sin  bT  + kpzN(zN  - cos  bT) 
z'  N - 2 cos  bTzN  + 1 


R / 

M 

C / CT/N 

T 

G 

T/N 

(E-l ) 


(E-?) 


(E-3) 


Figure  E-l.  Open-Lioop  System 
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J ] 

I 1 


so  that 


T/n  ^l*1*  8in  W ♦ k2*W(*,f  - cos  bT) 


*2N  - 2 cos  bT*1*  ♦ 1 


An*  sin  <%(T/N)  4-  Bn*[«  - cos  u^T/H)] 
z2  — 2 coa  ^(t/k)*  + 1 


Terms  due  to 
+ the  nodes  of 


(gm)T/n 


In  Eq.  E-1*, 


n = 0,  1,  2,  . . N - 1 


so  that  the  "steady  state"  output  waveform  is  given  by 


(N-1 

X)  ^An 

n=o 


sin  u^t  + Bn  cos  a^t] 


We  can  now  follow  the  same  procedure  as  in  Section  V in  order  to  define 
the  partial  fraction  coefficients  An  and  Bn  in  terms  of  the  system  GM. 
First,  multiply  each  side  of  Eq.  E-4  by  (z2  - 2 cos  colt(T/N)z  + l)  and 
evaluate  for  z = I^u^T/N).  as  in  Section  V,  the  only  surviving  term  on 
the  right-hand  side  of  Ea.  E-U  occurs  when  n = k.  Therefore,  we  obtain 
(letting  k revert  to  n): 

A , iB  - (rn) T/N -N-1  fk1  sin  bT  + k2(^w  - cos  bT)l 
An  + jBn  - (GM)  z — a>n(T/N) jx 


z2  - 2 cos  cd^T/iOz  + 1 
Z2N  _ 2 cos  bTz^  4 1 


z=  12K(t/N) 


APPENDIX  F 


NON -SYNCHRONOUS  SAMPLING 


Non-synchronous  sampling  is  a basic  tool  that  can  be  used  to  model 
distributed  computation  architecture,  data  skewness  in  the  A/D  and  D/A 
conversion  processes  as  well  as  the  internal  computational  delay  of  the 
digital  computer.  Ry  definition,  non-synchronous  sampling  occurs  when  all 
the  systems'  sampling  operations  are  repeated  at  the  same  rate  but  occur 
at  different  instants  of  time.  This  is  depicted  in  Pig.  F-1. 


Figure  F-1.  A Set  of  Non -Synchronously  Sampled  Signals 

In  Fig.  F-1  both  continuous  signals,  x1  and  Xp,  are  sampled  at. 
l/T  samples/second,  but  the  Xp  sampler  is  "out  of  sync"  with  the  x,  sampler 
by  T0  seconds.  The  sampling  operation  for  x'^  is  depicted  symbolically  in 
Fig.  F-2a  and  for  xj1?  in  Fig.  F-2b  (*  notation  on  Xp  indicates  that  an  "uncon- 
ventional" sampling  operation  is  being  carried  out). 

In  Fig.  F-Pb  we  indicate  that  the  non-synchronous  sampler  can  be  modeled 
by  a synchronous  sampler  if  we  precede  the  sampler  with  the  operation  W# 
followed  by  the  operator  W.  That  is, 


Figure  F-?.  Sampling  Notation 


x*  = W(W#x?)T  (F- 1 ) 

where 

W ^ e sTc>  , W # = esT°  (f -2) 

II'  we  proceed  according  to  Eq.  F-?,  one  advances  xp  by  To  seconds, 
samples  at  the  l/T  rate,  and  then  delays  (W#Xo)^  by  T0  seconds  to  obtain 
the  desired  time  sequence.  This  is  demonstrated  in  Fig.  F-3.  Note  that 
the  non-synehronous  sampling  operation  on  xp  can  be  modeled  in  terms  of  a 
scalar  factor;  thus  the  dimension  of  the  equivalent  single-rate  sampled 
signal,  (W#Xp)  , is  not  increased. 

This  model  readily  extends  to  the  case  vnere  x is  a vector. 


Figure  F-3.  Advance,  Sample,  Delay 
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APPENDIX  0 


CONVERSION  FROM  ANALOG  TO  DISCRETE 
STATE  EQUATIONS* 


The  development  of  Eq.  193  in  Section  III  is  outlined  using  a combined 
frequency  and  time  domain  approach  (page  648,  Ref.  20;  page  1h4,  Ref.  21). 
This  method  assumes  that  the  discrete  input  to  the  continuous  element  is 
reconstructed  using  a zero-order  hold  (ZOH). 

Consider  the  simplified  aircraft  model  of  Eq.  G-1: 


The  continuous  state  equations  exist  in  the  form; 

x = A x(t)  + B u(t) 
y = CT  x(t)  + D u(t) 


The  discrete  state  equations  are  given  by  (assuming  a ZOH): 


*k+1  = 

q'(T)Xk  + h(T)  U(k) 

Cp(T)  = 

<p(t)| 

| t.=T 

rT 

h(T)  = 

1 q'(a)B  da 

(G- 1 ) 


(G-2’1 

(G-3) 


(G-4 ) 


( G 


(G-6) 


* 

I 


*This  appendix  vas  prepared  by  Capt.  Dennis  G.  J.  Didaleusky. 
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The  continuous  state  transition  matrix,  qi(t),  is  an  nxn  matrix  describing 
the  fundamental  characteristics  of  the  continuous  aircraft  model  (Eq.  G-l) 
and  is  not  influenced  by  the  inputs  to  the  system.  The  state  transition 
matrix  is  obtained  via  the  s-doraain  by 


<p(t)  = i.  ’i [si  - a]-1  i 

For  Eq.  G-1, 


<T 


(t)  = £ 


( (s  + 2)?  + (6)2  ) 


s + 3 -37 

1 s + 1 


Rearranging  produces 


p(t)  = <£ 


-1 


(s  + 3)  + (l/6)(6)  \ -( 37/6) (6) 

( l/6)(6)  ! (s  + 2)  -(i/6)(6) 


I! 


(s  4 2)2  + (6)2 


(G-7) 


(G-8) 


(G-9) 


Equation  G-9  is  now  in  a form  that  enables  us  to  use  Table  1 in  Section  II 
directly.  Taking  the  inverse  Laplace  transform  of  Eq.  G-9  produces: 


— Pt, 

e cos  6t  + (l/6)e  ' sin  6t  1 

— (37/6)e  sin  6t 

( 1 /6)e  ^ sin  6t 

e ‘ * cos  6t  - (l/6)e  * sin  6t 

q'(t)  = 


For  a smpling  period  of  T = 0.1,  the  transition  matrix  becomes 


(G-10) 


t(t) 


T=0. 1 


0. 7t'277(rO10  J -2.8^078930^ 

1— 

0.0770)t83(0  j 0.^9 8079290 

? 13 


(G-n) 


Integrating  produces  (T  = 0.1): 


-4.490''7<'''94  1 -2.8*0789504 

h(T)  = ! (G-15) 

-0.212719^58  1 -0.401 520710 

v(t)  and  h (T)  can  be  obtained  alternatively  by  making  numerical  e\  \luations 

of  the  Taylor  series  expansions  l’or  eA‘  and  i'or  / 1 eA^“a^  p <ia  , The  dis- 

o 

Crete  state  equations  can  now  be  directly  written  using  Kqs.  0-4.  0-11,  and 
G-15. 

XlRO  0.7*  277(010  -2. 8*078950*  X1r 

Xpk+ , = 0.077048560  0.^98679290 


APPENDIX  H 


A PARTICULAR  DESIGN  PROCEDURE 


Several  closed-loop  designs  were  presented  in  the  text  without  detailed 
comment  on  the  manner  in  which  the  compensation  networks  were  synthesized. 

For  the  most  part,  a pole  assignment  technique  was  adequate  to  meet  the 
design  requirements  posed  in  these  simple  illustrative  examples.  The  use 
of  the  pole  assignment  procedure  will  be  outlined  using  the  "motivating" 
example  of  Section  IV-G. 

From  Fig.  36,  we  may  write  the  output  equation  as 

C = (GM)[I  + G^  g|  (GM)T]_1  G^R1  (H-l) 

One  design  objective  was  to  force  a closed-loop  response,  at  the  sampling 
instant,  dictated  by 


det[ 1 + G|  (GM)t]  oc  z + a (H-2) 

where  a is  determined  by  the  s-plane  pole  at  s = -0.5. 

The  second  requirement,  a unity  steady-state  response  to  a step  input, 
was  forced  by  applying  the  final  value  theorem  to  Eq.  H-1.  These  two  con- 
ditions are  sufficient  to  determine  the  free  parameters  of  G-j  and  Gj. 
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APPENDIX  I 


THE  LIMITING  FORM  FOR  A PARTICULAR  EXAMPLE 


The  steps  in  obtaining  the  limiting  form  given  in  Eq.  509  are  contained 
in  this  appendix.  Specifically,  we  show 


N 


/ 1 ~«T\ 
(±-e"_la 


s ( s + aV 


T/N 


i -«T 

1 -e  a 


UfcDn(T/N)|  Lim 

I N-^o. 


IT  s + a 


SSH 


(i-D 


For  the  example  we  wish  to  evaluate 


(GM) 


T/N 


N 


(1  -e  sT)a 


s(  s + a) 


T/N 


1 1 - e aT,/N  1 - z N 


z - e 


(1-21 


at  z = 14o)n(T/N)  = cos  a^(T/N)  + J sin  a^(T/N)  = e>n(T/N)  . 


(GM) 


T/N 


(1  - e~aT/N)0  - e~^*T) 
N[eJ^(T/N)  _ e~aT/N3[l  _ e-jo>n(T/N)3 


(1-3) 


For  the  sake  of  brevity,  let  1 /N  = u,  so  that  u -►  0 as  1)  -►  »,  and  consider 


F = 


_N 

D 


u 

^jtonTu  _ 


LLt^.1  . . . 

e aTu]M  - e *}a'hTu] 


(I-1*) 


Here  F represents  the  part  of  (GM/N)^^  which  is  "troublesome"  as  u 0. 
Clearly,  F becomes  indeterminate  as  u Oj  therefore,  use  L'Hopital’s 
rule  (for  the  first  time): 


dN/du 

dD/du 


-aTu,  _ -aTu 
1 — e ) + aTue 


(V)  «J“"TU  ♦ «Te-“TU  - T(Jl%  * 


(1-f) 


As  u -»-0,  the  form  remains  indeterminate;  therefore,  use  L'Hopital's 
rule  for  the  second  time: 


d^N/du^ 


aT<TaTU  , aTe~ aTu  + „ T-f  »T1g^~aTul 


d2D/a“2  (Jagr)2.**1"  - (aT)2e-aTu  ♦ !*(**♦  a)2e-Ta(J«n«) 


(J%T)S  - (aT)2  * T2^  + a) 


(1-6) 


T KJ^)2  — + (j(%)2  + gajoijj  +/?] 


2T  jc%[ju^  + a]  joJnTCjojn  + a] 


Finally,  reintroduce  the  (l  - e ^T)  term  to  obtain 


1 - e~Ju*T  a 
J<%T  + a 


, -sT 
1 — e a 

sT  s + a 


s=Mi 


This  verifies,  after  setting  a = 1,  Eq.  309  of  Section  V. 


(1-7) 


(1-8) 
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